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Abstract 

By virtue of the newly developed technique of integration within an ordered product (IWOP) 
of operators, we explore quantum optical version of classical optical transformations such as opti- 
cal Fresnel transform, Hankel transform, fractional Fourier transform, Wigner transform, wavelet 
transform and Fresnel-Hadmard combinatorial transform etc. In this way one may gain benefit 
for developing classical optics theory from the research in quantum optics, or vice-versa. We can 
not only find some new quantum mechanical unitary operators which correspond to the known 
optical transformations, deriving a new theorem for calculating quantum tomogram of density op- 
erators, but also can reveal some new classical optical transformations. For examples, we find the 
generalized Fresnel operator (GFO) to correspond to the generalized Fresnel transform (GFT) in 
classical optics. We derive GFO's normal product form and its canonical coherent state represen- 
tation and find that GFO is the loyal representation of symplectic group multiplication rule. We 
show that GFT is just the transformation matrix element of GFO in the coordinate representation 
such that two successive GFTs is still a GFT. The ABCD rule of the Gaussian beam propagation is 
directly demonstrated in the context of quantum optics. Especially, the introduction of quantum 
mechanical entangled state representations opens up a new area to finding new classical optical 
transformations. The complex wavelet transform and the condition of mother wavelet are studied 
in the context of quantum optics too. Throughout our discussions, the coherent state, the entan- 
gled state representation of the two-mode squeezing operators and the technique of integration 
within an ordered product (IWOP) of operators are fully used. All these confirms Dirac's assertion: 
" ...for a quantum dynamic system that has a classical analogue, unitary transformation in the 
quantum theory is the analogue of contact transformation in the classical theory". 

Keywords: Dirac's symbolic method; IWOP technique; entangled state of continuum variables; 
entangled Fresnel transform; Collins formula; Generalized Fresnel operator; complex wavelet 
transform; complex Wigner transform; complex fractional Fourier transform; symplectic wavelet 
transform; entangled symplectic wavelet transform; Symplectic-dilation mixed wavelet transform; 
fractional Radon transform; new eigenmodes of fractional Fourier transform 
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1 Introduction 



The history of quantum mechanics records that from the very beginning the founders of the quan- 
tum theory realized that there might exist formal connection between classical optics and quantum 
mechanics. For example, Schrodinger considered that classical dynamics of a point particle should 
be the "geometrical optics" approximation of a linear wave equation, in the same way as ray optics 
is a limiting approximation of wave optics; Schrodinger also searched for some quantum mechanical 
state which behaves like a classical 'particle', and this state was later recognized as the coherent state 
IUE1E2], which plays an essential role in quantum optics theory and laser physics; As Dirac wrote 
in his famous book < Principles of Quantum Mechanics > 0J : "• • • for a quantum dynamic system 
that has a classical analogue, unitary transformation in the quantum theory is the analogue 
of contact transformation in the classical theory". According to Dirac, there should exist a formal 
correspondence between quantum optics unitary-transform operators and classical optics transfor- 
mations. Indeed, in the last century physicists also found some rigorous mathematical analogies 
between classical optics and quantum mechanics, i.e. the similarity between the optical Helmholtz 
equation and the time-independent Schrodinger equation; Since 1960s, the advent of a laser and 
the appearance of coherent state theory of radiation field HTJ [2l [5j , quantum optics has experienced 
rapid development and achieved great success in revealing and explaining the quantum mechanical 
features of optical field and non-classical behavior (for instance, Hanbury-Brown-Twiss effect, pho- 
ton antibunching, squeezing, sub-Poissonian photon statistics) of photons in various photon-atom 
interactions J6j . The relationship between classical and quantum coherence has been discussed in 
the book of Mandel and Wolf ||6|| ; The Hermite-Gauss or the Laguerre-Gauss modes of a laser beam 
are described using the bosonic operator algebra by Nienhuis and Allen [7]. In addition, displaced 
light beams refracted by lenses according to the law of geometrical optics, were found to be the 
paraxial optics analog of a coherent state. Besides, phase space correspondence between classical 
optics and quantum mechanics, say for example, the Wigner function theory is inspected in the 
literature flS). 

On the other hand, classical optics, which tackles vast majority of physical-optics experiments and 
is based on Maxwell's equations, has never ceased its own evolving steps, physicists have endeavored 
to develop various optical transforms in light propagation through lens systems and various contin- 
uous media. The two research fields, quantum optics and classical optics, have their own physical 
objects and conceptions. From the point of view of mathematics, classical optics is framed in the 
group transform and associated representations on appropriate function space, while quantum op- 
tics deals with operators and state vectors, and their overlap seems little at first glance. It seems to 
us that if one wants to further relate them to each other, one needs some new theoretical method 
to "bridge" them. For example, what is the quantum mechanical unitary operator corresponding to 
the Fresnel transform in Fourier optics? Is there any so-called Fresnel operator as the image of clas- 
sical generalized Fresnel transform? Since generalized Fresnel transforms are very popularly used 
in optical instrument design and optical propagation through lenses and various media, it is worth 
of studying these transforms in the context of quantum optics theory, especially based on coherent 
state, squeezed state (9\ 11011 and the newly invented entangled state theory Hill 1121 [T3l [I4ll . 

Fortunately, the recently developed technique of integration within an ordered product (IWOP) 
of operators H151 [T6l I17H is of great aid to studying quantum optical version of classical optical 
transformations. Using the IWOP technique one may gain benefit for classical optics from quantum 
optics' research, or vice-versa. 

Our present Review is arranged as follows: in section 2 we briefly recall the classical diffraction 
theory H181 |19L this is preparing for later sections in which we shall show that most frequently 
employed classical optical transforms have their counterparts in quantum optics theory. In sec- 
tion 3 we introduce the IWOP technique and demonstrate that the completeness relation of funda- 
mental quantum mechanical representations can be recast into normally ordered Gaussian opera- 
tor form. Using the IWOP technique we can directly perform the asymmetric ket-bra integration 
A 4 " 1 ^ 2 jTlo dq | <?//■») {l\ i n tne coordinate representation, which leads to the normally ordered single- 
mode squeezing operator, this seems to be a direct way to understanding the squeezing mechanism 
as a mapping from the classical scaling q — > q/fj,. In section 4 with the help of IWOP technique 



and based on the concept of quantum entanglement of Einstein-Podolsky- Rosen [20] we construct 
two mutually conjugate entangled states of continuum variables, \r]) versus |£) , and their deduced 
entangled states (or named correlated-amplitude — number-difference entangled states), they are 
all qualified to make up quantum mechanical representations. It is remarkable that using the IWOP 
technique to performing the asymmetric ket-bra integration fi^ 1 J d 2 r] \r]/p,) (77 1 leads to the two- 
mode normally ordered two-mode squeezing operator, this implies that the two-mode squeezed 
state is simultaneously an entangled state. We point out that the entangled state \rj) also embodies 
entanglement in the aspect of correlative amplitude and the phase. We are also encouraged that 
the overlap between two mutually conjugate deduced entangled states is just the Bessel function — 
the optical Hankel transform kernel H21H , which again shows that the new representations in the 
context of physics theory match beautiful mathematical formalism exactly We then employ the 
deduced entangled states to derive quantum optical version of classical circular harmonic corre- 
lation. Section 5 is devoted to finding a quantum operator which corresponds to the optical Fresnel 
transform, with use of the coherent state representation and by projecting the classical sympletic 
transform z — > sz — rz* (|s| 2 — |r| 2 = 1) in phase space onto the quantum mechanical Hilbert space, 
we are able to recognize which operator is the single-mode Fresnel operator (FO) . It turns out that 
the 1 -dimensional optical Fresnel transform is just the matrix element of the Fresnel operator F 
in the coordinate eigenstates. Besides, the coherent state projection operator representation of FO 
constitutes a loyal realization of symplectic group, which coincides with the fact that two succes- 
sive optical Fresnel transforms make up a new Fresnel transform. Then in Section 6 based on the 
coherent state projection representation of FO, we prove ABCD rule for optical propagation in the 
context of quantum optics. In section 7 the quadratic operator form of FO is also presented and the 
four fundamental optical operators are derived by decomposing the FO. In section 8 we discuss how 
to apply the Fresnel operator to quantum tomography theory, by introducing the Fresnel quadrature 
phase FXF^ = Xp, we point out that Wigner operator's Radon transformation is just the pure state 
projection operator \x) s rs (x\, where \x) s — F\x) and \x) is the position eigenstate, so the prob- 
ability distribution for the Fresnel quadrature phase is the Radon transform of the Wigner function. 
Moreover, the tomogram of quantum state \tp) is just the squared modulus of the wave function 
S} r (x\ ip) . This new relation between quantum tomography and optical Fresnel transform may pro- 
vide experimentalists to figure out new approach for testing tomography. In addition, we propose 
another new theorem for calculating tomogram, i.e., the tomogram of a density operator p is equal 
to the marginal integration of the classical Weyl correspondence function of F^ pF. In section 9 by 
virtue of the coherent state and IWOP method we propose two-mode generalized Fresnel operator 
(GFO), in this case we employ the entangled state representation to relate the 2-mode GFO to clas- 
sical transforms, since the 2-mode GFO is not simply the direct product of two 1-mode GFOs. The 
corresponding quantum optics ABCD rule for two-mode case is also proved. The 2-mode GFO can 
also be expressed in quadratic operators form in entangled way. The relation between optical FT and 
quantum tomography in two-mode case is also revealed. In section 10 we propose a kind of integra- 
tion transformation, J -222.e 2% ( p ~ x ^ q ~ v ^h(p, q) = f (x, y) , which is invertible and obeys Parseval 
theorem. Remarkably, it can convert chirplet function to the kernel of fractional Fourier transform 
(FrFT). This transformation can also serve for solving some operator ordering problems. In section 
11 we employ the entangled state representation to introduce the complex FrFT (CFrFT), which is 
not the direct product of two independent 1 -dimensional FrFT transform. The eigenmodes on CFrFT 
is derived. New eigen-modes for light propagation in graded-index medium and the fractional Han- 
kel transform are presented. The Wigner transform theory is extended to the complex form and its 
relation to CFrFT is shown; The integration transformation in section 10 is also extended to the en- 
tangled case. In section 12 we shall treat the adaption problem of Collins diffraction formula to the 
CFrFT with the use of two-mode (3 parameters) squeezing operator and in the entangled state repre- 
sentation of continuous variables, in so doing the quantum mechanical version of associated theory 
of classical diffraction and classical CFrFT is obtained, which connects classical optics and quantum 
optics in this aspect. In section 13 we introduce a convenient way for constructing the fractional 
Radon transform, the complex fractional Randon transform is also proposed; In sections 14 and 
15 we discuss quantum optical version of classical wavelet transforms (WTs), including how to re- 
cast the condition of mother wavelet into the context of quantum optics; how to introduce complex 



wavelet transform with use of the entangled state representations. Some properties, such as Parseval 
theorem, Inversion formula, and orthogonal property, the relation between WT and Wigner-Husimi 
distribution function are also discussed. In section 16, we generalize the usual wavelet transform 
to symplectic wavelet transformation (SWT) by using the coherent state representation and making 
transformation z — > s(z — k) — r(z* — k*) (|s| 2 — \r\ 2 = 1) in phase space. The relation between 
SWT and optical Fresnel transformation is revealed. Then the SWT is extended to the entangled 
case by mapping the classical mixed transformation (z, z') — > (sz + rz'*,sz' + rz*) in 2-mode coher- 
ent state | z, z') representation. At the end of this section, we introduce a new symplectic-dilation 
mixed WT by employing a new entangled-coherent state representation \a,x). The corresponding 
classical optical transform is also presented. In the last section, we introduce the Fresnel-Hadamard 
combinatorial operator by virtue of the IWOP technique and \a, x). This unitary operator plays the 
role of both Fresnel transformation for mode ai ^ 2 and Hadamard transformation for mode ai ^ 2 , 
respectively, and the two transformations are combinatorial. All these sections are used to prove the 
existence of a one-to-one correspondence between quantum optical operators that transform state 
vectors in Hilbert space and the classical optical transforms that change the distribution of optical 
field. 



2 Some typical classical optical transformations 

Here we briefly review some typical optical transforms based on light diffraction theory. These 
transformations, as one can see in later sections, are just the correspondence of some representation 
transformations between certain quantum mechanical states of which some are newly constructed. 

It was Huygens who gave a first illustrative explanation to wave theory by proposing every point 
in the propagating space as a sub-excitation source of a new sub-wave. An intuitive theory mathe- 
matically supporting Huygens' principle is the scalar diffraction approximation, so named because 
optical fields (electromagnetic fields) actually are vector fields, whereby the theory is valid approx- 
imately. This theory is based on the superposition of the combined radiation field of multiple re- 
emission sources initiated by Huygens. Light diffraction phenomena has played an important role in 
the development of the wave theory of light, and now underlies the Fourier optics and information 
optics. The formulation of a diffraction problem essentially considers an incident free-space wave 
whose propagation is interrupted by an obstacle or mask which changes the phase and/ or amplitude 
of the wave locally by a well determined factor [22]. A more rigorous, but still in the scheme of 
scalar wave, derivation has been given by Kirchhoff who reformulated the diffraction problem as 
a boundary-value problem, which essentially justifies the use of Huygens principle. The Fresnel- 
Kirchhoff (or Rayleigh-Sommerfeld) diffraction formula is practically reduced to the Fresnel integral 
formula in paraxial and far-field approximation H181H9H that reads: 
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where U\ (x\ , y{) is the optical distribution of a 2-dimensional light source and U 2 (x 2 , y 2 ) is its image 
on the observation plane, A is the optical wavelength, k = ^ is the wave number in the vacuum 
and z is the propagation distance. When 
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is satisfied, Eq. Q]) reduces to 



U 2 (x 2 ,y 2 ) = 



exp (ikz) exp [i^j {x\ + y 2 )] 



iXz 

Ui (xi,yi)exp 



.2tt . . 
-1— (xix 2 + yiy 2 ) 
AZ 



dx\dyi 



(3) 



r 



which is named the Fraunhofer diffraction formula. 



The Fresnel integral is closely related to the fractional Fourier transform (FrFT), actually, it has 
been proved that the Fresnel transform can be interpreted as a scaled FrFT with a residual phase 
curvature H23H . The FrFT is a very useful tool in Fourier optics and information optics. This concept 
was firstly introduced in 1980 by Namias H24H but not brought enough attention until FrFT was 
defined physically, based on propagation in quadratic graded-index media (GRIN media) . Mendlovic 
and Ozaktas H251 126H defined the ath FrFT as follows: Let the original function be input from one 
side of quadratic GRIN medium, at z = 0. Then, the light distribution observed at the plane z = z 
corresponds to the a equal to the (zo/L)th fractional Fourier transform of the input fraction, where 
L = (7r/2)(ni/n 2 ) 1 / 2 is a characteristic distance. The FrFT can also be implemented by lenses. 
Another approach for introducing FrFT was made by Lohmann who pointed out the algorithmic 
isomorphism among image rotation, rotation of the Wigner distribution function [291, and fractional 
Fourier transforming H27B . Lohmann proposed the FrFT as the transform performed on a function 
that leads to a rotation with an angle of the associated Wigner distribution function, in this sense, 
the FrFT bridges the gap between classical optics and optical Wigner distribution theory. Recently, 
the FrFT has been paid more and more attention within different contexts of both mathematics and 
physics (24l[25j|26l|27l|28l. The FrFT is defined as 
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We can see that F is the identity operator and F v / 2 is just the Fourier transform. The most important 
property of FRFT is that F a obeys the semigroup property, i.e. two successive FrFTs of order a and 
j3 makes up the FrFT of order a + /3. A more general form describing the light propagation in an 
optical system characterized by the [A, B; C, D] ray transfer matrix is the Collins diffraction integral 
formula 113011 , 
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where AD — BC = 1 if the system is lossless. One can easily find the similarity between Collins 
formula and the FrFT by some scaling transform and relating the [A,B,C,D] matrix to a in the FrFT 

[13111 . Note that M = ( ^ ^ ^ is a ray transfer matrix describing optical systems belonging to the 

unimodular symplectic group. When treating the light propagation in optical elements in near-axis 
approximation, matrices M representing linear transformations are a convenient mathematical tool 
for calculating the fundamental properties of optical systems, which is the origin of the name of 
matrix optics. In cylindrical coordinates the Collins formula is expressed as [30, 
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where we have used the m-order Bessel function 
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When A = 0, ® reduces to the standard Hankel transform (up to a phase factor) 
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The compact form of one-dimensional Collins formula is 
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M is the parameter matrix [A,B,C,D]. Eq. (fl2l) is called generalized Fresnel transform ||33l I 
35l [36ll . In the following sections we will show how we find the quantum optical counterpart for 
those transformations of classical optics. For this purpose in the next chapter we introduce the IWOP 
technique to demonstrate how Dirac's symbolic method can be developed and be applied to quantum 
optics theory. Also, we briefly review some properties of the entangled state Hill 1121 [T3l [T4ll and 
reveal the connection between the mutual transform generated by these entangled states and the 
Hankel transform in classical optics. 



3 The IWOP technique and two mutually conjugate entangled 
states 

3.1 The IWOP technique 

The history of mathematics tells us that whenever there appears a new important mathematical 
symbol, there coexists certain operational rules for it, the quantum mechanical operators in ket-bra 
projective form (the core of Dirac's symbolic method) also need their own operational rules. The 
terminology "symbolic method" was first shown in the preface of Dirac's book <The Principle of 
Quantum Mechanics>: "The symbolic method, which deals directly in an abstract way with the quan- 
tities of fundamental importance- ■ ■, however, seems to go more deeply into the nature of things. It 
enables one to express the physical law in a neat and concise way, and will probably be increasingly 
used in the future as it becomes better understood and its own special mathematics gets devel- 
oped." [4] Then two questions naturally arise: How to better understand the symbolic method? How 
to develop Dirac's symbolic method, especially its mathematics? We noticed that Newton-Leibniz in- 
tegration rule only applies to commuting functions of continuum variables, while operators 
made of Dirac's symbols (ket versus bra, e.g., \q/fi) (q\ of continuous parameter q) in quantum 
mechanics are usually not commutative. Therefore integrations over the operators of type 
| ) ( | (where ket- and bra- state vectors need not to be Hermitian-conjugate to each other) 
can not be directly performed by the Newton-Leibniz rule. Thus we invented an innovative 
technique of integration within an ordered product (IWOP) of operators that made the in- 
tegration of non-commutative operators possible. The core of IWOP technique is to arrange 
non-commutable quantum operators within an ordered product (say, normal ordering) in a 
way that they become commutable, in this sense the gap between q-numbers and c-numbers 
is "narrowed". However, the nature of operators which which are within : : is not changed, 
they are still q-numbers, not c-numbers. After the integration over c-numbers within ordered 



product is performed, we can get rid of the normal ordering symbol after putting the inte- 
gration result in normal ordering. H37ll .The IWOP technique thus bridges this mathematical 
gap between classical mechanics and quantum mechanics, and further reveals the beauty and 
elegance of Dirac's symbolic method and transformation theory. This technique develops sym- 
bolic method significantly, i. e. makes Dirac's representation theory and the transformation 
theory more plentiful, and consequently to be better understood. The beauty and elegance 
of Dirac's symbolic method are further revealed. Various applications of the IWOP technique, 
including constructing the entangled state, developing the nonlinear coherent state theory, 
Wigner function theory, etc. are found; many new unitary operators and operator-identities 
as well as new quantum mechanical representations can be derived too, which are partly 
summarized in the Review Articles [11211 . 

We begin with listing some properties of normal product of operators which means all the bosonic 
creation operators are standing on the left of annihilation operators a in a monomial of a) and a. 

1. The order of Bose operators a and a) within a normally ordered product can be permuted. 
That is to say, even though [a, a^l = 1, we can have : aa) : — : a) a: — a^a, where : : denotes normal 
ordering. 

2. c-numbers can be taken out of the symbol : : as one wishes. 

3. The symbol : : which is within another symbol : : can be deleted. 

4. The vacuum projection operator |0)(0| has the normal product form 

|0)(0| =: e- aV . (13) 

5. A normally ordered product can be integrated or differentiated with respect to a c-number 
provided the integration is convergent. 



3.2 The IWOP technique for deriving normally ordered Gaussian form of the 
completeness relations of fundamental quantum mechanical representa- 
tions 

As an application of IWOP, (in the following, unless particularly mentioned, we take % = oj = m = 1 
for convenience.) Using the Fock representation of the coordinate eigenvector Q\q) — q\q), (Q = 
(a + a^)/V2) 

\q) =^ 1/4 e-^ + ^ at -^|0), (14) 

we perform the integration below 
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Substituting ([131) into ([15]) we see 

ft = f -%e~^>'-^: e -V e"*^ 5 -*. (16) 



Note that on the left of : e a+a : are all creation operators, while on its right are all annihilation 
operators, so the whole integral is in normal ordering, thus using property 1 we have 

S 1 = [°° dQ ; e -4(i+^)+v^(4+ a )-^( a + a+ ) 2 : _ ( 17 ) 



As a commutes with a< within : :, so a) and a can be considered as if they were parameters while 
the integration is performing. Therefore, by setting fx — e x , sechA = 1+ ^ 2 , tanhA = ^ 2 t\ , we are 
able to perform the integration and obtain 
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which is just the single-mode squeezing operator in normal ordering appearing in many references. 
It is worth mentioning that we have not used the SU(1,1) Lie algebra method in the derivation. The 
integration automatically arranges the squeezing operator in normal ordering. Using 
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Eq. (fl8l) becomes 
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This shows the classical dilation q — >• ^ maps into the normally ordered squeezing operator man- 
ifestly. It also exhibits that the fundamental representation theory can be formulated in not so 
abstract way, as we can now directly perform the integral over ket-bra projection operators. More- 
over, the IWOP technique can be employed to perform many complicated integrations for ket-bra 
projection operators. 

There is a deep ditch between quantum mechanical operators (g-numbers) theory and classical 
numbers (c-numbers) theory. The IWOP technique arranges non-commutable operators within an 
ordered product symbol in a way that they become commutable, in this sense the 'ditch' between 
g-numbers and c-numbers is "shoaled". However, the nature of operators are not changed, they 
are still q— numbers, not c-numbers. After the integration over c-numbers within ordered product is 
performed, we can finally get rid of the normal ordering symbol by using dl9D . 

When fi = 1, Eq. ([20]) becomes 

dq\q)(q\= f * : e^^M^)" : 

^:c^) 2 : =1, (21) 
> V 71- 

a real simple Gaussian integration! This immediately leads us to put the completeness relation of 
the momentum representation into the normally ordered Gaussian form 

dp\p)(p\= r ^:e-^- p > 2 : =1, (22) 
where P = (a — a*) /(i^/2), and \p) is the momentum eigenvector P\p) = p \p), 



\p) = 7r 4 exp 



1 o . /r + 1 ]2 



■1 



p 2 + i^f^ipa} -+ —a 



|0) • (23) 



In addition, we should notice that \q) and \p) are related by the Fourier transform (FT), i.e. (p\ q) = 
exp (—iqp) , the integral kernel of the Fraunhofer diffraction formula in 1-dimensional is such 



a FT, so FT in classical optics has its correspondence in quantum mechanical representations' trans- 
form. This enlightens us that in order to find more general analogy between unitary operators in 
quantum optics and transformations in classical optics we should construct new representations for 
quantum optics theory, and these are the bi-partite entangled state and many-particle entangled 
state. These ideal states can be implemented by optical devices and optical network H38H - In the 
following we focus on the bi-partite entangled state. 



3.3 Single-mode Wigner operator 

When we combine (I21D and (I22D we can obtain 

TT" 1 : e -(9-<2) 2 -(f-^) 2 : = A(q,p), (24) 

which is just the normally ordered Wigner operator since its marginal integration gives \q)(q\ and 
\p) (p\ respectively, i.e., 



dqA(q,p) = ^ r :e-^ 2 : = \p)(p\, (25) 
V 71 " 

dpA(q,p) = ^=: e -(^) 2 : =\q)(q\. (26) 



Thus the Wigner function of quantum state p can be calculated as W(q,p) =Tr[pA (q,p)]. On the 
other hand, the Wigner operator ( |24t can be recast into the coherent state representation, 

A (q,p) -> A (a, a*) = [ — \a + z) (a - z\ e az *' a ' z , (27) 



7T 

where \z) is a coherent state. In fact, using the IWOP technique we can obtain 

/(Pz 
: exp{— I jar] + (a + z) a) + (a* — z*) a 
TT 

+ az* — a* z — \a\ 2 } : 

= -: exp{-2(a-a)(a t -a*)} : , (28) 

7T 

which is the same as d24D . 



3.4 Entangled state \rf) and its Fourier transform in complex form 



The concept of quantum entanglement was first employed by Einstein, Rosen and Poldosky (EPR) 
to challenge that quantum mechanics is incomplete when they observed that two particles' relative 
position Qi — Q 2 and the total momentum Pi + P 2 are commutable. Hinted by EPR, the bipartite 
entangled state |»7) is introduced as [ ]39l I40H 



\v) = exp 



-\r]\ 2 +r)a\ -r?*4 + 44 



|00). 



(29) 



\r) = rji + ±772) is the common eigenstate of relative coordinate Qi — Q 2 and the total momentum 

Pi + Po 

(Qi-Q2)\r)) = V2r h \ V ), (Pi+P 2 )\ V )=V2 V2 \ V ), (30) 

where Q { = (cij + a^)/y/2, P 3 = (a,j — ap/(iy / 2), j = 1,2. Using the IWOP technique, we can 
immediately prove that \r/) possesses the completeness relation 

/ ^^(^ / ^2 :e -[u*-(aJ-«)][i,-(.i-a2)] : = 1, d 2 V = d m d V 2, (31) 



and orthonormal relation 

(v\v'}=^S(r ll -r ) ' 1 )6(ri 2 -r,' 2 ). (32) 
The Schmidt decomposition of \rj) is 



\ 7] ) =e -"i2m dx\q) 1 ® q-V2ft) e l ^ x , (33) 
The 1 77) state can also be Schmidt-decomposed in momentum eigenvector space as 

/oo 
dp|p>i® ^2?72-p) e- 1 ^"^. (34) 

The 1 77) is physically appealing in quantum optics theory, because the two-mode squeezing operator 
has its natural representation on basis ||40ll 



- - Wjl jj) (77 1 = e a l a 2 tanh A e ( a l Q i+ a 2 a 2 + l) Inscch A g -aia2 tanh A ^ _ g A p^-j 
7T/U 

The proof of (I35D is proceeded by virtue of the IWOP technique 



TT/i ,/ [ 2 \ [L- ) J 

* I a 2\tt t t 1 

+ V ai + a[a 2 + aia 2 — a[ai — a 2 a 2 > 



2/i J /i 2 / a\ \ I a 2 \ ( t 

V :eXP \^lM 

t t , , , , t , t 



tanh Ag(a J ai -f-a^ a2 + l) In scch aia2 tanh A _ (36) 



so the necessity of introducing \rj) into quantum optics is clear. S2 squeezes |ry) in the manifest way 

S 2 \v) = -\v/tA, M = e A , (37) 
M 

and the two-mode squeezed state itself is an entangled state which entangles the idle mode and 
signal mode as an outcome of a parametric-down conversion process H41H . 
We can also introduce the conjugate state of \rj) 114211 , 



10 = exp 



■l\s\ 2 +s4+?4-a\4 



|00), £ = £i+t£ 2> (38) 



which obeys the eigen-equations 

(Qi + Qa) 10 = V26 10 , (-Pi - Pi) 10 = V^ a 10 ■ (39) 

Because [(Qi — Q 2 ) , (-Pi — P%)\ — 2i, so we name the conjugacy between |£) and \q) . The complete- 
ness and orthonormal relations of |£) are 

/ d Ll\£ ) ^\= [ft :e -[e-(4+»2)][t-(»t+4)] : =1 (40) 

(£1 O = tt5(6 - - d 2 £ = dfcdfc, (41) 

respectively. |t;) and |£) can be related to each other by 

m = Up(^^), (42) 



2 V 2 

since £*?y — £77* is a pure imaginary number, Eq. Q42D is the Fourier transform kernel in complex 
form (or named entangled Fourier transform, this concept should also be extended to multipartite 
entangled states.) It will be shown in later sections that departing from entangled states \rj) and |£) 
and the generalized Fresnel operator a new entangled Fresnel transforms in classical optics can be 
found. 



1 o 



3.5 Two-mode Wigner operator in the \rf) representation 

Combining J31D and J40D we can construct the following operator 

_ . p-W"-( a l- a 2)}\<?-(ai-al)]-[y*-(al+a 2 )}h-{a 1 +al)] . 

= A(a,Q*)®A(ftf) = A (cr, 7) , (43) 

where 

a = a-0*, 7 = a + 13*. (44) 

Eq. A43D is just equal to the direct product of two single-mode Wigner operators. It is convenient to 
express the Wigner operator in the I77) representation as H43H 



A ( n, , ) = / ^ a rj) (a + V \ e^'-n* 7. (45) 



function. For example, noticing (?y| 00) = exp{— |7y| 2 /2}, the two-mode squeezed states' Wigner 



For two-mode correlated system, it prefers to using A (cr, 7) to calculate quantum states' Wigner 
function. I 
function is 

(00|4(/i)A(cr, 7 )5 2 (M)|00) 



n 2 exp -n 2 \<r\ 2 - |t| 2 / A* 2 



(46) 



4 Two deduced entangled state representations and Hankel trans- 
form 

4.1 Deduced entangled states 

Starting from the entangled state \rj = re 10 } and introducing an integer m, we can deduce new states 
11441. 

1 r 2 * 

\m, r) = — d8\ v = re w ) e - im8 , (47) 
2tt Jo 

which is worth of paying attention because when we operate the number-difference operator, 



D = a\ai - a\a 2 (48) 



on I rj) , using Eq.fl29l) we see 



D \ V ) = {r,a[ + V*<h) \v) = \v) ,V = \v\e ie , (49) 

so the number-difference operator corresponds to a differential operation i-^L in the (rj\ representation, 
this is a remarkable property of (r]\ . It then follows 







2tt 



D\m,r) = I —e~ lme [~i^-\ri = re w )\^m\m 1 r). (50) 



On the other hand, by defining 

K={ ai -a\){a\-a 2 ), (51) 
we see [D, K] = 0, and |m, r) is its eigenstate, 

K \m, r) — r 2 \m, r) , (52) 
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where K is named correlated-amplitude operator since K\rj) — \-q\ 2 \rj) . Thus we name \m,r) 
correlated-amplitude — number-difference entangled states. It is not difficult to prove completeness 
and orthonormal property of |m, r), 

x - r-OO 

V / d(r 2 ) \m,r) (m,r\ = 1, (53) 



ra— — oo 



(m,r| m',r') = 6 m m >— S(r - r') . (54) 
2r 



On the other hand, from |£) we can derive another state 



(55) 



which satisfies 

D 10 = (ok - 40) 10 = -i^ |£ = • (56) 
So D in (£ = r'e l¥> | representation is equal to i-g-. Consequently, 

D |s ' r ' ] = r ^ e ~ tse ir^o \* = rV "0 = s |s ' r ' } ■ (57) 



Note 



D, (a\ + a 2 )(ai + a\) 



and 

(4 + a 2 ){a x + 4) |s, r') = r' 2 \s, r') . (58) 



\s,r') is qualified to be a new representation since 



/.oo -i 

5] / d(r' 2 )| s ,r')( s ,r'| = l, ( s y | s ', r ") = 5 SiS , £ (/ - r") • 



(59) 



4.2 Hankel transform between two deduced entangled state representations 

Since |0 and \rf) are mutual conjugate, \s, r') is the conjugate state of |m, r). From the definition of 
\m, r) and |s, r') and ( I42t we calculate the overlap 112111 

( s , r '| g , r ) = / (£ = r 'e iv \ / d0 |r? = re w ) e - m9 

4tt Jo Jo 



^ /-27T /.27T 



e 



isip—im9 



exp [irr sin (0 — ip)] d8dip 



o Jo 

i p2ir p2-rr oo 

= — 2 / / e lslp ~ lme V J; (r/) e^ 9- ^ 
8?r -/o ;=-oo 

1 oo 1 

= 2 E s i,m5i, s Jt {rr 1 ) = -5 s , m J s (rr') , (60) 

(= — OO 

where we have identified the generating function of the s-order Bessel function J\ , 

oo 

e ixsint = 2 J ; (x)e lit , (61) 



/ — — oo 

and 



fc=0 



Eq. (1601) is remarkable, because J s [rr 1 ) is just the integral kernel of Hankel transform. In fact, if we 
define 

(m,r\ g) =g(m,r), (s,r'| g) = Q (s,r') , (63) 
and use J53D as well as ( 1601) , we obtain 



G(s,r')= ^2 / d(r 2 ) {s,r'\m,r) {m,r\ g) 

d(r 2 ) J s (rr') ff ( S ,r) = H[.g(s,0] 



171= — OO 
f'OC 



(64) 

which is just the Hankel transform of g (to, r) (or it can be regarded as a simplified form of the 
Collins formula in cylindrical coordinate, see (110D ). The inverse transform of ( 1641 ) is 

00 POO 

g(m,r) = (m,x\ ^ / d (r' 2 ) |s, r') (s, r'\ g) 

S =-co J ° 

d (r' 2 ) J g (rr') Q (m, r') = H" 1 (m, r')] . (65) 

Now we know that the quantum optical image of classical Hankel transform just corresponds to the 
representation transformation between two mutually conjugate entangled states (s,r'\ and \m,r) , 
this is like the case that the Fourier transform kernel is just the matrix element between the coor- 
dinate state and the momentum state, a wonderful result unnoticed before. Therefore the bipartite 
entangled state representations' transforms, which can lead us to the Hankel transform, was pro- 
posed first in classical optics, can find their way back in quantum optics. 



4.3 Quantum optical version of classical circular harmonic correlation 

From Eq. J47D we can see that its reciprocal relation is the circular harmonic expansion, 

oo 



(66) 



or correlated-amplitude — number-difference entangled state \m,r) can be considered as circular 
harmonic decomposition of \r) = re l °) . Let g (r,0), a general 2-dimensional function expressed in 
polar coordinates, be periodic in the variable 9, it can be looked as the wavefunction of the state 
vector \g) in the (77 = re l6 \ representation 



g(r,9) = ( V = re^ \g) 



using (1661) we have 



9 ( r > °) = 9m ( 



r e 



-irnO 



g m (r) = (to, r \g) 



(67) 



(68) 



g m (r) is the wavefunction of \g) in (m,r\ representation. By using d53D and noticing that It then 
follows from ( 1491 ) 



ia(a 1 a 1 -a 2 a 2 ) 1^ _ 



'ae n — 



rj = re 



77 = re l{e - a) 



(69) 



so e iQ ( aiai a 2 a2 ) behaves a rotation operator in \rj) representation, we see that the expectation 
value f e _M ( a i ai_a 2 a2 ) j n |^ j s 

n {g\ e - ia ^- a l a ^\g) =Tt{g\ J _JL e -MaUl-ala 2 ) |^ ^ |^ 

= / rdrdO / ( 5 r/^re^-^W^re 1 * \g) 
Jo Jo ' 

poo />2lT 

= rdr I g*(r,9-a)g(r,9)d9 = R a , (70) 
Jo Jo 

which is just the cross-correlation between g (r, 6) and an angularly rotated version of the same 
function, g* (r, 9 — a). On the other hand, using (153ft we have 

\g) = ^2 / d(r 2 ) |m,r) (m,r| .g) = ^ / d (r 2 ) |m, r) # m (r) . (71) 

m— — 00 m— — 00 

Substituting Q71D into (T70t and using the eigenvector equation (I5QD as well as J54D we obtain 

i? Q =7r £ / d(r^( m 'y\g^(r')e-H a ^-^) £ / d (r 2 ) |m, r) ,g m (r) 

m' = -oo *'° m=-oo ^° 

= 7r E E / d{r' 2 )g* m ,(r')e-" na d (r 2 ) g m (r) 5 m , n ,-6 (r - r>) 

m'=-oo m=-oo ^° ^° 

00 />OG 

= 2n ]T z~ lma \ r\g m (r) \ 2 dr, (72) 

m=-oo ^° 



from which we see that each of the circular harmonic components of the crosscorrelation undergoes 
a different phase shift —ma, so R a is not rotation invariant. However, when we consider only one 
harmonic component 

/>oo 

R aM = 2ire- lMa / r\g M (r) \ 2 dr, (73) 
Jo 

is extracted digitally, then from the phase associated with this component it is possible to determine 
the angular shift that one version of the object has undergone. When an optical filter that is matched 
to R a ,M °f a particular object is constructed, then if that some object is entered as an input to the 
system with any angular rotation, a correlation peak of strength proportional to / °° r\gM (r) \ 2 dr 
will be produced, independent of rotation. Hence an optical correlator can be constructed that will 
recognize that object independent of rotation H19H . 

So far we have studied the circular harmonic correlation in the context of quantum optics, we 
have endowed the crosscorrelation R a with a definite quantum mechanical meaning, i.e. the overlap 

ia(aa -a a ) 

between (g\ and the rotated state e ia \ aiai ° 2 a2 / \g} ; in the entangled state representation. 

Note that Fourier-based correlators is also very sensitive to magnification, however, the magnitude 
of Mellin transform is independent of scale-size changes in the input H19H . Now we examine when 
\g) is both rotated and squeezed (by a two-mode squeezing operator S 2 (A) = exp[X(a\al — aia 2 )]), 
then from Q37D and ( 1671 ) we have 

S 2 (A) | 5 > = / — \vM (V Is) = / — Wv) 9 (r, 9) , (74) 



1 r 



it follows the overlap between (g\ and the state e -«*(<»i<»i-<*3«2) 5 |^ ( 



ir(g\e 



. / t t > 

-ml ai — a 2 0,2 



5(A) |s) 



= (.9l 



/7 2 ri • ( t t \ 
'/ — ia ai — a 2 a2 J 



\V/H)9(r,0) 



/ ( ff ,/ = e^ ^r/^) g(r,6)d6 

11 In ' 



M Jo 

00 



r°° r dr f 27T 

/ — / g*(r/»,e-a)g(r,6)d9, 
Jo M Jo 



(75) 



which corresponds to the crosscorrelation arising from combination of squeezing and rotation (joint 
transform correlator) . On the other hand, from Q47D and (13 7D we see 



S2 (A) |m,r) 



1 

27T/i 7 



f /(9 



# \ „—imB 



7] = — 6 ) e 

A 4 



(76) 



and therefore 



m, - ) # m (r) 



^ /* / t t \ _ — _ r°° 

W a ,X = - E / d{r' 2 ){ m 'y\g* ml (r')e-< a ^- a ^) £ / d( 

= £ E E / d(r' 2 )g* m ,(r')e— d (r 2 ) 5m (r) (£ - f 



2tt 



E 



' m— — oo 



^5™ (r) 5m' e ) dr, 



(77) 



from which one can see that to achieve simultaneous scale and rotation invariance, a two-dimensional 
object g (r, 6) should be entered into the optical system in a distorted polar coordinate system, the 
distortation arising from the fact that the radial coordinate is stretched by a logarithmic transforma- 
tion (A = — In//), which coincides with Ref. [14511 ■ The quantum optical version is thus established 
which is a new tie connecting Fourier optics and quantum optics H46H . 

At the end of this section, using the two-variable Hermite polynomials' definition [47] 



i(m,n) 



H m ,n (£>£*) — E 



1=0 



m\n\ (-1) <f™"'£ 

n (m - iy. (n - iy. 



(78) 



which is quite different from the product of two single-variable Hermite polynomials, and its gener- 
ating function formula is 



00 ±mff7i 



E 



m.n—0 



mini 



H m ,n (£,, D = exp [-u' +14 + fC] , 



(79) 



and noting if m ,n(£>£*) = e^ m n ^i? mi „(r', r'), we can directly perform the integral in ([55]) and 
derive the explicit form of |s, r') , 

i r 27r 

|s,r'' 



1 

2"?r 7o 



d^exp{-/ 2 /2 + £(4 + £*4 - istp - a\a\} |00) 

fl-K oo fm fn 



z7r Jq ^— ' mini 



m,n—0 
/.2-7T 00 



H m ,n{r',r')e 



m, n) 



r ' 2/2 E 



i!f„+ s , „(/,/) |n + s,n) 



(80) 



which is really an entangled state in two-mode Fock space. Eqs. Q79D and J8QD will be often used in 
the following discussions. 

In the following we concentrate on finding the generalized Fresnel operators in both one- and 
two- mode cases with use of the IWOP technique. 



5 Single-mode Fresnel operator as the image of the classical Op- 
tical Fresnel Transform 

In this section we shall mainly introduce so-called generalized Fresnel operators (GFO) (in one- 
mode and two-mode cases both) H48H and some appropriate quantum optical representations (e.g. 
coherent state representation and entangled state representation) to manifestly link the formalisms 
in quantum optics to those in classical optics. In so doing, we find that the various transforms in clas- 
sical optics are just the result of generalized Fresnel operators inducing transforms on appropriate 
quantum state vectors, i.e. classical optical Fresnel transforms have their counterpart in quantum 
optics. Besides, we can study the important ABCD rule obeyed by Gaussian beam propagation (also 
the ray propagation in matrix optics) fl49H in the domain of quantum optics. 



5.1 Single-mode GFO gained via coherent state method 

For the coherent state \z) in quantum optics ||T1|2H 

\z) = cxp \za) — z*a\ |0) e 



z 
* 

z 



(81) 



which is the eigenstate of annihilation operator a, a\z) = z \z), using the IWOP and (113D . we can 
put the over-completeness relation of \z) into normal ordering 



±±\z){z\ = /^:e-(**-° >- >: =1 
the canonical form of coherent state \z) is expressed as 

|z} = |p,?> = exp[i(pQ-gP)]|0>= ( * 



(82) 



(83) 



where z = (q + ip) /V2. It follows that (p, q\ Q \p, q) — q, (p, q\ P \p, q) = p, this indicates that the 
states \p,q) generate a canonical phase-space representation for a state |W) , "J (p, q) — (p,q\ \f r ) • 
Thus the coherent state is a good candidate for providing with classical phase-space description 
of quantum systems. Remembering that the Fresnel transform's parameters (A,B,C,D) are ele- 
ments of a ray transfer matrix M describing optical systems, M belongs to the unimodular sym- 
plectic group, and the coherent state \p, x) is a good candidate for providing with classical phase- 
space description of quantum systems, we naturally think of that the symplectic transformation 





2 











in classical phase space may mapping onto a generalized Fresnel operator in 

Hilbert space through the coherent state basis. Thus we construct the following ket-bra projection 
operator 



dxdp 



A B q 
C D \p 



(84) 



as the GFO. In fact, using notation of \z) (coherent state), and introducing complex numbers s, r, 

s = \[A + D - i(B - C)] , r = -\[A- D + i(B + C)] , \s\ 2 - \r\ 2 = 1, (85) 



i n 



from 



we know 



A B 
C D 



s —r 
-r s 



= \sz — rz*) 



exp 



\sz — rz*\ 2 + (sz — rz*)a) 



10), 



(86) 



s —r 
—r* s* 



is still a symplectic group element, so J84D becomes H50H 



Fi(s,r) 



d 2 z 



(87) 



where the factor yfs is attached for anticipating the unitarity of the operator F%. Eq. G84D tells us 
that c-number transform (q,p) — > (Aq + Bp, Cq + Dp) in coherent state basis maps into F\ (s, r). 
Now we prove F\ (s, r) is really the FO we want. Using the IWOP technique and Eq. (1861) and (1131) 
we can perform the integral 



Fx (s, r) =y/s I — : exp 

7T 



|s| 2 \z\ 2 + sza) + z* (a — ra^) 



ex P ( _ 7TT at2 ) : ex P \ I ~7 ~ 1 ) aTfl } : ex P ( 7TT a 



2s* 



ex P ( -7TT at2 ) ex P i ( flTa + o ) ln T f - ex P ( 7TT a 



r s 

r* 



2 *2 + 

z H — — jz — a' a 



2s* 



2s* 



2s* 



(88) 



where we have used the mathematical formula |[5TTl 

J ^ ex P {C \z\ 2 +£z + V z* + fz 2 + gz* 2 } 
: exp • 



1 



C 2 - 4/5 



VC 2 - Wg 

with the convergent condition Re(£ ± / ± g) < 0, Re( ^ 2 ± ^ 4 ^ ) < 0. It then follows 

1 



(89) 



z| 2 + |z'| 2 rz* 2 rV 2 z*z' 



Then using 



exp 



-1/4 

7r 'exp 



2s* 



2s* 



2 2 

•- ^ + 

2 2 



2 



(90) 



(91) 



and the completeness relation of coherent state as well as Q85D we obtain the matrix element of 
Fi (s, r) (= Fi (A, B, C)) in coordinate representation (xi\, 



(x 2 \F 1 (s,r)\x 1 )= / — {x 2 \z)(z\F 1 (s,r) / \z') (z'\ Xl ) 



V2iriB 



exp 



2B 



(Ax 2 — 2x 2 x\ + Dx\ 



K.{x 2 ,x 1 ) , 



(92) 



which is just the kernel of generalized Fresnel transform K. (x 2 , xi) in (I12D . The above discussions 
demonstrate how to transit classical Fresnel transform to GFO through the coherent state and the 
IWOP technique. 



Now if we define g (x 2 ) = (2:2 1 g)> f ( x i ) 



txi 



/) and using Eq. J17D , we can rewrite Eq. (II 1 D as 



(x 2 \g)= / dx 1 (x 2 \F 1 (A,B,C)\x 1 )(x 1 \ f) 

J — CO 

= (x2l.F1 (A,B,C) I/) , 



(93) 



which is just the quantum mechanical version of GFO. Therefore, the 1 -dimensional GFT in classical 
optics corresponds to the 1-mode GFO F\ (A, B, C) operating on state vector |/) in Hilbert space, i.e. 
\g) = Fi (A : B : C)\f). One merit of GFO is: using coordinate-momentum representation transform 
we can immediately obtain GFT in "frequency" domain, i.e. 



(P2[F\ P1 ) 

1 



dxidx 2 {p 2 \ x 2 ) (x 2 \ F \xx) (xi \pi 



dx\dx 2 



y/2iriB J- 
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y/2*i (-C) 
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2B 



( o X\ . \ iD 2 

V 1 + ~A ' _ V + 2B X2 ~ IP2X2 



exp 



2 (-CO 



(Dpi - 2 P2Pl + Apj 



Obviously, F\ (A, B, C) induces the following transform 
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F^'iA.B.C) 1 J Fi (A,B,C) 
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(94) 
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5.2 Group Multiplication Rule for Single-mode GFO 

Because two successive optical Fresnel transforms is still a Fresnel transform, we wonder if the 
product of two GFO is still a GFO. On the other hand, we have known that the GFO is the image 
of the symplectic transform z — > sz — rz*, we expect that the product of two symplectic transforms 
maps into the GFO which is just the product of two GFOs. If this is so, then correspondence between 
GFT and GFO is perfect. Using ([87]), (z \z') = exp [-± (\z\ 2 + |z'| 2 ) + z*z'\ and the IWOP technique 
we can directly perform the following integrals 

F 1 (s,r)F 1 (s',r') 



where we have set 
or 

from which we see that it is just the mapping of the above (^4, B, C, D) matrices multiplication. 
Hence F\ (s, r) F\ (s', r') is the loyal representation of the product of two symplectic group elements 
shown in Q98D . The above discussion actually reveals an important property of coherent states, 
though two coherent state vectors are not orthogonal, but the equation 

v 7 ^ 7 / flfll. \ az - rz*) (z \s'z' ~ r'z'*) (z'\ = v 7 ^ 77 / — Is"* z - r" z*) (z\ (99) 



d 2 zd 2 z' 



: exp 



2s"* 



-,t2 



exp 



2s"* 



s ., , _| s "* z _ r " z *)( z | = i? 1 ( s " 5r ") 

7T 



(96) 



(97) 



s —r 
-r s 



= MM', \s"\ 2 - \r"\ 



= 1, (98) 



seems as if their overlap (z \s'z' — r'z'*) was a <5-function. The coherent state representation for 
GFOs' product may be visualized very easily, but it achieves striking importance, because it does not 
change its form when treating symplectic transform according to z — > sz — rz* . 

As a result of this group multiplication rule of GFO, we immediately obtain 

K, M " (x 2 ,x 1 ) = {x 2 \F 1 (s",r")\x 1 ) 

/OO 
dx 3 {x 2 \F 1 (s, r) \x 3 ) (x 3 \ F 1 {s', r') \ Xl ) 
OO 

/OO 
dx 3 K M {x 2 ,x 3 )K. M ' {x 3 ,x x ), (100) 

provided that the parameter matrices (s", r") satisfy J97l ) . Thus by virtue of the group multiplication 
property of GFO we immediately find the successive transform property of GFTs. 



6 Quantum Optical ABCD Law for optical propagation — single- 
mode case 

In classical optics, ray-transfer matrices, N = ( n ^ ] , AD — BC = 1, have been used to 



v C D , 

describe the geometrical formation of images by a centered lens system. For an optical ray (a 
centered spherical wavefront) passing through optical instruments there is a famous law, named 
ABCD law, governing the relation between input ray (n, ol\) and output ray (r 2 , a 2 ) , i.e. 

r - n I. (101) 



a 2 I \ Oi\ 

where n is the ray height from the optical axis, and a\ is named the optical direction-cosine, r\ joi\ = 
Ri specifies the ray's wavefront shape. Eq. (I101D implies 

R 2 = li = ARi+B . (102) 
" a 2 CRx+D 

This law is the core of matrix optics, since it tells us how the curvature of a centered spherical 
wavefront changes from one reference plane to the next. Besides, the multiplication rule of matrix 
optics implies that if the ray-transfer matrices of the n optical components are N±, N 2 , N 3 , ■ ■ ■ , N n , 
respectively, then the whole system is determined by a matrix = NiN 2 N 3 ■ ■ ■ N n . 

One of the remarkable things of modern optics is the case with which geometrical ray-transfer 
methods, constituting the matrix optics, can be adapted to describe the generation and propagation 
of Laser beams. In 1965 Kogelnik 115211 pointed out that propagation of Gaussian beam also obeys 
ABCD law via optical diffraction integration, i.e. the input light field / (x\) and output light field 
q (x 2 ) are related to each other by so-called Fresnel integration H19H g (x?) — K,(A,B,C;x 2 ,xi)f(xi)dxi, 
where 



K. (A, B, C; x 2 , xi) = ■. exp 
y/2iriB 



— (Ax\ — 2x 2 xi + Dx 
The ABCD law for Gaussian beam passing through an optical system is 1153 



2' 



Aqi+B 

92 = c^Td' (103) 

where qi (q 2 ) represents the complex curvature of the input (output) Gaussian beam, Eq. (I103D has 
the similar form as Eq. (I102D . An interesting and important question naturally arises H54B : Does 
ABCD law also exhibit in quantum optics? Since classical Fresnel transform should have its quantum 
optical counterpart? 



To see the ABCD law more explicitly, using Eq. J85D we can re-express Eq. J88D as 

F 1 (A,B,C) 



A + D + i(B-C) 
2 



: exp 



A - D + i [B + C) 



,t2 



A + D + i (B - C) 



- 1 



2 [A + D + i (B - C)] 

A - D — i(B + C) 2 
~ 2[A + D + i{B-C)] a 



(104) 



and the multiplication rule for Fx is F {A', B', C, D') F {A, B, C,D) = F {A", B", C", D") , where 

(105) 
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B" N 
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B>\ 
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I) 



Next we directly use the GFO to derive ABCD law in quantum optics. From Eq. (ll05D we see that 
the GFO generates 



Fx(A,B,C)\0) = 
if we identify 



■ exp 



then 



A - D + i {B + C) 
A + iB-i(C + iD) t 2 [A + D + i (B - C)} ' 

A- D + i(B + C) _ qi -i 
A + D + i(B-C) ~ qi +i' 

g] 

2 (qi + i 



10} , 



F l( W)| ) = J-W^e,p 



The solution of Eq.( fT07l ) is 



5i 



A + iB 



|0) 



(106) 



(107) 



(108) 



(109) 



Let Fi (A, 73, C) |0) expressed by (11081) be an input state for an optical system which is characteristic 
by parameters A',B',C',D', then the quantum optical ABCD law states that the output state is 



Fi (A , B , C ) F\ {A, B, C) |0) = 



-2/{C" + iD") 



q2 + i 



exp 



g2 - I |2 
2(q 2 + l) a 



|0) 



which has the similar form as Eq.( 1108D , where (C", D") is determined by Eq. (1105D . and 

A'qi + B' 



'72 



C'gi + D' 

which resembles Eq. (1103t . 
Proof: 

According to the multiplication rule of two GFOs and Eqs.( |104b - (I105D we have 



(110) 



(HI) 



Fi (A , B , C ) F\ (A, B, C) |0) 

/ 2 f A" - D" + i{B" + C") 

~ V A" + D" + i {B" - C") 6XP \ 2 [A" + D" + i (B" - C") 



,t2 



10} 



A' (A + iB) + B'(C + iD) - iC> (A + iB) - iD' (C + iD) 



x exp 



A' (A + iB) + B'(C + iD) + iC (A + iB) + iD' (C + iD) f2 l 
2 [A' (A + iB) + B' (C + iD) - iC" (A + iB) - iD' (C + iD)} a 1 1 " 



-2/ [C + iD) 



A' qi -B' + i (C'qx - D') ts 



A'gi —B' — i {C qi - D') 6XP ( 2 - B' - i (C'qx - D')] ° J 



(112) 



Using Eq.(fT09l) we see 2/(c+t - p) 



r ,,, - -2/ (C" + iD") . together using Eq. dTTTT) we can reach Eq. JTTOl ), 
thus the law is proved. Using Eq. J1Q8D we can re-express Eq. fllllD as 



<12 



A'(A + iB)+B'{C + iD) A" + iB" 



(113) 



C'{A + iB) + D'{C + iD) C" + iD' r 

which is in consistent to Eq.( |109D . Eqs. ( 1108D - (I113D are therefore self-consistent. 

As an application of quantum optical ABCD law, we apply it to tackle the time-evolution of a 
time-dependent harmonic oscillator whose Hamiltonian is 



H = \e-^P 2 + ^e^Q 2 , h=l, 



(114) 



where we have set the initial mass mo = 1, 7 denotes damping. Using u (t) — e^'^e (Q p + p Q) to 
perform the transformation 



u(t) Qu- 1 (t) = e~ 7t Q, 

u(t) Pu- 1 (t) = e 7t P - 7 e 7 *Q, 



then i^H^- = H |V> (t)) leads to • 



dt 



du- 1 (t) _ 1 



H^n = u(t) Hu- 1 (t) - m (i) = + ^ 2 Q 2 , 



(115) 



(116) 



where uj 2 



7 2 . H does not contain t explicitly. The dynamic evolution of a mass-varying 



harmonic oscillator from the Fock state |0) at initial time to a squeezed state at time t is 



W(t)) = u- 1 (t)\0) = e- 



-(QP+PQ), 



if we let A = D = 1,B 



0,C 



-7; and A' 
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(117) 
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~zt, according to Eq. flllOl) we directly obtain 



0, then qi = -±- v 



u- 1 (t)\0) = 



e -2 7 t _|_ j 7 _|_ 1 



exp 



-27* _ 1 



^_ fl t2 



2 (e- 2 ^* + 1 + ij) 



10) 



(118) 



so the time evolution of the damping oscillator embodies the quantum optical ABCD law. 



7 Optical operator method studied via GFO's decomposition 

Fresnel diffraction is the core of Fourier optics H191 l30l [331 [36ll . Fresnel transform is frequently 
used in optical imaging, optical propagation and optical instrument design. The GFT represents a 
class of optical transforms which are of great importance for their applications to describe various 
optical systems. It is easily seen that when we let the transform kernel K, (X2, x\) — exp (1x2X1), the 
GFT changes into the well-known Fourier transform, which is adapted to express mathematically 
the Fraunhofer diffraction. And if K, (x2,X\) — exp[i (x2 — xi) 2 ], the GFT then describes a Fresnel 
diffraction. In studying various optical transformations one also proposed so-called optical operator 
method 115511 which used quantum mechanical operators' ordered product to express the mechanism 
of optical systems, such that the ray transfer through optical instruments and the diffraction can be 
discussed by virtue of the commutative relations of operators and the matrix algebra. Two important 
questions thus naturally arises: how to directly map the classical optical transformations to the 
optical operator method? How to combine the usual optical transformation operators, such as the 
square phase operators, scaling operator, Fourier transform operator and the propagation operator 
in free space, into a concise and unified form? In this section we shall solve these two problems and 
develop the optical operator method onto a new stage. 



00 



7. 1 Four fundamental optical operators derived by decomposing GFO 

The GFO Fx (A, B, C) can also be expressed in the form of quadratic combination of canonical 
operators Q and P H56L i.e., 



iC 



2A 



Fx (A, B, C) = exp — Q z exp -- (QP + PQ) \nA exp -— P 



iB 



2A 



(119) 



where we have set h = 1, A ^ 0. To confirm this, we first calculate matrix element 

{x\F 1 {A,B,C) \ P ) =exp( l ^-x 2 ) (x\exp ( ~ (QP + PQ)lnA) |p)exp 
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iB 
2A 1 



1 



: exp 



V^kA 

where we have used the squeezing property 



iC 2 iB 2 ipx 
2A X ~2A P + ~ 



exp 



-- (QP + PQ)\nA 



\p) = -= \p/A) . 
VA 



It then follows from G120D and AD - BC = 1, we have 

/OO 
dp (x 2 \ Fx (A, B, C) \p) (p\ xx) = JC M (x 2 ,xx) 
-OO 



(120) 



(121) 



(122) 



Thus Fx (A, B, C) in (I119D is really the expected GFO. Next we directly use (I84D and the canonical 
operator (Q, P) representation (I119D to develop the optical operator method. 

By noticing the matrix decompositions [13111 



A B 
C D 



1 
C/A 1 



A 

A" 1 



1 B/A 
1 



and comparing d84b and J119D as well as using (11051) we know 

Fx (A, B, C) = Fx (1, 0, C/A) Fx (A, 0, 0) Fx (1, B/A, 0) : 

where 



(123) 



(124) 



^(1,0,0/^) = 



y/2 + iC/A 
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1 

C/A 1 



exp 1 2T 



(125) 



which is named quadrature phase operator; and 

y/2-iB/A 



Fx(l,B/A,0) 



2V2ir 
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exp ( -— A P 



(126) 



which is named Fresnel propagator in free space; as well as 

F x (A, 0, 0) = ^^J 4 1 / dxdp 



2V2ir 
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A- 1 



exp 



-^(QP + PQ)\nA 



(127) 



which is named scaling operator (squeezed operator l|9l 11010 . When A = D = 0,B = 1,C = — 1, 
from Q104D we see 



1 

-1 



exp 
exp 



a^a + — ) bxi 



,-k ( + 1 
-i— a 1 a H — 
2 V 2 



(128) 



which is named the Fourier operator, since it quantum mechanically transforms ||57l 
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(129) 



7.2 Alternate decompositions of GFO 

Note that when A = 0, the decomposition jll9D is not available, instead, from 
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(130) 



and (flQ4b , ([1051) and JTT91 ) we have 
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it then follows 
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Besides, when we notice 
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and 
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we have another decomposition for B^O, 
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and for C ^ 
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7.3 Some optical operator identities 



For a special optical systems with the parameter A = 0, C = -B 



we have 
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In particular, when A = D = Q,C = — B 1 , from 
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we have 
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it then follows that 

F 1 (A, B, C) = exp 

while from 
we obtain 
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so we have 
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In this section, based on a one-to-one correspondence between classical Fresnel transform in 
phase space and quantum unitary transform in state-vector space and the IWOP technique as well 
as the coherent state representation we have found a way to directly map the classical optical trans- 
formations to the optical operator method. We have combined the usual optical transformation 
operators, such as the square phase operators, scaling operator, Fourier transform operator and the 
propagation operator in free space, into a concise and unified form. The various decompositions of 
Fresnel operator into the exponential canonical operators are also obtained. 



8 Quantum tomography and probability distribution for the Fres- 
nel quadrature phase 



In quantum optics theory all possible linear combinations of quadratures Q and P of the oscillator 
field mode a and a) can be measured by the homodyne measurement just by varying the phase 
of the local oscillator. The average of the random outcomes of the measurement, at a given local 
oscillator phase, is connected with the marginal distribution of Wigner function (WF), thus the 
homodyne measurement of light field permits the reconstruction of the WF of a quantum system by 
varying the phase shift between two oscillators. In Ref. H58H Vogel and Risken pointed out that the 
probability distribution for the rotated quadrature phase Qg = [cfl exp(i8)+aexp(— i0)]/\/2, [a, a 1 ] = 
1, which depends on only one 9 angle, can be expressed in terms of WF, and that the reverse is also 
true (named as Vogel-Risken relation), i.e., one can obtain the Wigner distribution by tomographic 
inversion of a set of measured probability distributions, Pg (qg) , of the quadrature amplitude. Once 
the distribution Pg (qg) are obtained, one can use the inverse Radon transformation familiar in 
tomographic imaging to obtain the Wigner distribution and density matrix. The Radon transform of 
the WF is closely related to the expectation values or densities formed with the eigenstates to the 
rotated canonical observables. The field of problems of the reconstruction of the density operator 
from such data is called quantum tomography. (Optical tomographic imaging techniques derive two- 
dimensional data from a three-dimensional object to obtain a slice image of the internal structure 
and thus have the ability to peer inside the object noninvasively, the slice image is equated with 
tomogram.) The theoretical development in quantum tomography in the last decade has progressed 
in the direction of determining more physical relevant parameters of the density from tomographic 
data 115811591 l60ll6Tll62ll. 



8.1 Relation between Fresnel transform and Radon transform of WF 

In H631 [64l the Radon transform of WF which depends on two continuous parameters is introduced, 
this has the advantage in conveniently associating quantum tomography theory with squeezed co- 
herent state theory. In this subsection we want to derive relations between the Fresnel transform 
and the Radon transform of WF in quantum optics in tomography theory. 

By extending the rotated quadrature phase Qg to the Fresnel quadrature phase 

Q F = (s*a + ra ] + sa ] + r*a) / V2 = FxQF\, (146) 
where s and r are related to ABCD through J85l), 

s = \ [A + D - i (B - C)] , r = [A - D + i (B + C)] , \s\ 2 - \r\ 2 = 1, (147) 

we shall prove that the (D, B) related Radon transform of Wigner operator A (q,p) is just the pure 
state density operator \q) srsr (q\ (named as the tomographic density operator) formed with the 
eigenstates belonging to the quadrature Qp, ( \q) s r = F\ \q) , Q is the coordinate operator), 

/oo 
dq'dp'5 [q - (Dq' - Bp')] A (q',p') , (148) 
-OO 

D = -(s + s* +r + r*), B=— (s* -s + r* - r) , (149) 
2 2i 

Since F corresponds to classical Fresnel transform in optical diffraction theory, so Eq. J148D indicates 
that the probability distribution for the Fresnel quadrature phase is the Radon transform of WF [16511 . 

Proof: 

Firstly, from A88D we see 

F 1 (s,r)aFl (s,r) = s*a + ra\ (150) 



so from Q = . P = i a J 1 , indeed we have 
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FxQFl = ft-JL—pT = ( s *a + ra f + sa f + r* a) /V2 = Qj?. 

Secondly, we can derive the explicit form of \q) s r . Starting from s* +r* = D + iB, s* — i 
we set up the eigenvector equation 



QF\q). r = (DQ-BP)\q). r = q\q). r , 



it follows 



\q) S!r = F 1 (s,r)\q). 
In the coordinate and momentum representations we have 
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The normalizable solutions to (I154D and (I155D are 
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Using the Fock representation of \q) and \p) in Eqs.fO and (1231) . we obtain 
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Comparing Eq. G158D with Q159D we see 
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exp -^q" ^rpr . (160) 

On the other hand, according to the orthogonalization of \q) , sr (q' \q") s r — 5(q' — q") , we have 
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Thus combining Eq. (1 1601) and (11611) we deduce 
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It is easily seen that that \q) g make up a complete set (so |g) can be named as the tomography 
representation), 

|>00 

(165) 
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dq\q) s>rS ,r (q\ = 1- 
Then according to the Weyl quantization scheme H37H 

dpdqA(q,p)h(q,p), 

-oo 

where /i is the Weyl correspondence of iJ (Q, P) , 

/i(«/,p) = 27rT r [H(Q,P)A( ? ,p)], 
A (q,p) is the Wigner operator Il66l|6"7ll , 



(167) 
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and using J167D , (1168D and J158D we know that the classical Weyl correspondence (Weyl image) of 
the projection operator \q) s rs r (q\ is 
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Combining Eqs. (I150I) - (I153I) together we complete the proof. Therefore, the probability distribution 
for the Fresnel quadrature phase is the Radon transform of WF 



| (q\ Fl |V) | 2 = |., P (q\ V) | 2 - / dq' dp' 6 [q - (Dq' - Bp')} (^| A (q',p') ty) , 



(171) 



so we name \q) s rs r (q\ the tomographic density. Moreover, the tomogram of quantum state \ip) 
is just the squared modulus of the wave function s ^ r (q\ ip) , this new relation between quantum 
tomography and optical Fresnel transform may provide experimentalists to figure out new approach 
for generating tomography. 

The introduction of \q) s r also bring convenience to obtain the inverse of Radon transformation, 
using (11651) we have 



-igQt 



dq\q) sr s,r (q\e- ig « 



dqdpA (q,p) e -i9(Dq-Bp) 



(172) 



Considering its right hand-side as a Fourier transformation, its reciprocal transform is 

-| pOO pOO pTT 

A(q,p) = —j dq' J dg'\g'\ J d^) s<rS>r {cl\ 



x exp 



-w 



y/D 2 + B 2 



q cos ip — p sin if 



(173) 



where g' = g\jD 2 + B 2 , cos ip = -^=5=5. simp — ^ D ^ g2 -So once the distribution | s>r (q\ ip) | 2 are 
obtained, one can use the inverse Radon transformation familiar in tomographic imaging to obtain 
the Wigner distribution. 

By analogy, we can conclude that the (A, C) related Radon transform of A (q,p) is just the pure 
state density operator |p) s rs (p\ formed with the eigenstates belonging to the conjugate quadrature 
of Q F , 

/OO 
dq'dp'8[p-{Aq' -Cp')]A{q',p'), 
-oo 

A = — (s* — r* + s — r) , C = — (s — r — s* + r*) . 
2 2i 

Similarly, we find that for the momentum density, 

/oo 
dq'dp'6\ P -(Ap'-Cq')]A(q',p') 
-00 

where 

D + iBp 2 V2~ip + A + iCa) 2 ^ 

J- n ' 4- 



(174) 



(175) 



F l \P) = \P) S r = I \ eX P ' 



|0>. 



(176) 



y/A-iC 1 A — iC 2 A-iC A-iC 2 

As an application of the relation J148D , recalling that the Fi (r, s) makes up a faithful representation 
of the symplectic group [?], it then follows from G148D that 

^'(r'y^iM) I?) (q\F}(r,s)F?(r',s') = \q) s „ y , s „. r „ (q\ 

dq'dp'8 [q - ((B'C + DD') q' - (AB' + BD')p')] A (q',p') , (177) 

In this way a complicated Radon transform of tomography can be viewed as the sequential operation 
of two Fresnel transforms. This confirms that the continuous Radon transformation corresponds 
to the symplectic group transformation 11631 16411 , this is an advantage of introducing the Fresnel 
operator. The group property of Fresnel operators help us to analyze complicated Radon transforms 
in terms of some sequential Fresnel transformations. The new relation may provide experimentalists 
to figure out new approach for realizing tomography. 

8.2 Another new theorem to calculating the tomogram 

In this subsection, we introduce a new theorem, i.e., the tomogram of a density operator p is equal 
to the marginal integration of the classical Weyl correspondence function of F^pF, where F is the 
Fresnel operator. 

Multiplying both sides of Eq. J170D by a density matrix p and then performing the trace, noting 
the Wigner function W(p, q) — Tr [pA(p, q)] , one can see 



f f dqdpS [q - (Dq- Bp)] A (q, p) p 

J J —OO 



Tr 



= Tr \q) 



W\P = 



\q) s , r = (q\F^pF\q) 



dq'dp6[q - (Dq- Bp)]W{p,q). 



(178) 



The right hand side of Eq. Q178D is commonly defined as the tomogram of quantum states in (B, D) 
direction, so in our view the calculation of tomogram in (B, D) direction is ascribed to calculating 

(q\ F^pF \q) = S. (179) 

This is a concise and neat formula. Similarly, the tomogram in (A, C) direction is ascribed to 

(p\F^ P F\p). 

According to the Weyl correspondence rule 



H(X,P)= dpdxt)(p,x)A(p,x), (180) 

J J — CO 

and the Weyl ordering form of A(p, q) 

A(p,q)=[6(q-Q)6(p-P)], (181) 

where the symbol ; ; denotes Weyl ordering, the classical correspondence of a Weyl ordered operator 
; f)(Q, P) ] is obtained just by replacing Q — > q, P — > p in h, i.e., 



'MQ,P)[= dpdqt ) (p,q)A(p,q), (182) 

J J — oo 

Let the classical Weyl correspondence of F^pF be h{p, q) 

F^pF= ff dpdqh(p,q)A(p,q) 1 

J J —OO 

then using ( fT79l ) and J168D we have 

E=(q\F^ P F\q) 

= (q\ JJdpdq%(p,q)A(p,q)\q) 

= J J dpdqhip, q) J + ™ ^e*» (X q+ |) (q- | q) 

= J I dpdq H p, q) £" ^6 (q-q + l)s(q-q- *-) 

= ^JJ dpdq-h(p,qy 2 ^'-^S(2q- 2q) = (183) 

Thus we reach a theorem: 

The tomogram of a density operator p is equal to the marginal integration of the classical Weyl 
correspondence h(p,q) of F^pF, where F is the Fresnel operator, expressed by 



Tr 



I>\'1K,,, i'/l = / (184) 



or 

dx 

Tr 



In this way the relationship between tomogram of a density operator p and the Fresnel transformed 
p's classical Weyl function is established. 



o 1 



9 Two-mode GFO and Its Application 



For two-dimensional optical Fresnel transforms (see (O) in the x — y plane one may naturally think 
that the 2-mode GFO is just the direct product of two independent 1-mode GFOs but with the 
same (A : B : C, D) matrix. However, here we present another 2-mode Fresnel operator which can 
not only lead to the usual 2-dimensional optical Fresnel transforms in some appropriate quantum 
mechanical representations, but also provide us with some new classical transformations (we name 
them entangled Fresnel transformations) . 



9.1 Two-mode GFO gained via coherent state representation 

Similar in spirit to the single-mode case, we introduce the two-mode GFO F 2 (r, s) through the 
following 2-mode coherent state representation Il50l 

f d 2 zid 2 Z2 

F 2 (r,s)=s J — \szi + rz 2 ,rz 1 + sz 2 ) {Z\, z 2 \ , (186) 

which indicates that F 2 (r, s) is a mapping of classical sympletic transform (zi,z 2 ) — > (szi + rz 2 ,rz\ + sz 2 ) 
in phase space. Concretely, the ket in ( 11861 ) is 

\sz\ + rz\ , rz\ + sz 2 ) = \szx + rz 2 ) 1 <8> \rz\ + sz 2 ) 2 , ss* — rr* = 1, (187) 

s and r are complex and satisfy the unimodularity condition. Using the IWOP technique we perform 
the integral in ( 1186D and obtain 



s| 2 (|zi| 2 + \z 2 \ 2 ) - r*sz lZ2 - rs*z* lZ * 2 



F 2 (r,s)—s J —d 2 zid 2 z 2 : exp[ 

+ (sz\ + rz 2 ) a\ + {rz\ + sz 2 ) a\ + z\a\ + z 2 a 2 — a\a\ — a\a 2 ] : 



1 ' r ,t„t 



— exp [ —a[a 2 ) : exp 



1 



— -111 a\ai + a\a 2 



r 

exp I a±a 2 



exp ^-^aja^ exp[^a{ai + a\a 2 + ij In (s*) 1 ] exp ( — ^—a\a 2 \ . (188) 



vs* / v / \ s* 

Thus F 2 (r,s) induces the transform 

F 2 (r, s) oiFa -1 (r, s) = s* ai - ra\,F 2 (r, s) a 2 F~ x (r, s) = s*a 2 - ra\. (189) 

and F 2 is actually a general 2-mode squeezing operator. Recall that J37D implies the intrinsic relation 
between the EPR entangled state and the two-mode squeezed state, which has physical implemen- 
tation, i.e. in the output of a parametric down-conversion the idler-mode and the signal-mode 
constitute a two-mode squeezed state, meanwhile are entangled with each other in frequency do- 
main, we naturally select the entangled state representation to relate F 2 (r, s) to two-dimensional 
GFT. Letting \g) = F 2 (r, s) |/} , and then projecting onto the entangled state (77'! defined by Q29D 
and using the completeness relation (13 ID ot\rj), we obtain 

g{ v ') = {r 1 , \g) = {r, , \F 2 {r,s)\f) 

^ (t/I F 2 (r, s) \n) ( V \f) = J d 2 V 4 r - s) (r,', V ) f ( V ) . (190) 
Then using the overcompleteness relation of the coherent state, and 

F 2 (r, s) \ Zl ,z 2 ) = 1 exp I -h\ Zl \ 2 + \z 2 \ 2 + \z[f + \z' 2 f) 



s* [ 2 

-—z[*z 2 - —ZlZ 2 + — {z'^zi + z 2 z 2 ) 
s* s* s* 



O ' 1 



we can calculate the integral kernel 

4^ (v',v) = -W\F 2 (r,s)\ V ) 

d 2 Z\d 2 z 2 d 2 z[d 2 z 2 , , , 



(v'\ 4>4) ( z i>4l F 2 (r,s) \z 1 ,z 2 ) (zi,z 2 \rj) 



1 f d z\d 2 z 2 d 2 z[d 2 z' 2 



exp 



-(|z 1 | 2 + |z 2 | 2 + |4| 2 + |4| 2 )-i(b/| 2 + H 2 ) 



x exp 



- — ZlZ 2 +Z 1 Z 2 + TJZi + — Zl + z 2 z 2 ) + — *i z 2 + + Tj z 1 -r] z 2 -ri z 2 



1 



(— r — s + r* + s*) 7r 



exp 



(-s + r*) \t]\ 2 - (r + s) |r/'| 2 + W '* + 77*77' 1 



-7" — s + r* 



(W\ 2 + \vf 

(191) 



Using the relation between s, r and (A, B, C, D) in Eq. (|85D Eq. J191D becomes 



(192) 



where the superscript M only means the parameters of are [A, -B;C, D], and the subscript 2 
means the two-dimensional kernel. Eq. Q192D has the similar form as (I12D except for its complex 
form. Taking 771 = x±, rj 2 = x 2 and r/[ = x[, r/ 2 = x' 2 , we have 



^2 (v'iV) =^2 (x'i,x' 2 ;xi,X2) =1C 1 (x^xfj) (gi/Cj (x 2 ,x' 2 ). 



(193) 



This shows that F2 (r, s) is really the counterpart of the 2-dimensional GFT. 

If taking the matrix element of F 2 (r, s) in the |£) representation which is conjugate to \rj), we 
obtain the 2-dimensional GFT in its 'frequency domain', i.e., 

d 2 r)'d 2 r) ,, 



(e\F 2 (r,s) |0 



1 



8iBn 
1 



(0| 77') (a\F 2 (r,s)\r,) ( v \ 

9 ^2 vl ? */J ex P 



2i(-<7)7T 



exp 



2 (-CO 



-^(0,0, (194) 



where the superscript N means that this transform kernel corresponds to the parameter matrix 

N = [D, —C, —B, A]. 

The two-mode GFO also abides by the group multiplication rule. Using the IWOP technique and 
01861 ) we obtain 



F 2 (r, s)F 2 (r',s') 

, f d 2 Z\d 2 z 2 d 2 z[d 2 z' 2 



ss 



: exp{-\s\ 2 (\ Zl \ 2 + \z 2 \ 2 ) -r*sz lZ2 



1 



* * * *- n / |2 . 1 / |2 1 I / / 1 I i*\2 ,111*, 1 1 |2i 

-rs z 1 z 2 ~-[\z 1 \ +\z 2 \ +\sz 1 +rz 2 \ +\r Zl +sz 2 \\ 

+ (szi + rz 2 ) a\ + {rz{ + sz 2 ) a\ + z'*ax + z 2 a 2 

+ z{ (s'z[ + r'z 2 ) + z 2 (r'z 1 * + s'z 2 ) - a\a x - a\a 2 } : 



■ exp 



2.s 



a\a\ + a 2 a 2 



exp 



2s h 



-a\a 2 



F 2 (r",s"). 



(195) 



where (r", s") are given by Eq.([97l) or ( |98l ). Therefore, Q195D is a loyal representation of the multi- 
plication rule for ray transfer matrices in the sense of Matrix Optics. 



9.2 Quantum Optical ABCD Law for two-mode GFO 

Next we extended quantum optical ABCD Law to two-mode case. Operating with F 2 (r, s) on two- 
mode number state |m, n) and using the overlap between coherent state and number state, i.e. 



(z 1 ,z 2 \m,n) 



\Jm\n\ 



exp 







-5(1*1"- 


- M 2 ) 







and the integral formula H57II 



d 2 z 



(196) 



(197) 



we can calculate 

F 2 (r, s) | to, n) 

d 2 z\d 2 z 2 



\sz\ + rz 2 ,rz{ + sz 2 ) (z\,Z 2 \m,n) 



d 2 zid 2 Z2 



z 1 z 2 exp 



s [ zi + z 2 



\Jm\n\ J 7T 2 

x exp —sr*z\z 2 — rs* z\z\ + (szi + rz 2 ) a\ + (?z\ + sz 2 ) a\ |00) 



I i2m+2 / i T 

\s\ ymmi 



1 



,t. 



t„t 



ra\ a 



sr z 2 ) exp ( — \z 2 \ +—a 2 z 2 - 



1"2 



|00) 



s* n + 1 y/mln\ 



exp 



ra\a 2 



|00> 



(198) 



where H m ^ n (e, e) is the two variables Hermite polynomial Il63l 1 6411 , shown in d78D and Q79D . 
Using Eqs.((85]) and ( fT09l ), we recast Eq.( fT98l ) into 

-2/(C + iD) ( ql+iC-iD\ n 



F 2 (r,s) \m,n) = 



(qi + i) \Jra\n\ \ qi + i C + iD 

2a\/(C + iD) 2a\/{C-iD) 



exp 



~^a\4) \m. (199) 



qi+i q x + 1 

Noticing the multiplication rule of F 2 (r, s) in Eq. J195D , which is equivalent to 

F 2 (A 1 , B', C) F 2 (A, B, C) = F 2 {A", B", C") , (200) 

where (A', B', C) , (A, B, C) and (A", B", C") are related to each other by Eq.(fT05]). 

Next we directly use the GFO to derive ABCD rule in quantum optics for Gaussian beam in two- 
mode case. According to Eq. Q198D and Eq. (I200D we obtain 

F 2 (A',B',C')F 2 {A,B,C)\m,n) 



exp 



r"a\a\ 



s''* n+1 VmM 

Similar to the way of deriving Eq. Q199D , we can simplify Eq. Q201D as 

F 2 (A',B\C')F 2 (A,B,C)\00) 

-2/ (C" + iD") ( q^+i C" - iD"^ n 



|00), 



(201) 



( ?2 + i)vfi! V q 2 + iC" + iD" 

2a\/{C" + iD") 2a\/{C"-iD") 



?2 + i 



q* +i 



exp 



92 - i t t 
a a., 

qi +i 1 2 



1 00), 



(202) 



where the relation between q 2 and qi are determined by Eq. tllllD which resembles Eq. (ll03D . this is 
just the new ABCD law for two-mode case in quantum optics. 
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9.3 Optical operators derived by decomposing GFO 
9.3.1 GFO as quadratic combinations of canonical operators 

In order to obtain the quadratic combinations of canonical operators, let first derive an operator 
identity. Note Q l = (ai + a\)/\/2, Pi = (a, - a|)/(v / 2i), and Eq.([30]),([3l]) we can prove the operator 
identity 

e |[(Qi-Q2) 2 + (Pi+P2) 2 ] = J ^Z e #[(Qi-Q 2 ) 2 +(Pi+P2) 2 ] |^ ( ?? | 



1-A 



: exp 



2A 
1-A~ 



-K, 



(203) 



where we have set 



K + = ^[(Q 1 -Q 2 ) 2 + (P 1 +P 2 ) 2 }. 



When B = 0, A = 1,C ^ C/A, D = 1, and using Eq. tfTSSl ) we see that 



(204) 



Ft (1,0, C/A) 



exp 



iC/A 

2 - iC/A ' \ 2 - iC/A 



(a\ai + a 2 a 2 — a\a\ — aia^ 



(205) 



which is corresponding to the square phase operator in single-mode case. In a similar way, using 
and Q41D we can derive another operator identity 



[(Qi+Q2) 2 +(Pi-p 2 ) 2 ] = f d2 £ 4[(Qi+Q2) 2 +(Pi-P2) 2 



10 <£| 



1 - A 



exp 



2A 



-K_ 



(206) 



where 



K- = J[(Qi + Q 2 ) 2 + (P 1 -P 2 ) 2 ]. 



It then follows from Eqs. (fT88l) and ([2061) 

2 



F 2 (1,B/A,0) 



2 + iB/A 
2 

2 + iB/A 



exp 
exp 



iB/A 
2 + iB/A 

2iB/A 
2 + iB/A 



(a\a\ + a\a 2 + a\a\ + aia^ 



(207) 



. iB 
exp<j —£ K - 



(208) 



which is corresponding to Fresnel propagator in free space (single-mode case) . In particular, when 
B = C = 0, and D = A' 1 , Eq. ( flBBl becomes 



F 2 (A, 0,0) = sechA: exp — a\a 2 tanh A + (sech A — 1) \a\a\ + a} 2 a 2 )+ a i°2 tanh A 



(209) 



where = tanh A, A = e A . Eq. J209D is just the two-mode squeezing operator, 

F 2 (A, 0,0) = exp [t (QiP 2 + Q 2 Pi)lnA] = exp [-2K \aA] , 

K = ~(QlP2 + Q2Pl), 



(210) 
(211) 



o r- 



which actually squeezes the entangled state |£) (its conjugate state is \r])), 

F 2 (A, 0,0) \0= f^j W/A) m = ± \Z/A) . 



(212) 



Using the decomposition (I123D of the matrix and combining equations (1205D . J208D and (I210D 
together, we see that 



F 2 (A, B, C) = F 2 (1,0, Cj A) F 2 (A, 0, 0) F 2 (l,B/A, 0) 

{iC 1 ( iB 

— K + > exp {-2i^o In A} cxp j --r^- 

This is the two-mode quadratic canonical operator representation of F 2 (A, B, C). 
To prove Eq. (|213l ), using (I2T21 and ([42]) we see 

(771 F 2 (A, B, C) |0 = exp (g h| 2 - g |£| 2 ) („| 1 ^ (£'| 

1 exp H 2 - i? |^ i/A) 



A 



2A 



2A 



1 fiC . |2 lt|2 ', 
2A eXP Ul |7?l "2A lfl )eXP 



It then follows 



(r/| F 2 (A, B, C) It?) = f° ^ (r/| F 2 |£) (£| r?) 



2iB 



exp 



2~B 



(213) 



(214) 



(215) 



which is just the transform kernel of a 2-dimensional GFT and the definition given in (I213D is true. 

Note that the quadratic combinations in Eqs. d204l ), Q207D and (I211D of the four canonical op- 
erators (Qi, Q 2 ; Pi, P 2 ) obey the commutative relations [K + ,KJ\ = 2K , [K ,K±] = ±K±, so 
F 2 (A,B,C) involves a SU(2) Lie algebra structure (this structure is also compiled by Q 2 /2, P 2 /2 
and —i (QP + PQ) /2 that have been used in constructing F± (A, B, C)). 



9.3.2 Alternate decompositions of GFO and new optical operator identities 

When A = D = 0,B = 1, C = -1, from Eq. (fl88l we see 

F 2 (0, 1, —1) = exp — (a\ai + a\a 2 + 1^ In i 
= exp 



7T / + + 

-i— ya[ai + a 2 a 2 + 1 



(216) 



which can also be named the Fourier operator, since it induces the quantum mechanically transforms 
11571 

PQ,T = P u PP t P = -Qi. (217) 

it then follows that 



(218) 



On the other hand, in order to obtain the decomposition of F 2 (A, B, C) for A = 0, similar to 
deriving Eq. (1132t . we have 



F 2 (A, B, C) = exp 



iB 
~D 



exp [2Kq In D] exp 



iC 
~D 



/V4 



, for D ^ 0. 



(219) 



o r 



While for B ^ or C 7^ 0, using Eqs. (ll33D and J134D we have another decomposition of F 2 (A, B, C), 
i.e., 



and 



F 2 {A, B, C) = exp 



F 2 (A,B,C) =exp 



— K- 
B 



cxp [—2Kq In B] J 7 exp 





r iA ' 
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~c K - 
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-2K Q In -1 


J 7 exp 


-f-1 



, B^Q, 



,C^0. 



(220) 



(221) 



In addition, noticing Eqs. ( 11431 ) and ( I141D , we can rewrite Eqs. (12201) and (I221D as follows 



and 



F 2 (A,£,C) =exp 



F 2 (A, S, C) = exp 



- (D - 1) 



(A - 1) JT_ 



cxp [— iBK—] exp 



exp [iCif+j exp 



^(A-1)K + 



C 



( D - 1) K. 



(222) 



(223) 



respectively. 

Next, according to some optical systems used frequently in physical optics, we derive some new 
entangled optical operator identities. For a special optical system with the parameter A = 0, C = 
—B^ 1 , ( 1137D which corresponds to the Fourier transform system, we have 



exp 



exp [2Kq In D] exp 









-4d k \ 


= exp 






_-B K+ _ 



exp[-2i4T ln B] T. 



(224) 



In particular, when A = D = 0, C = —B 1 , Eq. fll39D corresponding to the ideal spectrum analyzer, 
we have 



B 



cxp [-iBK-] = cxp [-2K In B] T. 



exp [— iBKJ\ exp 
WhenB = 0, D = A -1 , 

(A \ _ / 1 
\C A- 1 )-{ Cj A 1 

which corresponds to the form of image system, another operator identity is given by 

~iC. 



(225) 
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WhenC = 0, A = D~ 1 , 
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(226) 
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which corresponds to the far foci system, 
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1 B/D 
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(227) 



When D = 0, C = — £> 1 , corresponding to the Fresnel transform system, 



A B 

-B- 1 



B 
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1 



we have 

exp 
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tK 4 



exp [— 2Kq In A] exp 



"T 



-1 J \ A/B 1 
exp [— 2Kq In B] J- exp 
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IT 



-K, 



(228) 



The GFO can unify those optical operators in two-mode case. Various decompositions of the GFO 
into the exponential canonical operators, corresponding to the decomposition of ray transfer matrix 
[A, B, C, D] , are also derived. In our derivation, the entangled state representation is of useness in 
our research. 



9.4 Quantum tomography and probability distribution for the Fresnel quadra- 
ture phase — two-mode entangled case 

In section 8 we have found that under the Fresnel transformation the pure position density \q) (q\ 
becomes the tomographic density \q) s rs r (q\, which is just the Radon transform of the Wigner op- 
erator A (q,p) . In this section we want to generalize the above conclusion to two-mode entangled 
case. Here we shall prove 

F2 \v) (v\ F 2 = \v) s ,rs,r (v\ = * j <*W<5 {m - Da 2 + B 7l ) 8 { m - Dai - B l2 ) A {a, 7) , (229) 

i.e., we show that \r]) g rs (77 1 is just the Radon transform of the entangled Wigner operator A (a, 7) . 

Similar in spirit to the single-mode case, operating F 2 (r, s) on entangled state representation 
we see 
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where we have used the integration formula 

— exp K |z| 2 + £z + r;z* ) = 



-i e -^,Re(C) <0. 

Noticing the completeness relation and the orthogonality of \rj) we immediately derive 

l»?>..r..r fal ^ L -,r (v\ v') s , r = ttS ( V - rf) 5 {rf - r/*) , 



(230) 
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a generalized entangled state representation \ rj) s r with the completeness relation (I233D . From (I231D 
we can see that 



a i l^f 

«2 \v) t 



V 



D — iB 



D + iB D + iB 



4 ) l»7>.,r> 



D-iB 



4 |»7). r> 



'> D + iB ' D + iB' 
so we have the eigen-equations for \rj) as follows 

[D (Qi -q 2 )-b (p 1 - p 2 )} \ v ) Sir = V2 Vl \r,) B , . 

[B (Qi + Q 2 )+D (P 1 + P 2 )] \ v ) s r = ^2 V2 \ v ) s r . 
We can also check Eqs. Q234l )-( |237D by another way. 



(234) 
(235) 

(236) 
(237) 



o n 



9.4.1 \rj) r _ {rj\ as Radon transform of the entangled Wigner operator 



For two-mode correlated system, we have introduced the Wigner operator in Q45D . According to the 
Wely correspondence rule H37B 



H 



(a\, a 2 ;a 



i,a 2 



d 2 jd 2 ah (cr, 7) A (cr, 7) 



where h (cr, 7) is the Weyl correspondence of H (a\,a\; a\,a^ , and 

h (cr, 7) = 47r 2 Tr H (a\, a\; 01, a 2 ^j A (cr, 7) 
the classical Weyl correspondence of the projection operator \rf) sr {rj\ can be calculated as 

\v) r ,sr,s(v\&(^j) 

d 2 rf 



(238) 



(239) 



47r 2 Tr 



= 4tt^ 



7T 
72 



r, s (»7| cr - r/) (cr + rf\ Tj) expfaV - 7/'* 7) 



= ^ 2 J^T (V\ 4 \a - rf) (cr + r/| F 2 |ry) expfoV - r/* 7 ). 
Then using Eq. ( 11921 ) we have 

4vr 2 Tr [|»j}, jr , jr (r?| A (cr, 7)] = irS (r? 2 - Dcr 2 + 5 (77! - Da x - B l2 ) , 
which means the following Weyl correspondence 

l*?>. rs.r (V\=K d 2 1 d 2 a5 (r? 2 - Da 2 + B71) <5 (r?i - £><ti - £ 72 ) A (cr, 7) , 



(240) 



(241) 



(242) 



so the projector operator |t?) s rs r (ry| is just the Radon transformation of A (cr, 7), D and B are the 
Radon transformation parameter. Combining Eqs. (12301) -( |242l ) together we complete the proof 
Q229D . Therefore, the quantum tomography in two-mode entangled case is expressed as 

kr (v\ V>) I 2 = I (v\ |V>) I 2 = / d 2 1 d 2 a5 ( m - Da 2 + B 7l ) <5 ( m - Dui - £ 72 ) (^| A (a, 7) |^) . 



(243) 

where (ip\ A (a, 7) \ip) is the Wigner function. So the probability distribution for the Fresnel quadra- 
ture phase is the tomography (Radon transform of the two-mode Wigner function). This new relation 
between quantum tomography and optical Fresnel transform may provide experimentalists to figure 
out new approach for generating tomography. 

Next we turn to the "frequency" domain, that is to say, we shall prove that the (A, C) related 
Radon transform of entangled Wigner operator A (cr, 7) is just the pure state density operator 

I0 s , rs ,r (f I > Le '> 

F2 10 (CI Fl = \i) s rs r (£| = 7T J 6 (6 - Aa 1 -C l2 )5 & - Aa 2 + C lx ) A (a, 7) d 2 ad 2 7 , (244) 

where |£) is the conjugated entangled state to |tj). 

By analogy with the above procedure, we obtain the 2-dimensional Fresnel transformation in its 
'frequency domain', i.e., 

K»(Z>,Z) = l{Z>\F 2 (r,s)\Z) 

IT 

d 2 rjd 2 a , ,, , 



(Civ') (v'\F 2 (r,s)\r,) (r,\ 



1 



8iBn 
1 



d 2 ad 2 r] 



■ exp 



Kg'" (a, V ) 



2i(-C)TT 



exp 



2(-C) 



(245) 



where the superscript N means that this transform kernel corresponds to the parameter matrix 
N = [D, —C, —B, A]. Thus the 2D Fresnel transformation in its 'frequency domain' is given by 

*(£') = / ^(^,6^(6^- (246) 

Operating F 2 (r, s) on |£) we also have 

1 f D + iB , , 2 fa\ Cat A + iC + + 1 , nArn 

®" =A^C C ^{-HA-JC) W + T^C + Hc- A-JC^) ■ (247) 

or 

\0 S r = "T~~ — * exp{- f+ r * Id 2 + + 4A; 4^44) 100) . (248) 

s* — r* 2 (s* — r* ) s* — r* s* — r s* — r* 

Noticing that the entangled Wigner operator in (£| representation is expressed as 

A (*> 7) = / ^ |7 + (7 - £| exp(£V - (249) 
and using the classical correspondence of |£) rs (£| which is calculated by 

/ l (<T, 7 )=47r 2 Tr [|£) s , rs , r (e|A(<7, 7 )" 

= 4 1 ^< 7 -£|F 2 |£) (e | F t| 7 + 0ex P (ra-a*0 

= tt5 (6 - Aa t - C l2 ) 5 (f 2 - Aa 2 + C71) , (250) 

we obtain 

I0 i>r ,,r {£1=*/* (6 - - C 7 2 ) 5 (6 - Aa 2 + C 7l ) A (a, 7 ) dVd 2 7 , (251) 

so the projector operator |£) (£| is another Radon transformation of the two-mode Wigner op- 
erator, with A and C being the Radon transformation parameter ('frequency' domain). Therefore, 
the quantum tomography in s r (£| representation is expressed as the Radon transformation of the 
Wigner function 

| (£| Ft |V>) | 2 = | a , r (f | V) I 2 - tt J d 2 -yd 2 aS (ft - An - C 72 ) 5 (6 - Aa 2 + C 7 i) (VI A fa 7) IV) , 

(252) 

and s , r (£| = (C|Ft. 

9.4.2 Inverse Radon transformation 

Now we consider the inverse Radon transformation. For instance, using Eq.( |242D we see the Fourier 
transformation of \r)) Sirs r (r)\ is 

d2 V \v) s ,r S ,T (v\ exp(-iCi?7i - ^2%) 
= 7T / d 2 jd 2 aA fa 7) exp [-<i (£><n + B l2 ) - i( 2 (Da 2 - B 7l )\ , (253) 

the right-hand side of (I253D can be regarded as a special Fourier transformation of A fa 7), so by 
making its inverse Fourier transformation, we get 

-1 />00 />00 />7T 

A ( (J '7) = 7 ^ I / A-ilnl/ dr 2 |r 2 | / d0i£0a 
f°° d 2 n 

x / —\n). r .M K (ri,r 2 ,d u 9 2 ), (254) 

J-00 7T ' ' 



where cos 9\ = cos 



VB 2 +D 2 

K (ri,r 2 ,9 1 ,Q 2 ) = exp 
x exp 



n = Ci V-B 2 +D 2 ,r 2 = C2VB 2 + D 2 and 



-VT\ 



-IT 2 



VB 2 + D 2 

m 

VB 2 + D 2 



o\ cos f 1 — 72 sin ti\ 
a 2 cos 62 + 71 sin 62 



(255) 



Eq. (I254D is just the inverse Radon transformation of entangled Wigner operator in the entangled 
state representation. This is different from the two independent Radon transformations' direct 
product of the two independent single-mode Wigner operators, because in J231D the \rj) sr is an 
entangled state. Therefore the Wigner function of quantum state can be reconstructed from the 
tomographic inversion of a set of measured probability distributions \ 8iT (rj 

1 



i.e., 



W, 



(2*)< 



dri |n| 



dr 2 \r 2 \ / dflidfla 
10 



/°° r\ 2 n 
^| s , r (r, \iP)\ 2 K(r u r 2 ,e 1 ,62). 
-00 T 



(256) 



Thus, based on the previous section, we have further extended the relation connecting optical 
Fresnel transformation with quantum tomography to the entangled case. The tomography repre- 
sentation s . r = (rj\ f\ is set up, based on which the tomogram of quantum state \ip) is just the 
squared modulus of the wave function S}7 . (rj\ tp) . i.e. the probability distribution for the Fresnel 
quadrature phase is the tomogram (Radon transform of the Wigner function) . 



10 Fractional Fourier Transformation (FrFT) for 1-D case 

The fractional Fourier transform (FrFT) has been shown to be a very useful tool in Fourier optics and 
information optics. The concept of FrFT was firstly introduced mathematically in 1980 by Namias 
H24H as a mathematical tool for solving theoretical physical problems [68], but did not brought 
enough attention until Mendlovic and Ozaktas 11251 ?] defined the a-th FrFT physically, based on 
propagation in quadratic graded-index media (GRIN media with medium parameters n(r) = n\ — 
n 2 r 2 /2). Since then a lot of works have been done on its properties, optical implementations and 
applications EH [73 E0 E3 ■ 

10.1 Quantum version of FrFT 

The FrFT of 0-order is defined in a manner, i.e., 

TtffV J e<f ~ e) f°° f - x2 + y 2 , ixx '\ <i w rim 
Fe / (x)} = \ - — — / exp <^ -1— — - + ^— \ f (x) dx, (257) 

V 2tt sm J_oo 1^ 2 tan ti sin J 

where the exponential function is an integral kernel. In order to find the quantum correspondence 
of FrFT, multiplying the function exp |— i 2 tan 9 + sSfe } ^ tne / dy\y) and bra J dx (x\ from left 
and right, respectively, where \y) and \x) are coordinate eigenvectors, X\x) = x\x), and then using 
dl4D and the IWOP technique to perform the integration, we obtain 

dxdy |y) exp + S} (x| 

= sin 6e lf> : exp { (e l9 - l) a f a} : 

= \J sin tie 10 exp { i9a f a } , (258) 



where we have used the operator identity in the last step of Eq. J258D 



exp {fa 1 a} =: exp { (e f - l) a f a} : . (259) 
From the orthogonal relation (x' \x) — S (x — x') , we know that Eq. (|258D indicates 



"" e) exp|-^ + S^^l^ ta k) (260) 



27rsin6' 2 tan 9 sin( 



which implies that the integral kernel in Eq. d257D is just the matrix element of operator exp {iOc^a} 
in coordinate state (exp {iOcJa} called as Fractional Fourier Operator 117311 - Therefore, if we consider 
/ (x) as (x \f), the wave function of quantum state |/) in the coordinate representation, from Eqs. 
( [2571 ) and ( fSSol it then follows 



OO 



Fe[f(x)]= dx(y\e^ a \x)f(x) = (y\e Wa ' a \f)=g(y), (261) 

J — OO 

which suggests 

\g) = e lQa ' a \f) ■ (262) 
From Eqs. d26l| ) and Q257D one can see that the FrFT in Eq. J257D corresponds actually to the rotating 
operator (e l8afa ^J transform in Eq. (l26lD between two quantum states, which is just the quantum 
version of FrFT. 

In fact, using quantum version of FrFT, one can directly derive various properties of the FrFTs. An 
important feature of the FrFT is that they are composed according to Tg'Te = J-e<+e (the additivity 
property) . Without losing generality, we examine 

Fe+o' [f {x)\ ee r (y\ A 9 ^)^ \x) f (x) . 



(263) 

According to the completeness relation of coordinate eigenvector, dx' \x') (x'\ — 1, Eq. (12631) 
yields 

poo 

Tow [/ (*)] - / dx (y\ e ^ f V eVa \x) f (x) 



— OO 
OO 



dx' (y\ e l6ata \x') / dx (x'\ e w ' a ' a \x) f (x) 

i J — OO 

/OO 
dx' (y\ e Wa " a \x') F e [/ (x)} = TaT v [f (x)} , (264) 
-OO 

which is just the additivity of FrFT. 

In particular, when |/) is the number state, |/) = \n) = |0), its wavefunction in coordinate 
representation is 

/ (x) = (x \n) = r =L= e- x2 / 2 H n (x), (265) 

the FrFT of (x \n) is 

To [(x \n)\ = (y\ e Ba ' a |n) = e mB (y \n) , (266) 

or 

H n (x)] = e me e~y 2 ' 2 H n {y), (267) 



To 



which indicates that the eigenfunction is Hermite-Gaussian function with the corresponding eigen- 
value being e m0 . 



A O 




Figure 1: A thick lens as a kind of fractional Fourier transform device. 



10.2 On the Scaled FrFT Operator 

In studying various optical transformations the optical operator method is proposed H74II as map- 
ping of ray- transfer ABCD matrix, such that the ray transfer through optical instruments and the 
diffraction can be discussed by virtue of the commutative relations of operators and the matrix alge- 
bra. The square phase operators, scaling operator, Fourier transform operator and the propagation 
operator in free space have been proposed in the literature, two important questions thus naturally 
arise: 1. what is the scaled FrFT (SFrFT) operator which corresponds to the SFrFT's integration 
kernel (24j 

1 \i(x 2 +x' 2 ) ix'x 1 , nrm 



\f2mfe sin cf> [ 2/ e tan0 / e sint 

where f e is standard focal length (or a scaled parameter); 2. If this operator is found, can it be 
further decomposed into simpler operators and what are their physical meaning? Since SFrFT has 
wide application in optical information detection and can be implemented even by using a thick lens 
H76L so our questions are worth of paying attention H75H - 

Let start with a thick lens (shown in Fig.l) which represents a transfer matrix H76B 



•A- D \ I 1 Kj „ 



C V ) I _[( n _i) Mi 1 _(±^Mi 



(269) 



where n is the reflective index; I is the thickness of thick lens; R\ and R2 denotes the curvature 
radius of the two surfaces of the lens, respectively. When we choose R 1 = R 2 = R, then Eq. ( 12691 ) 
reduces to 

1 _ (l-l/n)i 1 











I) 


■( 



-[(—1)1-^] 1-^ 



(270) 



By defining 1 - il^lMi = C os<t>, I = f e sin0, and ^ = " (1 t ; to 1 s0) , / = nf e sm(f>, we can recast ( [2701 
into the simple form 

According to J104D we immediately know that the operator of SFrFT is 

Fi (i,B,C) = exp <^ — a 1 

f / + l\ , /2sec< 
x exp < I a 1 a + — I In 



2/ V V" 

fi{fe~ l// e ) tan 2 
x exp <^ — }> , (272) 

F= [2 + i(/ e + l// e )tan0]. 
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(273) 



Noting that the matrix ( ^ p can ^ e decomposed into 

1 \ / COS0 \ / 1 ietanc 

-^tan0 1 J ^ cos0 J ^ 1 

according to the previous section we have 

Fx (A,B,C) = F fl,0,-yt&xi<p \ F (cos 0, 0, 0) P (1, f e tan <f>, 0) 

= exp (|^Q 2 ) exp {-1 (QP + PQ) lncos^j exp (A^p^) , (274) 



where Q = (a + at) P = (a - a*) / (y^i) and exp (-^^Q 2 ] , exp {-§ (QP + PQ) lncosf} 
and exp ^_Lk4pL£p 2 ^ are the quadrature phase operator, the squeezing operator and the free prop- 
agation operator, respectively. On the other hand, from ^ ^ p ^ = (^C A J ' an< ^ ^^H) 
we see 

Pf 1 (i,B,C) = exp (-T^Q 2 ) exp (-1 (QP + PQ) In p) exp (^§P 2 ) , (275) 
it then follows 

P {A, B,C) = exp (^^P 2 ) exp Q (QP + PQ) In cos 0) exp (|^Q 2 ) ■ (276) 

Using the canonical operator form Q274D or Q276D of Pi (A,B,C) we can deduce its matrix element 
in the coordinate states \x) (its conjugate state is \p)) 

<z'| P (AS,C) |«) = i , exp { f ^f ^ - -^-"l ) , (277) 
y/2nif e sin ^ y 2/ e tan / e sm / J 

which is just the kernel of SFrFT, thus we name Pi (A,B,C) SFrFT operator. Noticing that the 
Q 2 /2, P 2 /2 and | (QP + PQ) construct a close SU(2) Lie algebra, we can put Eq. d276D into a more 
compact form, i.e., 

Pi (A,B,C) = exp L-itk fp* + <pj j , (278) 

Eqs. Q274D , J276D and Q278D are different forms of the same operator of SFrFT. Especially, when 
f e = 1, Pi (A, B, C) — > exp {— i(j>a^a) , which is the usual FrFT operator. 

Using J277D the SFrFT of / (x) — (x\ f) , denoted as Fj [f (x)} , can be expressed as an matrix 
element in quantum optics context, 



[/ (*)] = J dx W\ p i (AB,C) \x) (x\ f) = (x'\ Pi (A, B,C) I/) 



(279) 



The above discussions are useful since any unimodular matrix y ^ ^D^j Can d ecom P ose d ^ nt0 
EZJ 

A B \ _ / 1 \ / ra \ f coscf) / e sin< 
C D ) ~ \ -V 1 ) \ m- 1 ) I -sin0// e cos0 



(280) 



where the parameters m,V,(p are all real, 



2 , 2 -S 2 B AC + DBl fl , 
m 2 =A 2 + — , tan0=— , P = -f^. (281) 
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Correspondingly, the operator of FrFT is given by 

Fi (A, B, C) — Fi (1,0, -V) Ft (m, 0, 0) Fi (A 6, C) 

= exp (-~PQ 2 ) cxp (-1 (QP + PQ) lnm) cxp (> + ^) } , (282) 

where Fi(l,0, — V) = exp [-^VQ 2 ] is the quadratic phase operator. Thus the general Fresnel 
transform can always be expressed by SFrFT as follows 



g(x')= I dy(x'\ F x (1,0, -V) Fx (m,0,0)\y) I dx (y\F x {A,B,C)\x) f (x) 

= J dx"dy(x'\exp(-^TX 2 ^j \mx") (x" \y) j dx {y\F x {A,B,C)\x) f (x) 
= exp(-^V 2 ) J J^5(^-y) J dx (y\F x (A,B,C)\x)f(x) 
-Vx ,2 \ I dx[— F x (A,B,C)\x)f(x) 
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2 J J \m 

i.e., the output g (x 1 ) is the SFrFT of the input / (a;) plus a quadratic phase term exp (— §Pa/ 2 ) . 

10.3 An integration transformation from Chirplet to FrFT kernel 

In the history of developing optics we have known that each optical setup corresponds to an optical 
transformation, for example, thick lens as a fractional Fourier transformer. In turn, once a new 
integration transform is found, its experimental implementation is expected. In this subsection we 
report a new integration transformation which can convert chirplet function to FrFT kernel H78H , 
as this new transformation is invertible and obeys Parseval theorem, we expect it be realized by 
experimentalists . 

The new transform we propose here is 

^-e 2 ^-^-^h(p,q)^f(x,y), (284) 

which differs from the usual two-fold Fourier transformation J ^^^■e ipx+iqv f(x, y). In particular, 
when h(p,q) = 1, Eq. J284D reduces to 

*d? e 2i(p-*)(,-i0 = r dqS ( g _ y ) e -*ri(9-») = i, (285) 

ft J-00 

so e 2 ^ p ^ x ^ q ^ v ^ can be considered a basis funtion in p — q phase space, or Eq. Q284D can be looked as 
an expansion of / (x, y) with the expansion coefficient being h(p, q) . We can prove that the reciprocal 
transformation of Q284D is 

dxdy e . 2l(p . x)(q . y)f{ ^ y) = ^ y (2g6) 

7T 

In fact, substituting A284D into the left-hand side of ( 12861 ) yields 

dp ' dq '-h(p',q') [[ dxd y c H(p'~x)(o'-y)-(p-^-y)} 



ft .1.1 ft 

2i(p'q'-pq) 



= 11 dp'dq'h(p',q')e 2 

J J — OO 

x5(p-p')6(q-q') = h(p,q). (287) 



This transformation's Parseval-like theorem is 



dpdq 



d P d( l 2i[(-y'p-x'q) + (py+xq)] 



x 5 (x ~ x') 8 (p - p') 



dxdy 



\f(x,y) | 5 



(288) 



Now we apply Eq. J284D to phase space transformation in quantum optics. Recall that a signal 
%j) (q)'s Wigner transform ||62l [66l [67l [79ll is 

du 



(289) 



Using Dirac's symbol H80B to write tp (q) = (q\ ijj) , \q) is the eigenvector of coordinate Q, the Wigner 
operator emerges from Q289D , 



1 

2^ 



due~ lpu 



q -2/\ q+ 2 



A(p,q), h=l. 



(290) 



If h (q,p) is quantized as the operator H (P, Q) through the Weyl-Wigner correspondence H37H 

H{P,Q)= ff dpdq A (p, q) h (q, p) , (291) 

J J — oo 



then 



h (q,p) = / due 



q+ 2 



H(Q,P) 



q -2 



(292) 



this in the literature is named the Weyl transform, h (q, p) is the Weyl classical correspondence of the 
operator H (Q, P). Substituting ( I292D into ( 12841 ) we have 



dpdq 



2i(jp-x)(q-y) 



h(p, q) 



dpdq 



2i(p-x)(q-y) 



x (q + 



OO /» OO 



H(Q,P) 



Using (y + u\ = {u\ e iPy and (v 7 ^)" 1 e~ ix ' 

dueT lxu (y + u\ = 



due~ lpu 

J —OO 

H(Q,P) q 
dq du (q 

— oo J — OO 

xs ( q - y -l) 

due- lxu {y + u\H{Q,P)\y). (293) 

(j>=x \ u) , where (p\ is the momentum eigenvector, and 

due- ixu (u\ e lPv 



u 
2 

2ix(q-y) 



2ir / du{p =x \ u) (u\e lPy 



(294) 



a r 



then Eq. Q293D becomes 

f r dpdq e 2i( P - x )( q -y) h ^ q) = ^ {p=xl (Q, P) \ y ) e **v, (295) 

thus through the new integration transformation a new relationship between a phase space function 
h(p, q) and its Weyl-Wigner correspondence operator H (Q, P) is revealed. The inverse of J295D , 
according to (12861) . is 

f r ^y e -^ P -^-y) {p=x \H (Qj P) \ y) e i*y = h{pj q y (296) 

J J -co y 7T/2 

For example, when H (Q, P) = e ^ p2+< ^ 2 ^ 1 ^ 2 , its classical correspondence is 

e /(p= + Q 2 -i)/ 2 = exp ( P 2 + , 2 )} . (297) 

Substituting Q297D into Q295D we have 



= V^{p= x \e^ p2+Q2 -^ /2 \y)e^. (298) 
Using the Gaussian integration formula 

J j°° ^ e 2 i ( P -x)(?-y) e -A( P 2 + 9 2 ) 

1 i-\(x 2 +y 2 ) 2i\ 2 \ 

in particular, when A = -itan (f - f ) , with ^ = = 1 - J^, Eq. ((2991) becomes 

2 f [°° d P d 1 c 2i(p-x)(Q-y) 

- ia + lJJ-oo * 

exp{i(p 2 + g 2 )tan(^-|)} 

= ^J— ex P { - 4^ ) (300) 

Vie^ M sina [ 2 tan a sua a J 

where expjitan — ^) (p 2 + g 2 )} represents an infinite long chirplet function. Comparing (I300D 
with J2981 ) we see ie~ ia = e f , / = i (f - a) , it then follows 



ie" iQ + 1 
x 



1 



1 i ( x2 +V 2 ) \ 
ex P \ ~n7 1 - ? ) ( 301 ) 



\/27rie~ iQ sin a 1 2tana sina 



where the right-hand side of (1301D is just the FrFT kernel. Therefore the new integration transfor- 
mation (1284D can convert spherical wave to FrFT kernel. We expect this transformation could be 
implemented by experimentalists. 

Moreover, this transformation can also serve for solving some operator ordering problems. We 
notice 

- exp[2i (p -x)(q- y)} 

TT 

dv / v \ 

[q-y- 2 J ex P {* (p ~ x ) u i ' ( 302 ) 
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so the transformation Q284D is equivalent to 



h(p,q) 



dpdq 



2i{p-x){q-y) 
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OO 
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(303) 



For example, using Q290D and (I302D we have 

dpdq 



A(p, g) -> 



Z7T Z 
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q u\ / q u 

/OO j /*oo 
-oo 27r J-oo 



g u 



g it 

' (/+ 2"2/\ y+ 2 + 2 
° |y) (y + U | = |y> <y| 

.^o Z7T 



= 5 (y - Q) 6 (x - P) , 



UU c iu(P-u) 
27T 



so 



dp'dg' 



A (g',p') e 21 (p-p')(9-9') = 5 (g - g) 5 (p - P) , 



(304) 



(305) 



thus this new transformation can convert the Wigner operator to 6 (q — Q) S (p — P) . Similarly, we 
have 

- ff dp'dg'A (g>') e - 2i ( p - p ')( q - q ') =5{p-P)8{q-Q). 



Then for the Wigner function of a density operator p, W^,(p, q) = Tr [pA(p, q)] , we have 

dp'dq' 



-Tr [pA(p',q')}e 2i ( p ' p ')( q " 1 ') 



Tr [ p 8(q-Q)6(p-P)} 
dudv 



Air 2 



■Tr 



p e i(i-Q) u e i{p~ p ) v 



(306) 



we may define Tr [pe l( - q Q^ u e^ p as the Q — P characteristic function. Similarly, 

r °° dp'dq' 



-Tr 



[pA{p',q')} e - 2i (p-p')i<!-l') 



= Tr [pS (p-P)S(q-Q)} 
dudv 



4tt 2 



-Tr 



p e i(P~ P)v e i(q-Q)u 



(307) 



we name Tr [pe l ( p p ) t 'e^ 9 "^"l as the P — Q characteristic function. 



11 Complex Fractional Fourier Transformation 



In this section, we extend 1-D FrFT to the complex fractional Fourier transformation (CFrFT). 



11.1 Quantum version of CFrFT 

According to Ref. H81L based on the entangled state \rj) in two-mode Fock space and its orthonormal 
property, we can take the matrix element of exp — ia (a\ai + a^a^) in the entangled state \rj) , 

JC F (r/,77) — (r/| exp -ia {a\a x + a\a^ \rj) , (308) 
as the integral transform kernel of CFrFT, 



F a [/] (r?) 



Aw, - 



l)a\ 



a i 



(309) 



and the completeness 



Using the normally ordered expansion of e Aa i Ql =: exp 

relation of the coherent state representation, \zi) = exp < — ~ \zi\ 2 + ^a|| |0) i , we calculate that 

K. F (77', i]) is 



7T^ 

d 2 z 1 d 2 z 2 . . , . . 
x / s z 2 ) (21,22] 77) 



exp 



1) (4 



a\ + 0,1(12 



7T^ 

e i(a-§) 

2 sin a 



exp 



*(|gf + M?) 
2 tan a 



2 sin a 



(310) 



which is just the integral kernel of the CFrFT in H82L Thus we see that the matrix element of 
exp — ia (a\ai + a 2 a,2j between two entangled state representations and \r)') corresponds to 
CFrFT. This is a new route from quantum optical transform to classical CFrFT transform. Let r\ = 

X2 + iyz, if =x\ + iyi, Q309D becomes 



Fa [f] (X2,V2) = 



2 sin a 



exp 



+ yl) 



2 tan a 



dxidyi 
x / exp 



i (x\ + y\) . ( Xl x 2 + ym) 



sma 



2 tan a 

In fact, letting / (77) = (rj\ /) , then using Eqs. (l3"TT) and (I31QD we have 

(7/I exp —ia (a\a\ + a^azj \f) 



(311) 



7T 

»(«»-§) 



(?/| exp 



-ia I a\ai + a, 2 a,2 



\V) (V\ I) 



1,2 n v - , , 



2 sin a 



■ exp 



i(W\ 2 + M 2 ) i(ri'*V + V*v') 



2 tan a 



2 sin a 



/(»/)• 



Thus the quantum mechanical version of CFrFT is given by 



Fa [/] (V) = (»/| exp -ia fajai + 4<z 2 J \f) 



(312) 



(313) 



The standard complex Fourier transform is J r 7r / 2 - -7-0 is the identity operator. 



11.2 Additivity property and eigenmodes of CFrFT 

We will show later that this CFrFT can help us to reveal some new property which has been over- 
looked in the formulation of the direct product of two real FrFTs 118311 . The definition Q323D is of 



a r\ 



course required to satisfy the basic postulate that T a Tp [/] (r)') = J- a +p [f (n)] ( rrie additivity prop- 
erty). For this purpose, using Eq. (l323D and Eq. (l31D we see 

T a+P [/ (77)] = (r,'\ e -i(«+«(4«x+4« 3 ) I/) 

r°° d 2 n" t i t ^ 

J-oo T 

^L( v »\ e -*P(4a 1+ ala 2 ) | » j / j 
7T 



— oo 



> I" 

= J- a J> [/ (77)] . (314) 

This derivation is clear and concise by employing the |?7) representation and quantum mechanical 
version of CFrFT. 

On the other hand, the formula (I313D can help us to derive CFrFT of some wave functions easily. 
For example, when |/) is a two-mode number state \m, n) = a\ m a^ 71 / V mini |00), then the CFrFT of 
the wave function (77! m, n) is 

T a [(77I m,n)} = ( ?? '| e '(«+«H«i+4^) | m>n ) 

= e *(<*+/3)(™+«) (r)'\m,n). (315) 

To calculate (77'! m, n), let us recall the definition of two-variable Hermite polynomial H m , n (£,£*) 
(??), we can expand (77' | as 

(7/| = (00| V i m+n ^H m , n (-«/*,»»/) eH" I /2 , (316) 
* — ' mini 

m.n—O 

thus 

(77'! m,n) = —i==H m , n (-177'*, 777') e~\ v I /2 . (317) 
v mini 

As a result of J317D we see that equation A315D becomes 



#™,n (-«?*, «7) e" |,)|2/2 l = e^ + W ro +™>£r m ,„ (-*»/*,«/) e-|"'| 2/2 . (318) 



If we consider the operation J" Q as an operator, one can say that the eigenfunction of JF a (the 
eigenmodes of CFrFT) is the two-variable Hermite polynomials ff m ,„ with the eigenvalue being 
e »(a+/3)(m+n) j jj^js j s a new p r0 p er ty of CFrFT. Since the function space spanned by H m , n (n, 77*) is 
complete, 



and 



J -V^'i^n (77,77*) [H m , n (n,r]*)}* = Vm\nlm'\n'l5 m , m >6 n , n ,, (319) 

OO 1 

V ^H m , n (n,n*)[H m , n (n\n'*)]*e-M 2 =n5(n-n')S( V *-r 1 '*), (320) 



TO!rt 
m,n— 



one can confirms that the eigenmodes of CFrFT form an orthogonal and complete basis set 118411 - 
Note that the two variable Hermite polynomial H m n (77, 77*) is not the direct product of two inde- 
pendent ordinary Hermite polynomials, so CFrFT differs from the direct product of two FrFTs. 



11.3 From Chirplet to CFrFT kernel 

In this subsection, by developing Eq. Q284D to more general case which can be further related to 
the transformation between two mutually conjugate entangled state representations |£) and I77), we 
shall propose a new integration transformation in £ — 77 phase space (see Eq. (I321D below) and its 
inverse transformation. We find that Eq. (I321D also possesses some well-behaved transformation 
properties and can be used to obtain the CFrFT kernel from a chirplet H85H . 



1- r\ 



11.3.1 New complex integration transformation 

Corresponding to the structure of phase space spanned by |£) and \rj) and enlightened by Eq. J284D . 
we propose a new complex integration transformation in £ — rj phase space 

f^l e (A-^v'-^)-(v-UC-^) T{r]i tf= D t v M ) . (321) 
When J"(£,7?) = 1, (132TT ) becomes 



cP^cPrj 



^-m)(.v*- v *)-{v-^KC-h*) 



J d 2 £S (£ - n) 6 (C - if) e^-^-^tt-ri = 1, (322) 

so e^~^ v ~ v _ ^*) can be considered a basis function in £ — r\ phase space, or Eq. (I321D 

can be looked as an expansion of D (y, (j,) in terms of e^ -1 *'^ ~ v '^-"'it ~v ) ) with the expansion 
coefficient being J-(rj, £). 

We can prove that the inverse transform of (I321D is 



/ 



7T^ 

In fact, substituting (13211) into the left-hand side of (13231) yields 

^*f) lisp. 

X e («'-M)('7'*-«'*)-('7'-«')(f'*-A'*)+(«*-A»*)('?-«')-('7*-«'*)«-M) 
2~ 



d £d T) , gs e (t' v '*-ri'£'*+ev-V*i) 
7T 2 ' 

!^!^!^ e (')*-'?'*)M+(')'-'?)M* e (c'*-?*)''+(c-?') ,y * 

= / d 2 e'd 2 ? 7 'j-(7/,e')e (5 '' , '^'"' ? '* +r " _ ' , * ?) 



x 5« - ,,) 5( 2 ) 5 K - O = HV, ■ (324) 
This Parseval-like theorem for this transformation can also be demonstrated, 



^ 2 



x exp [(//V - i/» + (/iV* - i///*)] 

2~ ex P f ~ &)V + (p- P>)V 

xexp[(^-^*K + (^'-^)r] 
d 2 iid 2 v 



-\D{v,n)Y II d^d z v' 

x exp [(/xV — ^V) + (//z/* — f'//*)] 
x^ 2 ' (ji-n')8^S{i/ -u) 



D{v,ix)\\ (325) 



11.3.2 Complex integration transformation and complex Weyl transformation 

In Ref. H86H for correlated two-body systems, we have successfully established the so-called entan- 
gled Wigner operator, expressed in the entangled state (rj\ representation as Q45D , 

A (a, 7 ) -> A(r h £) = / \r, - a) (r, + a\ e^"^^, (326) 

the advantage of introducing A (77, £) can be seen in Ref. H87L The corresponding Wigner function 
for a density matrix p is 

W P (V, = J (V + "I P \v - °) (327) 
If F(rj, £) is quantized as the operator F (Qx, Q 2 , Pi, Pa) through the Weyl-Wigner correspondence 

P(Qi,Q 2 ,Pi,P 2 ) = f d 2 r!d 2 &(r),0&(ri,0, (328) 

then using (I327D we see 

T( V , = 47r 2 Tr [F (Qi, Q 2 , Pi, P 2 ) A(r?, 0] 

= 4 / ^e CT «*- CT *«(7 ? + a|P(Q 1 ,Q 2 ,P 1 ,P 2 )|7 7 - ( r), (329) 

which is named as the complex Weyl transform, and iF(r], £) is the Weyl classical correspondence of 
F (Qi, Q 2 , Pl, Pa). Substituting J3291 into J32TT > we get 

y ^ ^ e (e-rt( 1 ,*-,«)-(,-,)(f- (1 «)_ F ^^ 

d 2 ^ 1 ! JZ-Mn'-v')-(v-v)(e-n*) 
72 



/j2 
(77 + ct| F (Qi, Q 2 , Pi, P 2 ) |ry - a) 

x ( ?7 + CT|P(Q 1 ,g 2 ,P 1 ,P 2 )| ? 7- ( T) 



>2 

= 4 / — e" v -" ff * (^ + 2 CT |P(Qi,Qa,Pi,P2)k). (330) 

7T 



Using d30D . we have 



/- + 2n = ,2a ,. M , <j J_ [i/i (P x - P 2 ) - v 2 (Qx + Q a )] } • (331 ) 



v = V\ + iv 2 . 

As a result of (f33Tb and ie M * <T_A1<T * = (£=^| 2cr) , we see 



4 / (fae^ 17 -^" (i/ + 2cr| 

= % j d 2 o (£ =M | 2a) (2c7|exp{-j= (Pi - P 2 ) - i/ 2 (Qj + Q 2 )}} 
= 2vr (f =|t | exp{-J= [u t (P 1 - P 2 ) - v 2 (Q, + Q 2 )}} 

= 2ir(Z =fl \e i ^ 2Ul -> llU2 ' ) . (332) 



Using (I332D . we convert Eq. (I330D as 



d 2 £d z rj 



{$-»)(ri*-v*)-(n-vW-n*) 



F(V,0 



= 2tt (£ =p | F (Qi, Q 2) Pi,P 2 ) H e *(^w-^/*0. 
The inverse of Q333D , according to Q323D , is 



(333) 



2(ffxd 2 v (e _^ )(v _ u) _ (ri *_ u * )($ _ fl) 



(^\F(Q 1 ,Q 2 ,P 1 ,P 2 ) 



(334) 



Thus through the new integration transformation, a new relationship between a phase space func- 
tion J-(rj, £) and its Weyl-Wigner correspondence operator F (Qi, Q2, Pi, P 2 ) is revealed. 
For example, from the following Weyl-Wigner correspondence 



(ef + iy 



■ exp 



4->| 2 + l£! 2 ) 
e J + 1 



exp{/[^ + + X_-l]}, 



(335) 



(iT + and X_ are defined in Eqs. d204T) and ( |2Q7b ) and d334l) we have 



(e/ + 1) 
x exp 



' /_1 (M a + iei a ) 



_ef + l 

= 2tt (^\F(Q 1 ,Q 2 ,P 1 ,P 2 )W) e i(^-^0. 
Using the Gaussian integration formula 



(336) 
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^ 2 £A ( £ - M )fa , -^ , )-fa- i /H£ , -M*V-A(hl 2 +|£| 2 ) 



1 



1 + A 2 



exp 



A(| M | 2 + |H 2 ) , A 2 ^*-/^) 



1 + A 2 



1 + A 2 



(337) 



in particular, when A = — itan (? — #), with v^— r . — - , 

" ' V4 2 / 1 A^ + l 2 tan a ' A^ + l 2 2 sin 



txt - ^tx, ttxt - I - ^tb>Eq. 033 becomes 



d £d r] (£ _ Mv *_ u * ) _ (n _ u)(e _ IJr * ) 



7T^ 

i sin a 
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i(\f,\ 2 + |zf) pW-tfv 



2 tan a 



2 sin a 



itan(J-|)(H 2 + |e| 2 ) 



(338) 



7r a > 

4 2 / 



|?y| 2 + |£| 2 )] represents an infinite long chirplet function. By taking / 



where exp [i tan 

£(| — a) in ( 13361 ), such that ie~ ta — , and comparing with Q338D we obtain 

(^1^(^1,^2,^1,^2)1^) 



2ir sin a 



exp 



*(ImI + W\ ) M* 7 * - M*^ 



2 tan a 



2 sin a 



(339) 



where the right-hand side of ( 13391 ) is just the CFrFT kernel whose properties can be seen in Ref. 
H87H . (One may compare the forms (13106 and ( 13391 ) to see their slight difference. For the relation 
between them we refer to Ref. H851[87l ). Dragoman has shown that the kernel of the CFrFT can be 
classically produced with rotated astigmatic optical systems that mimic the quantum entanglement. 
Therefore the new integration transformation (I321D can convert spherical wave to CFrFT kernel. We 
expect this transformation could be implemented by experimentalists. 



r- o 



11.4 Squeezing for the generalized scaled FrFT 



In some practical applications it is necessary to introduce input and output scale parameters [?, ?] 
into FrFT, i.e., scaled FrFT. The reason lies in that two facts: (1) the scaled FrFT may be more useful 
and convenient for optical information processing due to the scale parameters (free parameters) 
introduced into FrFT; (2) it can be reduced to the conventional FrFT under a given condition. In this 
subsection, by establishing the relation between the optical scaled FrFT and quantum mechanical 
squeezing-rotating operator transform in one-mode case, we employ the IWOP technique and the 
bipartite entangled state representation of two-mode squeezing operator to extend the scaled FrFT 
to more general cases, such as scaled complex FrFT and entangled scaled FrFT. The properties of 
scaled FrFTs can be seen more clearly from the viewpoint of representation transform in quantum 
mechanics. 

11.4.1 Quantum correspondence of the scaled FrFT 

The scaled FrFT H76H of ct-order is defined in a manner such that the usual FrFT is its special case, 

{! W ] . Jj^± r e x P {-, ^;+»> 2 + j-i} , w C340) 

y AitjivsYna J_ 00 ( ztana /ji/smaj 

where the exponential function is an integral kernel. In a similar way to deriving the quantum cor- 
respondence of FrFT in (12601) . and using the natural repression of single-mode squeezing operator 
Si in coordinate representation [ 88 J , 



V7* . 

we have 

x 2 /ri 2 + y 2 lv 2 

exp ' 



dx 

oo 



- ) (x\ , (341) 

fj, 



2 tan a [iv sin a 

= \f— 2-Ki[ive la sin a (y\ S[ (u) exp {zaa^a} Si (/Li) \x) , (342) 

which implies that the integral kernel in Eq.( |340D is just the matrix element of operator S\ (v) exp {iOa^a} Si (fi) 
in coordinate states. From Eq.( 1342D it then follows that 



T a [f (x)] = / dx (y\ SI (v) e* aa ' a Si (/i) \x) f (x) 

J — OO 

= {y\Sl(u)e iaaia S 1 ( f x)\f)=g(y), (343) 

which suggests 

\g)=S\(u)e^ a Si(p)\f). (344) 
From Eqs. Q343I ) and Q340I ) one can see that the scaled FrFT in Eq. d343| ) corresponds actually to 
the squeezing-rotating operator (^S\ (v) e l8a *' a Si transform in Eq.( |343D between two quantum 
states. 

11.4.2 The Scaled CFrFT 

On the basis of quantum mechanical version of one-mode scaled FrFT, we generalize it to two-mode 
case, i.e., we can introduce the integral 

T c a [/ ( V )} = (r/| Si (u) e -(4ai+4« 2 )s 2 I/) 

^ (rf\ Si (u) e H°-W+"W)s 2 fr) \rj) f (n) , (345) 



r* a 



where / (77) = (rj \ f). Using the natural expression of the two-mode squeezing operator S 2 (??), and 
noticing that (rf\ e l6 { a \ a ^+ a l a ^) |^) i s just the integral kernel of CFrFT (I3TU1) , we can reform J3451 ) as 

a ) fd 2 v 



2/«/sina J 7r 



xexp \mJ^±Km + lM^l±m I (346) 
2 tan a 2/i^sina 

It is obvious that Eq. ( 13461 ) is just a generalized CFrFT with squeezing parameters, we name it the 
scaled CFrFT. Thus we link a two-mode squeezing-rotating operator transform to the scaled CFrFT 
of complex functions. 



11.4.3 Entangled scaled FrFT 

On the other hand, recall that the entangled state \r/) can be Schmidt-decomposed as [89 



\v) = 



dx \x) x <g> 



r — V^rji ) ( 

2 



' \/2xTj 



(347) 



(348) 



we see that 

{x\,x' 2 \ rf) = e-<^6 (v^r/i +x' 2 - x'^j e iV ^ x '^, 
<r?| x 1 ,x 2 ) = e 4 " 1 "^ (V2 m +x 2 - a?i) e - *^ 1 " 2 . 
Using Eq. (131ft we have 

K (x' 1 ,x' 2 ,x 1 ,x 2 ) = S* (1/) e ia (4^+4^)s 2 ( p ) \ Xl ,x 2 ) . 

d -^- r?') (77'! (i/)e fa (°i<*+«J»0s a 0) |r?) (77 1^,^) , (349) 



where \x\,x%) — \x-y) ® \x 2 ) . On substituting Eqs. J347D and Q348D into Eq. Q349D , we can derive 

. A' 2 + A 2 . 



K(x 1 ,x 2 , xx,x 2 ) 




2tt sin a 
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2 tan a 



sin a 



2 tan a 



sma 



(350) 



where A, 



x l x 2 . K 

spiv ' ^ 



1 K u 



V2 



V2 • 



From Eq. d350D one can see that a new 2-dimensional (2D) scaled FrFT can be composed of one ID 
scaled FFT in its space domain and the other in its "frequency" domain, while the transform variables 
being the combination of two coordinates as shown in Eq.(??), so Eq. d350D is quite different from the 
direct product two ID scaled FrFTs that are both in 'space domain' are indicated in Eq.( |343D . Note 
that the new 2D scaled FFT is still characterized by only 3-parameter. Therefore, for any function 
/ (xi,x 2 ) — (xi,X2 I/) we can define an entangled scaled FrFT, i.e., 



(351) 



Fa[f{ x U x 2)\= / K (x[,x 2 ,xi,x 2 ) f (xi,x 2 )dxidx 2 

J — 00 

= (a/ x ,4| S\ (v) e H a ^+4^)s 2 ( M ) j/) . 



p<p 





a' ; 
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S W S 2 (rj') S 3 (?J) 

Figure 2: The Fresnel diffraction through ABCD optical system. 

Next we examine the properties of these scaled FrFTs in the quantum optics context. Without 
losing generality, for the additivity property, we consider the scaled CFrFT, 



TS +f3 [/ („)] =J^- W\ 4 H e ^)(4^+4^)s 2 ( M ) \V) I (V) 
Inserting the completeness relation of \rj) into Eq.( |352D yields 



(352) 



7T 



(rf\ St (i/) e H a W+4^)s 2 (r) S\ (r) e ^W+^) s 2 ( M ) |rj) / (ry) 



^ (r/| 0) e fa ( * ai+0 * 0a )s , 2 (/4t) IV) 
^ (VI 5| („) e «(«4ax+«4« a ) 5a ( T ) |^ r ) ^ [/ („)] = Ta T, [/ („)] 



(353) 



which is just the additivity property. It should be pointed out that the condition of additive operator 
for the scaled FrFTs is that the parameter v' of the prior cascade should be equal to the parameter 
(j,' of the next one, i.e., fx' = v' . For other scaled FrFTs, the properties can also be discussed in the 
similar way (according to their quantum versions). 

To this end, we should emphasize that different scaled FrFTs correspond to different quantum 
mechanical squeezing operators or representations. That is to say, it is possible that some other 
scaled FrFT can be presented by using different quantum mechanical squeezing operators or repre- 
sentations. 



12 Adaption of Collins diffraction formula and CFrFT 

The connection between the Fresnel diffraction in free space and the FrFT had been bridged by 
Pellat-Finet H90B who found that FrFTs are adapted to the mathematical expression of Fresnel diffrac- 
tion, just as the standard Fourier transform is adapted to Fraunhofer diffraction. In previous sections, 
a new formulation of the CFrFT and the Collins diffraction formula are respectively derived in the 
context of representation transform of quantum optics. In this section we inquire if the adaption 
problem of Collins diffraction formula to the CFrFT can also be tackled in the context of quantum 
optics. We shall treat this topic with the use of two-mode (3 parameters) squeezing operator and 
in the entangled state representation of continuous variables, in so doing the quantum mechanical 
version of associated theory of classical diffraction and classical CFrFT is obtained, which connects 
classical optics and quantum optics in this aspect. 

For Gaussian beam, the ABCD rule is equally derived via optical diffraction integral theory — the 
Collins integral formula. As shown in Fig. 2, if / (rj) represents the input field amplitude at point rj 



i r 



on Si, and g (77') denotes the diffraction field amplitude at point 7/ on S2, then Collins formula in 
complex form takes the form (??) . Next we shall examine adaption of the Collins formula to the 
CFrFT by virtue of the entangled state representation in quantum optics II S3 II . 



12.1 Adaption of the Collins formula to CFrFT 

Using the completeness relation of \r)) , we can further put Eq. fll90D into 



9(v') = W\U 2 (r,s)\f) = (r ] '\U2(r,s)fi 2 1 
and taking 77' = J A*i = as well as writing 



d 2 v 



ViV) (Mi»7 I/) 



(354) 



9{rf) 



>B a 
DK 



g) = G(a), f(»i V ) = F( v ) 



(355) 



where K and L are two constants to be determined later, then according to Eqs. ( 11901 ) and d354D 
we have 

cPri 



G(a)=^ 



(77'! U 2 (r,s)\ f i 1 f 1 )F( V ) 



1 



2iAL 2 



exp 



2K 2 



d 2 V 



■ exp 



i|?7| i (cr* 77 + arj*) 



2L 2 



2LKVAD 



F(rj). 



Comparing Eq. (I356D with Eq. J192D leads us to choose 



L 2 = tana, K = ^sm2a/ (2 AD). 



(356) 



(357) 



Then Eq. (13561) becomes 

G{a) 



cos a 



AD — cos 2 a . .2 
% — k 




F{rj) 



so Eq. Q358D is a standard CFrFT up to a quadratic phase term exp 
Eq.( 1357D and = rf, it can also be written as 



■ AD— cos a _ 2 
I : — 7. (T 

sm 2a 1 1 





AD — cos a |2 




' i 1 'i 2 " 


exp 


i : — 7. K 7 

sin 2a 1 1 _ 


= exp 









(358) 
, according to 

(359) 



which represents a quadratic approximation to a sphere wave diverging from a luminous point at 
distance 

R = An AB 2 (360) 
AD — cos^ a 

from 52. Let 5 be the sphere tangent to 52 with radius R (see Fig.2). A point on 5 is located by its 
projection on 52, this means that coordinates on 52 can also be used as coordinates on 5. Therefore, 
the quadratic phase term can be compensated if the output field is observed on 5 but 52- Then, after 
considering the phase compensation, the field transforms from Si to 5 is 

f*ns rv 

Gs (cr) = [F] . (361) 

In this way, the field amplitude on 5 is the perfect a — th FFT-C of the field amplitude on Si . 



12.2 Adaption of the additivity property of CFrFT to the Collins formula for 
two successive Fresnel diffractions 



The most important property of FrFT is that T a obeys the additivity rule, i.e., two successive FrFT 
of order a and j3 makes up the FFT of order a + p. For the CFrFT, its additivity property is proven in 
Eq. J314D . For Collins diffraction from Si to S' (see Fig.l), the additivity means that the diffraction 
pattern observed on S'(fj) (the sphere tangent to S3 with radius i?') and associated with F a+ fj 
should be the result of a first diffraction phenomenon (associated with J" Q ) on S (with 7/), fol- 
lowed by a second diffraction phenomenon (associated with J 7 ^) from S to S". This is a necessary 
consequence of the Huygens principle. Next we prove that such is indeed the case. 

Firstly, let us consider the field transform from Si (with 77) to S" (see Fig.2) described by the ray 
transfer matrix [M ,B\ C',D']. Similar to deriving Eq. ( 13601) , after the squeezing transform and the 
phase compensation, 

2MB' 

exp 



R' 



A'D' ~ cos 2 a' 



1 I -|2 

7 M 



R 



thus we can obtain the expression of CFrFT for Collins diffraction from Si to S' (not S3), 

cos a' 



G s > (1') 



A' 



(362) 



(363) 



where V = y w'W an< ^ Vi 



L' 2 = tana', K' = ^/ sin 2a'/ (2A'D'). 



(364) 



Eq. (l363D is the same in form as Eq. (13611) but with primed variables. Using Eqs. (13621) and Q190D one 
can prove that the transform from S 1 to S' is (see Eqs. Q190D , d358D -( |360l )) 



9S' in) 



\U 2 (r',s')\f)=G S ' {<?'). 



(365) 



In Eq.( 1365t we have taken the phase compensation term d362D into account. 

Secondly, let us consider the second diffraction from S to S3 determined by the ray transfer 
matrix [A", B", C", D"]. For this purpose, using the group multiplication rule of F 2 (r, s), we can 
decompose the diffraction from Si to S' into two parts: one is described as the matrix [A, B, C, D] 
from plane Si (with rj) to S 2 (S) (with 77O, the other is [A" , B" , C" , D"] from plane S 2 to S 3 (S') 
(with ff), then the total matrix from Si to S3 is 



A' 


B> N 




( A" 


B" \ 






a 


D' 


)- 


v C" 


D" I 


(c 





(366) 



Using Eq. d31D and the group multiplication rule of F 2 (r',s'), we can further put Eq. (13651) into 
another form 



Gs> W) = exp 



exp 



/4 2 



(fj\U 2 (r",s")U 2 (r,s)\f) 



d 2 a 



(r?| U 2 (r", s") ji 2 / — 1/4(7) (/4«7| U 2 (r, s) M ? / \^ v ) ( M W I/) 



d 2 rj 



d 2 a 



■ exp 



R' 



\U 2 {t" iS ")\^cj) 



^(n 2 a\U 2 (r,s)\»' 1 ri)F(r,) 



B 



d 2 a 



■ exp 



R' 



U 2 (r",s")\^ 2 a)G(a), 



(367) 



1- n 



where y! 2 = J j^-^ and we have made a reasonable assumption that fi\ = ni(so f (^[t]) = F (77)), 
which means that there are same scaled variants for the input field amplitudes on Si of the diffrac- 
tions from Si to S and from Si to S'. 

In order to examine the second diffraction domain from S to S' (not S3), we need to translate the 
output field amplitude G (a) observed on plane S2 to the field amplitude observed on sphere plane 
S, i.e., putting G (a) into Gs (c) (see Eq. (l361D ) by taking the phase compensation (see Eq. (l356l) ) 
into account. Thus the field transform from S to S' is (Gs (cr) — > Gs> (c')) 

B f d 2 a I i B\a\ 2 i\fj\ 2 \ „ „ , , 



B 1 f do f iB' / D" I 



CX P i T^7 1 777777 T77 I \® 



DK 2 2iB"J ir ^\K' 2 D'\2B" R 
iB\a\ 2 (A" 1 \ i(ao'* +a'o*) 



,2 



DK 2 \2B" RJ 2B n KIK ^ 



G s (a) . (368) 



Comparing Eq. J368D with Eq. J192D leads us to choose 



„„ T ., T . DD' B" /sin2a'sin2a 

svn.B = B"K'K\ = — \ — — , (369) 

Xl BB' 2 V A' ABB' ' 



and noticing that n' x = m yields A ,P' a , = Ai ? na , thus we have 



, B" cos a' sin a B" sin a' cos a 

^ = ~B sin/3 ' A= B^ sin/3 ' (37 ° j 



Combining Eqs. J354D and Q370D it then follows (letting a' = a + B) 

B 1 iB 1 ^ cosa' 



Dif 2 2iB" sin 2a zB" iA' cos a 2 sin /3 ' 

and 

B' { D" 1 \ 1 

= - cot /3, 



(371) 



if' 2 D' V2B" i?V 2 

H f A" 1 \ 1 

= - cot/3, (372) 



2 V 2B" RJ 2 



where we have used Eqs. (13601) . (13621) and (13661) . 

Substitution of Eqs.([36D, d369l , d37TT ) and ( 15721 ) into Eq. d368l yields 



G S . (a') = 




iA' 2 sin B 

<<*+P)F f a [ F ] (</) . (373) 



A' 

The first equation of Eq. d373D indicates that it is just a CFrFT of Gs (cr) from S to S'. Comparing 
Eq. (13731 with Eq. (13631 . we see 

T?T a [F] (a') = F a+f3 [F] (a 1 ) . (374) 

Thus we complete the study of adaption of CFrFT to the mathematical representation of Collins 
diffraction formula in quantum optics context. 



r - r\ 



13 The Fractional Radon transform 



Optical tomographic imaging techniques derive two-dimensional data from a three-dimensional ob- 
ject to obtain a slice image of the internal structure and thus have the ability to peer inside the 
object noninvasively. The mathematical method which complete this task is the Radon transforma- 
tion. Similarly, one can use the inverse Radon transformation to obtain the Wigner distribution by 
tomographic inversion of a set of measured probability distributions of the quadrature amplitude 
11911 19211 . Based on the Radon transform H93H and the FrFT we can introduce the conception of 
fractional Radon transformation (FRT) which combines both of them in a reasonable way. We notice 
the well-known fact that the usual Radon transform of a function / (r) can be proceeded in two 
successive steps, the first step is an n— dimensional ordinary Fourier transform, i.e. performing a 
usual FT of / (r) in n-dimensional k space, 



F[kj=F (te) = f(r) e'^^df, (375) 

where k = te, e is a unit vector, t is a real number. Its inverse is 

/ (r) = J F (tj e 2 ™ tF dk. (376) 
Letting s = t\ and rewriting Q375D as 

F (te) = J ds J drf (r) e - 27Tls S (s - k ■ f) = J d\ e - 2T,ltx J f (r) 5 (A - e • r) dr, (377) 
one can see that the integration over has been defined as a Radon transform of / (f), denoted as 

J f(r)5(X-e-r)df= f R (X,e). (378) 
So F (te) can be considered as a 1— dimensional Fourier transform of Jr (A, e) , 

F(te) = / d\e~ 2 ^ tx f R (\,e). (379) 



Its inverse transform is 

/>oo 

f R (X,e) = / F(te)e 2 ^ tx dt, (380) 



this ordinary 1— dimensional Fourier transform is considered as the second step. Combining result 
of J3751 ) and J380T ) we have 

/fl(A,e)= r [ f(r)e~ 2 * m F e 2 * ux drdt. (381) 

J — oo J 

i. e. two usual FTs make up a Radon transform, The inverse of J381D is 

f R (A, e) e 27Tllp e- 2nitx dkd\ = f (r) . (382) 



By analogy with these procedures we can make two successively FRFTs to realize the new fractional 
Radon transformation [?] . The ri-dimensional FrFT of / (r) is defined as 

$ a k [/] = (Ca)" / ex P -\l - / W & = F * , * = te. (383) 

a ' K J I 2 tan a sin a / 

r ia 1 1/2 

where a is named as the order of FrFT, C Q 



r r\ 



Firstly, we perform an 1— dimensional inverse fractional Fourier transform for F a (te) in i-space, 



fR,a (A, e) = [Cq] 1 n $- a ,t[F a (te)} 



[Ca] 1 n C- a l exp - 

J —CO \ 



i (A 2 + t 2 ) iXt 



2 tan a sin a 



F a (te) dt, 



(384) 



[Ca] 1 " was introduced for later's convenience. Then substituting (I383D into (I384D we have 



/fl,o: (A, e) = [C a ] n C- a (C a ) n I exp 
'i(r*-X 2 ) 



i(t 2 + X 2 ) i\t i^+t 2 ) ite- 



exp 



2 tan a 



2 tan a sin a 2 tan a sin a 
5 (X — e- f) f (r) dr, 



f (r) drdt 
(385) 



which completes the n-dimensional fractional Radon transformation. Especially, when a = tt/2, 
Q385D reduces to the usual Radon transform (I378D . Now we examine if the additive property of 
FrFT is consistent with (I385D . According to the additive property of FrFT 3 a $ a = $ a +p, and (13831) 
we see 



F a+P (k = te) = ff^ff^ [/] 
= (Cp) n (C, ^" 

z£ • f ik ■ £ 



i[r 2 +^ 2 



exp 



^ 1 fc 2 



2 tan a 2 tan /3 



x exp 



sin a sin /3 



~ $B+a,k ifi 



(386) 



The corresponding one-dimensional inverse FrFT should be 



27ri sin (a + /3) 

zA/i i/ii 



C-pC- a 



exp 



i (A 2 + /i 2 ) i (t 2 + /i 2 ) ' 
2 tan /3 2 tan a 



x exp 



sin /3 sin a 
,i(a+/3) 1 1 1 "")/ 2 



F a +f) (te) dtd^ji 



2-iri sin (a + f3) 

( i(r$-X 2 ) 

exp 



C-{ a +p) / F Q+/ s (te) exp 



i (t 2 + A 2 ) iA< 



2 tan (a + /3) 



2 tan (a + ^) sin (a + j3) 
8 (A - e • r) / (r) dr = f R>a+ p (A, e) , 



(387) 



which coincides with (I385D . From (I385D and Q387D , we can confirm that the transform kernel of 
ath FrFT is 



exp 



2 tan a 



5(X-e-r) 



(388) 



For example, one can calculate the fractional Radon transform of the n— mode Wigner operator 
to obtain some new quantum mechanical representations. Finally we give the inversion of the 
fractional Radon transformation, From (13851) we have 



• 2 \n/2 



(2ir sin 2 a 
which is an extension of J382D . 



fn, a (A, e) exp 



(A 2 - r 2 ) iXt ite ■ 



2 tan a sin a sin a 



dkdX = f(r), 



(389) 



r i 



In summary, based on the Radon transform and fractional Fourier transform we have naturally 
introduced the n-dimensional FRFT, in Ref. Zalevsky and Mendlovic [94] also defined 2-dimensional 
FRFT, but in different approach. We have identified the transform kernel for FrFT. The generalization 
to complex fractional Radon transformation is also possible 1951 . 



14 Wavelet transformation and the IWOP technique 

In recent years wavelet transforms H961 197H have been developed which can overcome some short- 
comings of the classical Fourier analysis and therefore has been widely used in Fourier optics and 
information science since 1980s. Here we present a quantum optical version of classical wavelet 
transform (WT) by virtue of the IWOP technique. 



14.1 Quantum optical version of classical WTs 



A wavelet has its energy concentrated in time to give a tool for the analysis of transient, nonsta- 
tionary or time-varying phenomena. (It is a wavelet because it is localized and it resembles a wave 
because it oscillates.) Mathematically, wavelets are defined by starting with a function ip of the real 
variable x, named a mother wavelet which is required to decrease rapidly to zero as \x\ tends to 
infinity, 



ip {x) dx = 0, 



(390) 



A more general requirement for a mother wavelet is to demanded ip (x) to have vanishing moments 
/f° x l ip (x) dx = 0, I = 0, 1, 2..., L. (A greater degree of smoothness than continuity also leads to 
vanishing moments for the mother wavelet) . The theory of wavelets is concerned with the represen- 
tation of a function in terms of a two-parameter family of dilates and translates of a fixed function. 
The mother wavelet ip generates the other wavelets of the family ip( p . s ), (^ is scaling parameter, s is 
a translation parameter, s G R), the dilated-translated function is defined as 



ip 



ii 



while the wavelet integral transform of a signal function / (x) e L 2 (R) by ip is defined by 

1 



WW (m, s) 



dx. 



We can express Q392D as 



WW O, s) 



U(fM,s)\f). 



(391) 



(392) 



(393) 



where (ip\ is the state vector corresponding to the given mother wavelet, |/) is the state to be 
transformed, and 



I dx 



(394) 



is the squeezing-translating operator H96[|981l99l , (x\ is the eigenvector of coordinate operator. In 
order to combine the wavelet transform with quantum states transform more tightly, using the IWOP 
technique we can directly perform the integral in Q394D (Q = (a + aP)/\/2, \x > 0) 



U(jl,8) 



7T/i J_ 
2[L 



dx: exp 



2^ 2 



XS i- r J 



V2xa~ — 



exp 



2(1 + ^ 2 



- + V2a^ + V2fia ) - V2-a) - -— r - Q 

fx 2\x l 



(395) 



r o 



This is the explicitly normal product form. Let p, = e , sechA = , tanhA = At 2+1 , using the 
operator identity e 9Q+a = : exp [(e 9 — 1) a^a] : , Eq. d395[ ) becomes 



t/ (/i, s) = exp 
x exp 
x exp 



2(1 + M 2) 



— tann A = seen A 

2 y/2 



a} a + — ] In scch A 



— tann A H = seen A 

2 V2 



(396) 



In particular, when s = 0, it reduces to the well-known squeezing operator, 



— \ (x| = ex P[^ (° 2 



For a review of the squeezed state theory we refer to (9j ■ 

14.2 The condition of mother wavelet in the context of quantum optics 



(397) 



Now we analyze the condition Q390D for mother wavelet from the point of view of quantum optics. 
Due to 



\x) dx=\p = 0), (398) 
where \p) is the momentum eigenstate, we can recast the condition into quantum mechanics as 

%ji (x) dx = -> (p = 0| ip) = 0, (399) 



which indicates that the probability of a measurement of \ip) by the projection operator \p) (p\ with 
value p = is zero. Without loss of generality, we suppose 



I M 



= G(at)|0)=X>»a tn |0>, 



(400) 



71=0 



where g n are such chosen as to letting obeying the condition Q390D . Using the coherent states' 
overcompleteness relation we have 



(p = 0| t) = (p = oi \z) (z\Y,g^" |0) 

(*)" 



e -\z\ z *n 



E 



(401) 



Eq. J401D provides a general formalism to find the qualified wavelets. For example, assuming g 2n = 
for n > 3, so the coefficients of the survived terms should satisfy 



and becomes 



.9o + .92 + 3.g4 + 15g 6 = 0, 

(ffo + .9 2 a f2 + 54fl f4 + .9ea t6 ) |0) 



(402) 
(403) 
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Figure 3 : Traditional Mexican hat wavelet. 



Projecting it onto the coordinate representation, we get the qualified wavelets 

$ (x) = Tr^^e-^' 2 [g + 92 (2x 2 - l) + g 4 (4a; 4 - 12s 2 + 3) 

+g 6 (8a; 6 - 60a; 4 + 90x 2 - 15)] , (404) 

where we have used (x \n) = (2 n n\y/%)~ 1 ^ 2 H n (x) er* 2 1 2 , and H n (x) is the Hermite polynomials. 
Now we take some examples. 

Case 1: in ( 1401 D by taking g a = |, g% = — 5, <?2n = (otherwise), we have 

l^)M = ^(l-« t2 )|0>, (405) 

it then follows 

V'm (x) = \ (x\ (1 - at 2 ) |0) = \ (x| (|0) - V2 |2)) = Tr" 1 ^"* 2 / 2 (l - x 2 ) , (406) 

which is just the Maxican hat wavelet, satisfying the condition e~ x2 / 2 (l — x 2 ) dx — 0. Hence 
i (l — a^ 2 ) |0) is the state vector corresponding to the Maxican hat mother wavelet (see Fig. 3). 
Once the state vector (ip\ corresponding to mother wavelet is known, for any state |/) the matrix 
element m U s) \ f) is just the wavelet transform of f(x) with respect to ■ 

Case 2: when go = —2, g 2 = — 1, 54 = 1 and g 6 = 0, from G404D we obtain (see Fig. 4) 

ip 2 (x) = 2ir- 1/A e~ x2/2 (2x 4 - 7x 2 + l) , (407) 

which also satisfies ip 2 {%) dx = 0. Note that when g = —1, g 2 = —2, g 4 = 1 and g e = 0, we 
obtain a slightly different wavelet (see Fig. 5) . Therefore, as long as the parameters g 2n conforms 
to condition Q402D , we can adjust their values to control the shape of the wavelets. 

Case 3: when g = 1, g 2 = 2, 34 = 4 and g 6 = —1, we get (see Fig. 6) 

^ 3 (a;) = 7r~ 1/4 e~ 2;2/2 (-8a; 6 + 76a; 4 - 134a; 2 + 26) , (408) 



r a 



Figure 4: Generalized Mexican hat wavelet tp 2 (x)when g — —2, g 2 = — 1, .94 = land g§ = 0. 
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Figure 5: Generalized Mexican hat wavelet when go = —1, g 2 = —2, g± = land g 6 = 0. 
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Figure 6: Generalized Mexican hat wavelet ip^ (.T)when g = 1, g 2 = 2, #4 = 4and g 6 = — 1. 



and J_ ^3 (x) cfa = 0. From these figures we observe that the number of the nodes of the curves at 
the a>axis is equal to the highest power of the wavelet functions. 

To further reveal the properties of the newly found wavelets, we compare the wavelet trans- 
form computed with the well-known Mexican hat wavelet tpx (x) and that with our new wavelet 
ip2 {x). Concretely, we map a simple cosine signal cos7ra; by performing the wavelet transforms with 
ijji [T (x — X)], i = 1,2, into a two-dimensional space (X,T), where X denotes the location of a 
wavelet and a its size. The resulting wavelet transforms by ipx (x) (=ipM (x)) and ip 2 (x) are 



where 2 / V3 and 1 /VSO are the normalization factors for tpi and i\) 2 respectively, the wavelet integral 
Qi (X, T) are also called wavelet coefficient which measures the variation of cos7ra; in a neighbor- 
hood of X, whose size is proportional to 1/T. The contour line representation of £l\ (X : T) and 
fl 2 (X : T) are depicted in Fig. 7 and Fig. 8, respectively, where the transverse axis is X-axis (time 
axis), while the longitudinal axis (T-axis) is the frequency axis. 

It is remarkable that although two overall shapes of the two contour lines look similar, there exist 
two notable differences between these two figures : 1) Along T-axis fli (X,T) has one maximum, 
while n 2 (X,T) has one main maximum and one subsidiary maximum ("two islands"), so when 
ip 2 scales its size people have one more chance to identify the frequency information of the cosine 
wave than using tp^. Interesting enough, the "two islands" of f2 2 (X, T) in Fig. 8 can be imagined 
as if they were produced while the figure of (X, T) deforms into two sub-structures along a- 
axis. 2) Near the maximum of fl 2 (X, T) the density of the contour lines along a-axis is higher than 
that of fli (X, T), which indicates that the new wavelet ij) 2 is more sensitive in detecting frequency 
information of the signal at this point. Therefore, t/> 2 (x) may be superior to ipM (x) in analyzing 
some signals. Finally, we mention that there exist some remarkable qualitative similarities between 
the mother wavelets presented in Figs. 3 through 6 and some of the amplitude envelopes of higher 
order laser spatial modes and spatial supermodes of phase locked diode laser arrays H1001I1011I102I1 . 
which are due to spatial coherence. 
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Figure 7: Contour line representation of Oi (X, a) . 
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Figure 8: Contour line representation of Vt2 o) • 
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14.3 Quantum mechanical version of Parseval theorem for WT 

In this subsection, we shall prove that the Parseval theorem of ID WT |f98l 1991 [10311 : 

% \ dsW f h (M, *) Wlfa (M, *) = / h (x) ft (x) dx, 



(411) 



where ip (x) is a mother wavelet whose Fourier transform is ip (p) ,C^ = 2n dp < oo. In the 

context of quantum mechanics, according to Eq. Q393D we see that the quantum mechanical version 
of Parseval theorem should be 



% r ds M U fa s) I/O (f 2 \ C/t fa s) = (/ a |/ x ) . 

A 1 J-oo 

and since if) (x) — (x\ so ijj (p) involved in is (p\ tp) , (p| is the momentum eigenvector 

1 



ip fa) = (p\ tp) 



'2tt J- 



dxiP (x) e~ lpx . 



(412) 



(413) 



Eq. J412D indicates that once the state vector (ip\ corresponding to mother wavelet is known, for any 
two states |/i) and |/ 2 ), their overlap up to the factor (determined by Eq. J426D ) is just their 
corresponding overlap of WTs in the (/x, s) parametric space. 

Proof of Equation Q412D : In order to show Eq.( 1412D , we calculate 



C/t fa s ) |p) = -± 



|m| j-oo 

-iEl 

e v 



dx \ x) 



i _ 



A* 

dx \x) e l ^ x 



p) 



which leads to 



dsU^ fa s) (jj| U fa s) = 2ttS (p - p') 
where we have used the formula 

— e-(p-p') =S(p-p'). 

, 27T 



P_\ /P 

A* / \ H 



(414) 



(415) 



(416) 



Inserting the completeness relation dp \p) (p\ = 1 into the left side of Eq. fl411D and then using 
Eq. J415D we have 



L.H.S of Eq.d4lTT) 



dfj, 



' — oo A* J — 

J — oo /* J —oo 

h + h, 



dsdpdp'r ( P ) i, fa) (h\ [/t ( M , s ) | P ') (p| t/ (m, S ) \h) 



dp%j)* (p) ip fa) (f 2 



where 



h = 2tt f°°% r d P r(p)^( P )(f2 

JO M ./-oo 

= 2^ r dp [ [ X |V (Ml 2 —1 (/a |P> (p| A) 

./-oo [Jo A* J 



p\/p 

A*/ \A* 



P\/P 
A*/ \A* 



A) 



(417) 



(418) 



and 



h = 2n j % ^ dp\^{p)\ 2 (f, 

J — oo A* J— oo 



2vr / dp 



P\/P 
A*/ \P 



2 ^P 



I^(-pp)I — 



/l) 

(/a Ip> (Pi fi) 



Further, we can put Eqs. 0418D and Q419D into the following forms, 



A* J-oo 



p\/p 

P/ \ A* 



/i 



IVKpp)!' 



= G*/ dp(h\p)(p\f 1 ) + 2ir dp 

JO J -oo 

p oo pQ 

= C^ dp(f 2 \p)(p\fi)+2n dp 

JO J—oo 
poo pO 

= cv / dp (h \ P ) ( P \ fx) + c; / dp (h \ P ) (p| /i) 



PP 



I^(P')| 2 ^' 
P . 



(/a |P) (Pi /i> 
(/ 2 IP) (Pi A) 



and 



where 



h = cU dp (f 2 \ P ) ( P \ h) + cv / dp (h \ P ) ( P \ h) 



Jo M Jo P 



|^(p')| 2 ^ = 27r / |r( -,,)| 



thus when the definite integration satisfies the admissible condition, i.e., 



dp 



P Jo 



dp 



h&oor-f = i^(-p)i- 



which leads to 



2vr |jV,(p)| 2 ^j=2C . 



Eq. 041 7D can be transformed to 

/oo 
dp (/a |p) (p| A) = R-H.S of Eq.lHTTT). 
-oo 



where 



C,i, — 2ir 



iv>(p)r 



■dp < oo, 



/q P 

thus the theorem is proved. Especially, when |/i) = |/ 2 ) , Eq. J412D becomes 



I — oo M" J—oo 

which is named isometry of energy. 



ds\(ip\U(fx,s)\f 1 )\ 2 = 2C^(f 1 |/0 



14.4 Inversion formula of WT 



Now we can directly derive the inversion formula of WT, i.e. we take (/2I = (x\ in Eq. (1412D , then 
using Eq.( |393D we see that Eq. J412D reduces to 



OO 7 />OG 



- 00 M J - :\_ 

Due to Eq. J394D we have 



% I dsW^h (p, s) (x\ (m, a) |V) = 2C^ (x |/0 



(x\tf(p,s) = -=(x\ / dx'\x') 



It then follows 



which means 



00 j />oo 



( i f I: 



oo r" </ — 00 



1 



dsW,pfi (fi, s) 
dp, 



1 x — s 



1 / X ~ s 



2CV i-00 P 2 \/|p[ 
this is the inversion formula of WT. 



ImI \ m 

dsip 



p vImI \ ^ 

V) = (x |/i) , 
W^/i (P> s) . 



(428) 



(429) 



(430) 



(431) 



14.5 New orthogonal property of mother wavelet in parameter space 

Form the Parserval theorem d411D of WT in quantum mechanics we can derive some new property 
of mother wavelet 1110411 . Taking = \x), |/ 2 ) = \x') in (141 ID one can see that 



dp 



, ds^(^-)r (^^-) =2C^5(x-x'), (432) 
00 P IpI J-00 \ P / V P / 

which is a new orthogonal property of mother wavelet in parameter space spanned by s). In a 
similar way, we take \fi) = |/ 2 ) = \n) , a number state, since (n \n) = 1, then we have 

/OO 1 />OG 
4/ ds|(^|C/(M,s)|«)| 2 = 2C^, (433) 
-00 P J —00 

or take = |/ 2 ) = |z) , \z) — exp ( — |z| 2 /2 + zat) |0) is the coherent state, then 



f°° ^ f°° ds \ y,\ [/ ( Mj s ) | z ) p = ac^. (434) 

./ — 00 P ./ —00 

This indicates that is -independent, which coincides with the expression in ( 14261 ) . Next, we 
consider a special example. When the mother wavelet is the Mexican hat d4Q6t , we have 



(p) ^ (P I^m) = \ ((P |0) - V2(p |2)) = TT-^Ve"^ 2 



(435) 



where 

(p |n) = j-!^ e - p2/2 g„ (p) • (436) 

Here if„ (p) is the single-variable Hermit polynomial [110511 ■ Substituting Eq. d435D into Eq. ( 14261 ) we 
have 

C f 2n r ^MLdp = v^F. (437) 
Jo P 



r-7/"l 



Thus, for the Mexican hat wavelet J406D , we see 



m 2 H 



= 2^8 (x ~ x') . 



(438) 



Eq. J438D can be checked as follows. Using Eq. ( 14061 ) and noticing that tpw (x) — tpM (—%), we can 
put the left hand side of Eq. fl438D into 

L.H.S.of (|438l) 

/>oo />oo 

= 2 / du dsipM {ux — s) ip* M (ux 1 — s) 

JO J-oo 



0, 



x 7^ x 1 



\ J °° du — > oo , x = x' 
where we have used the integration formulas 



R.H.S.of (|438|) . 



r°° 2 3 

/ (1-s 2 ) exp(-s 2 )ds = -v^F, 

J — oo ^ 



(439) 



(440) 



and 



e - S V2-(b- S ) 2 /2 dfi 



(441) 



Next, we examine if the Morlet wavelet obey the formalism G432D . The the Morlet wavelet is 
defined as lfl06l [T07l [T08H 



Substituting ( I442D into the left hand side of ( I432D yields 



/2 



(442) 



30 <fc 

oo ImI 



x - s , 



0. 



"1 2(l + e-/ 2 -2e-3/V4j ir ^^ 0O! 
Thus the Morlet wavelet satisfies Eq. Q432D . 



.x = x' 



(443) 



14.6 WT and Wigner-Husimi Distribution Function 

Phase space technique has been proved very useful in various branches of physics. Distribution func- 
tions in phase space have been a major topic in studying quantum mechanics and quantum statistics. 
Among various phase space distributions the Wigner function F w (q,p) H66II67B is the most popularly 
used, since its two marginal distributions lead to measuring probability density in coordinate space 
and momentum space, respectively. But the Wigner distribution function itself is not a probability 
distribution due to being both positive and negative. In spite of its some attractive formal properties, 
it needs to be improved. To overcome this inconvenience, the Husimi distribution function Fh (<?',??') 
is introduced [ 109], which is defined in a manner that guarantees it to be nonnegative. Its definition 
is smoothing out the Wigner function by averaging over a "coarse graining" function, 



Fh(q,p, k) 



dq'dp'F w (q',p')exp 



-K (q' - q) 



2 (p'-p) 



(444) 



where k > is the Gaussian spatial width parameter, which is free to be chosen and which deter- 
mines the relative resolution in p-space versus g-space. 

In the following, we shall employ the optical wavelet transformation to study the Husimi distribu- 
tion function, this is to say, we shall show that the Husimi distribution function of a quantum state 
can be obtained by making a WT of the Gaussian function er x / 2 , i.e., 



(V>| A h (q,p,K)\ij) = 



where 



dxip* 



(kq + ip), h = Vk, 



(445) 



(446) 



and (-01 Ah (q,p) \ip) is the Husimi distribution function as well as Ah {q,p, is the Husimi operator, 



A h {q,p,n) = — — : exp- 

1 ~\~ K, 



-k (q — Q) 2 {p-py 



1 + K 



1 + K 



(447) 



here : : denotes normal ordering; Q and P are the coordinate and the momentum operator. 

Proof of Eq. (|445D . According to Eqs. fl392D and Q393D , when |/) is the vacuum state |0), e^ 2 / 2 
7T 1 / 4 (x |0), we see that 



-1/4 



dx , „ / x — s 



e-^^dx^ (ip\U(^s)\0) 



From Eq. Q396D it then follows that 

U{n,s) |0) = scch 1/2 A exp 



-* 2 o f « ^ « t2 ^' 
?r- ■= seen A tann A 

L2(l + /i 2 ) V2 2 



10) 



Substituting Eq. (l446D and tanh A = cosh A = into Eq.( |449D yields 



(448) 



(449) 



p I ipq , 

e -2^ + 7J+T[/ = = -y^g - ip/VK) |0) 



2y^ 
1 + K 



1/2 



exp 



2 (1 + k) 2(1 + k) 

1- K 



+TTT ( K 9 + *P) + on , X 2 f 1°) = IP' 



1 + k v "' ' ' 2(l + «) 
then the WT of Eq. Q448D can be further expressed as 



e 2 K ^»+i 



fi.7' 



-oo (K7r) 



/./, 



-x 2 /2 



(450) 



(451) 



Using normally ordered form of the vacuum state projector |0) (0| =: e a a : ,and the IWOP 
method as well as Eq.(??) we have 



\p,q) KK (p,q\ = - exp 



-k (q — Q) 2 {p-Pf 

1 + K 1 + K 



= A/j {q,p,K) 



(452) 



Now we explain why A^ (q, p, n) is the Husimi operator. Using the formula for converting an operator 
A into its Weyl ordering form H110H 



A = 2 J ^-{-P\ A W) ;exp{2(/3*a-a t /3 + a t a)};, 
d 2 /3 = dfiidfo, P = /3i+ifr, 



(453) 



r-7 O 



where the symbol ; ; denotes the Weyl ordering, \ j3) is the usual coherent state, substituting Eq.( 1452D 
into Eq. J453D and performing the integration by virtue of the technique of integration within a Weyl 
ordered product of operators, we obtain 



\P,Q) KK (P,«l = 2 ' exp 



-K (q - Q) 



2 (P-Pf 



(454) 



This is the Weyl ordering form of \p, q) KK (p,q\ . Then according to Weyl quantization scheme H37H we 
know the classical corresponding function of a Weyl ordered operator is obtained by just replacing 
Q "/•/' •/''• 



exp 



n(q-Q) 



2 (P-Pf 



exp 



-«(« - q'Y 



(p-pT 



(455) 



and in this case the Weyl rule is expressed as 

\p, q) KK (p, q\ = 2 / dq'dp' '5(q'-Q)6 (p' - P) \ exp 

= 2 / dq'dp A w (q',p') exp -k (q - qf 



l /s2 (g ~ 

-n(q-q) 



(456) 

where at the last step we used the Weyl ordering form of the Wigner operator A w (q,p) H111H 

A w (q,p)=[S(q-Q)5(p-P)[. (457) 

In reference to Eq. J444D in which the relation between the Husimi function and the WF is shown, 
we know that the right-hand side of Eq. 04561 ) should be just the Husimi operator, i.e. 



\P,l) KK (p,q\ = 2 / dq'dp' A w (q',p')exp 
= &h(q,p, k), 



-k {q' - q) 



2 (p'-p) z 



(458) 



thus Eq. Q445D is proved by combining Eqs. d458D and J451D . 

Thus the optical WT can be used to study the Husimi distribution function in quantum optics 
phase space theory [112]. 



15 Complex Wavelet transformation in entangled state repre- 
sentations 

We now turn to 2-dimensional complex wavelet transform (CWT) H113II . 
15.1 CWT and the condition of Mother Wavelet 

Since wavelet family involves squeezing transform, we recall that the two-mode squeezing operator 
has a natural representation in the entangled state representation (ESR), cxp[A (a\a\ — a^J] — 

jt J^oo ^} P = &X ' t ^ us we are natu rally led to studying 2-dimensional CWT in ESR. Us- 

ing ESR we can derive some new results more conveniently than using the direct-product of two 
single-particle coordinate eigenstates. To be concrete, we impose the condition on qualified mother 
wavelets also in \rj) representation, 

^(l)=0, (459) 



where ip (i]) = {r]\ 4>) ■ Thus we see 



\rj) = exp{-44} |00> = |f = 0) , (460) 
and the condition J459D becomes 

(f = 0| V) = 0. (461) 

Without loss of generality, assuming 

oo 

|V>)= 53 K n , m a\ n al m |00) , (462) 

n,m— 

then using the two-mode coherent 1 2:12:2) state we can write ( 14611 ) as 

<£ = 0| V) = <£ = 0| / ^ Z2 |ziz 2 ) (ziz 2 | ^ /v n , m ai"4 m |00) 

n,m— 

= E g w ,m / ^ ^ 2 ^l^2 m exp [~|zi| 2 - |z 2 | 2 - z lZ2 ] 

n,m—0 

n 3 ra=0 
oo 

= 53n!*r„,„(-l) n = 0, (463) 



n=0 



this is the constraint on the coefficient K n . n in (I463D , i.e., the admissibility condition for \ip). Thus 
Eq. ( 14621 ) is in the form: 

OO 

=J2nlK n , n \n,n) . (464) 

n=0 

To derive the qualified mother wavelet ?A (?7) = (v\ "0) from |^), noticing Eq.( [316t and d464t we 
have 

oo 

4, fa) = e -H 2 /2 53 Kn , n H n , n (77* 5 7?) (-1)" 

oo 

= e-|"l 2 / 2 53n!^ n L„(|r;| 2 ), (465) 

n=0 

where L n (x) is the Laguerre polynomial. In this case, we may name the wavelet in Eq. ( 14651 ) as the 
Laguerre-Gaussian mother wavelets, analogous to the name of Laguerre-Gaussian modes in optical 
propagation. For example: 

(1) When taking if ,o = h, K\ t x — \, K, hn = for n > 2, so we see 

h% = ~(l + al4)|00>, (466) 



which differs from the direct-product state ^1 — a\ 2 ^j |0) 1 <g> ^1 — a} 2 ^j |0) 2 . It then follows from Eq. 
( 14651 ) that 

^ (i?) = \ ( V \ (|00) + |11» - e-^"l 2 {l - i \ V \ 2 }, (467) 

which differs from e~v x +a V 2 (l — a; 2 ) (l — y 2 ), the direct product of two ID Mexican hat wavelets 
(see also the difference between Figs. 9 and 10). 
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Figure 9: The Laguree-Gaussian mother wavelet ipi (77). 




Figure 10: 2D Mexican hat mother wavelet (Hermite Gaussian mother wavelet). 



Figure 11: Laguerre-Gaussian mother wavelet ^2 (v) 



(2) when A .o = 1, K\ t \ = 3, K 2 ,2 = 1, A nj „ = for n ^ 3, we have (see Fig. 11) 

4> 2 (<n)= (6-7|77| 2 + |r/| 4 )e-^"l 2 . (468) 

(3) when A' ,o = 1, ^1,1 = 1, ^2,2 = 3, A 3i3 = 3, AT n>n = for n ^ 4, the mother wavelet tp 3 (77) 
(see Fig. 12) reads 

V> 3 (r?) = fu-Sll^ + Wl^-H^c-il"! 3 . (469) 



From the figures we can see that as long as the coefficients K n<n satisfy condition J465D , we can 
construct arbitrary complex mother wavelet by adding or reducing the number of coefficients, or by 
adjusting the value of them. And since only K rhm (m = n) survive in all the coefficients, the mother 
wavelets obtained are all circularly symmetric on the complex plane. 

Moreover, the CWT of a signal function F (77) by * is defined by 

W^F (/j, k ) = IJ^Lf (r,) r ■ (470) 

Using the (77 1 representation we can treat it from the quantum mechanically, 

W4F (//, «) = - ( ^ (VI ^— -) (v\ F) = (*| U 2 ( Ml k) \F) , (471) 



/J, J 7T /X 

where 



U 2 (ji,K)= - / 



?? K \{v\,» = e\ (472) 



/' 

is the two-mode squeezing-displacing operator. Using the IWOP technique we can calculate its 
normally ordered form, 



U2 (/i, k) = sechA: expl^ajaJ, — ai<22^ tanh A + (scchA — 1) (a\ai + a\a2^ 
+ ~ {v*a\ - ca,i) scc h ^ + Y + ~2 ( K * ai - K ° 2 - ^ ' K ' 2 ) J' : ' 



(473) 



Figure 12: Laguerre-Gaussian mother wavelet ip 3 (77) 



When k = 0, it reduces to the usual normally ordered two-mode squeezing operator. Once the 
state vector m (^\ corresponding to mother wavelet is known, for any state \F) the matrix element 
m (^| U 2 (Mj k ) \F) is J ust tne wavelet transform of F(r}) with respect to m (^\ ■ Therefore, various 
quantum optical field states can then be analyzed by their wavelet transforms. 



15.2 Parseval Theorem in CWT 

In order to complete the CWT theory, we must ask if the corresponding Parseval theorem exists 
H114H . This is important since the inversion formula of CWT may appear as a lemma of this theorem. 
Noting that CWT involves two-mode squeezing transform, so the corresponding Parseval theorem 
differs from that of the direct-product of two ID wavelet transforms, too. 
Next let us prove the Parseval theorem for CWT, 

"If J ~^ W ^9i (M, «) W; 92 (fx, k) = C^J -^g* 2 (77) gi (r,) , (474) 
where k = k± + in 2 , and 

% 'Jo iei 

ip (£) is the complex Fourier transform of rp (r?) , ip (£) = (£| rp) = ^Jl (£| ^ ^| According to 
J471D and Q472D the quantum mechanical version of Parseval theorem should be 



CL=Al #|^(0| 2 . (475) 



I" % I — W U2 (m, «) (sal ^ (m, «) IV) = (32 \gi) ■ 

JO ^ J n 



(476) 



Eq. (l476D indicates that once the state vector {ip\ corresponding to mother wavelet is known, for 
any two states \gx) and \g 2 ), their overlap up to the factor C$ (determined by (I475D ) is just their 
corresponding overlap of CWTs in the {p, k) parametric space. 

Next we prove Eq. (14741) or (14761) . In the same procedure as the proof of Eq. 04121 ). We start with 

A 



calculating \j\ (ji, k) |£) . Using (02]) and Q472D . we have 



4 \ e i«l«2-S2Kl) ) (477-) 



it then follows 



d 2 K 



Ul<Ji,K)\?) (t\U 2 (jJl,K) 



JL f ^ e *[(ei-*o<«»+(&-&M 

[I 2 J 7T 

4t 



0W| 

M / \M 



M/ \/* 

Using the completeness of |£) and Q478D the left-hand side (LHS) of J476D can be reformed as 

LHS of Eq. (14751) 



(478) 



o M° 



rt, 



x(Z\U 3 (ji,K)\gi) (S2|C^(A»,«)|0 (C | V>) 
= /^{4 j ( 00 ^|^(^)| 2 }(. 92 10 « 



(479) 



where the integration value in {..} is actually £— independent. Noting that the mother wavelet ip (rj) 
in Eq. J465D is just the function of \r/\ , so ip (£) is also the function of |£| . In fact, using Eqs. d465D and 
d42l ), we have 



^(0=e-Val€l ^ n ,„ff„ >n (|£|,|£|) 

where we have used the integral formula 

d 2 z 



C,\z\ 2 +tz+ V z* 



i e -^,Rfi(C) <o. 



So we can rewrite d479D as 

lhs of (Szu = c; | ^ ( 52 io (e isi) = c; < 52 m 



where 



^=4r^i^(^)i a =4r^i^(oi s 



(480) 



(481) 



(482) 



(483) 



Then we have completed the proof of the Parseval theorem for CWT in 04761 ) . Here, we should 
emphasize that 0476D is not only different from the product of two ID WTs, but also different from 
the usual WT in 2D. 

When \g 2 ) = \rj) , by using d472t we see (rj\ Xj\ (//, k) \ip) = j^ip ( ZZ 7T £ ) > tnen substituting it into 
( 14761 ) yields 

9i (v) = ttt ( "T / — W i>9i (Mi k ) V> '' 



(484) 



which is just the inverse transform of the CWT. Especially, when \gx) = \g 2 ) , Eq. (1476D reduces to 



or 



% ! — MU 2 { t i,K)\g 1 )\ 2 = C^{g 1 \ 9l ), 

M J 7T 



(485) 



which is named isometry of energy. 
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15.3 Orthogonal property of mother wavelet in parameter space 

On the other hand, when \g±) = \ri) , \g 2 ) = \rj'), Eq.( l476D becomes 

which is a new orthogonal property of mother wavelet in parameter space spanned by (/j, n). In a 
similar way we take \gi) = \g 2 ) = \m, n) , a two-mode number state, since (to, n |to, n) — 1, then we 
have 

f°°% [ — MU 2 (n, K )\m,n)\ 2 = C^ (487) 
Jo M J 7r 

or take \g±) = \g 2 ) — \z±,z 2 ) , \z) = exp ^— |z| 2 /2 + za'j |0) is the coherent state, then 

"' / ' 1 d -mU 2 {^K)\ Zll z 2 )\ 2 = C' r (488) 



o M 

Next we examine a special example. When the mother wavelet is ipi (77) in ( 14671 ), using J42D we 
have V (0 = \ l£| 2 e"5l«l 2 , which leads to C; = J °° |£| 3 d\£\ = |. Thus for Vi (-?), we see 

2 jT ^ / ^ (^) (^) = ^ (V - ,') ■ (489) 

Eq. Q489D can be checked as follows. Using ( 14671 ) and the integral formula 

^ ' < uy 2 \ 



(x 2 +y 2 ) 4 

we can put the left-hand side (LHS) of ( 1489D into 



<* 4 -^y 2 + y% R e (^)>0, (490) 



IfflrfflB-f ^-^ , (49D 

where = |r/' — «| 2 , j/ 2 = |ry — k| 2 . 
When 77' — rj, x 2 — y 2 , 

™.,«=/-A 7 =frf?.» (492) 

On the other hand, when 77 7^ 77' and noticing that 

X 2 = (T][ - Klf + (772 - K 2 f , 

y 2 = (vi - «i) 2 + (% - « 2 ) 2 , (493) 



which leads to dx 2 dy 2 = 4 1 J| dKidK 2 , where J (x, y) = 
0491D reduces to 



ki -m K 2- m 



As a result of J4931 , 



LHsof(m = -4/ ^ Ag ^~jr v* -°» (494) 



ctedy xy[x 4 — 4x 2 y 2 + y 4 ) 
H [7j (x 2 +7/ 2 ) 

where we have noticed that J (a;, y) is the funtion of (x 2 ,y 2 ) . Thus we have 



LHS of ( 14891 = I °° ' 77 , V = RHS of J489T ). (495) 
1 0, rjf=r]. 



15.4 CWT and Entangled Husimi distribution 



Recalling that in Ref.|[ll5], the so-called entangled Husimi operator A^ (a, 7, k) has been intro- 
duced, which is endowed with definite physical meaning, and it is found that the two-mode squeezed 
coherent state |er, 7) representation of A^ (er, 7, k) , Ah (ct, 7, ft) = \<r, 7, k) (a, 7, k\. The entangled 
Husimi operator Ah (a, 7, k) and the entangled Husimi distribution Fh {a, 7, k) of quantum state \ip) 
are given by 



and 




respectively, where F w (cr',7') = (ip\ A w (cr',7') IV") is two-mode Wigner function, with A w (cr',7') 
being the two-mode Wigner operator. Thus we are naturally led to studying the entangled Husimi 
distribution function from the viewpoint of wavelet transformation. In this subsection, we shall 
extend the relation between wavelet transformation and Wigner-Husimi distribution function to the 
entangled case, that is to say, we employ the CWT to investigate the entangled Husimi distribution 
function (EHDF) by bridging the relation between CWT and EHDF. This is a convenient approach 
for calculating various entangled Husimi distribution functions of miscellaneous two-mode quantum 
states. 



15.4.1 CWT and its quantum mechanical version 



In Ref. ||113ll , the CWT has been proposed, i.e., the CWT of a complex signal function g (77) by ip is 
defined by 

77 — 2 



A* J TT 



(498) 



whose admissibility condition for mother wavelets, J ^"ip {rj) — 0, is examined in the entangled 
state representations {r}\ and a family of new mother wavelets (named the Laguerre-Gaussian 
wavelets) are found to match the CWT 111 1311 . In fact, by introducing the bipartite entangled state 
representation (rj = 771 + irj 2 \ ,we can treat Q497D quantum mechanically, 



W^g (p, z) 



— w\ — 

7T fl 



( V \ g) = (i>\U 2 (n ) z) 



(499) 



where z = z\ + iz 2 € C, < fi e R, g (77) = (r]\ g) and ip {rj) — (rj\ tp) are the wavefunction of state 
vector |ff) and the mother wavelet state vector \ip) in {rj\ representation, respectively, and 



TT ( \ - 1 f d ^ 

A* J n 



rj — z 



(500) 



is the two-mode squeezing-displacing operator. Noticing that the two-mode squeezing operator has 
its natural expression in (77 1 representation (1361) . which is different from the direct product of two 
single-mode squeezing (dilation) operators, and the two-mode squeezed state is simultaneously an 
entangled state, thus we can put Eq. (15001) into the following form, 



U 2 (/J, z) = S 2 (aO 2) (z) , 
where D {z) is a two-mode displacement operator, D {z) \rf) = \rj — z) and 

d 2 77 



(501) 



T>{z) 



\r] - z) (77I 



exp 



IZ\- 



V2 

Dx{-z/2)D 2 {z*/2) 



iz 2 - 



Qi + Q 2 



V2 



(502) 



n r\ 



It the follows the quantum mechanical version of CWT is 

W*g ( M , C) = M S 2 (/*) 2) (^) | 5 ) = M S 2 (/*) £>i i-z/2) D 2 (z*/2) \g) . (503) 

Eq. (l503D indicates that the CWT can be put into a matrix element in the (ri\ representation of the 
two-mode displacing and the two-mode squeezing operators in Eq. (??) between the mother wavelet 
state vector and the state vector \g) to be transformed. Once the state vector (tp\ as mother 
wavelet is chosen, for any state \g) the matrix element U 2 {fi, z) \g) is just the wavelet transform 
of g(rj) with respect to {ip\ . Therefore, various quantum optical field states can then be analyzed by 
their wavelet transforms. 



15.4.2 Relation between CWT and EHDF 

In the following we shall show that the EHDF of a quantum state \ip) can be obtained by making a 
complex wavelet transform of the Gaussian function e _ ' ?) ' / 2 , i.e., 



A h (ct,7,k) \ip) = e 



where fi = e x = 0c, z = z\ + iz 2 , and 



-h\ 2 



A 



-\vf/2 



rj — z 



COShA r * / * m 

Zi = — — 7 - 7 - k {a + a ) , 

1 + K 

i cosh A 

Z-2 = — — [7 + 7 + K (<T - (J )\ 



1 + K 

Ah (a, 7, k) is named the entangled Husimi operator by us, 



A h (<7,7,/c) 



4k 



~2 ■ exp 



(ai + <4 - 7) (a\ +a 2 - 7*) 



1 + K 



(504) 

(505) 
(506) 



(507) 



(ijj I Ah (u,7, ft) |^>) is the Husimi distribution function. 
Proof of Eq. (15041 ). 

When the state to be transformed is \g) = 1 00) (the two-mode vacuum state), by noticing that 
(77 1 00) = e~^'/ 2 , we can express Eq.( 1499D as 



1 r A p -\ n 



[l J 7T 



— ] =^|C/ 2 (/x,^)|00). 



(508) 



To combine the CWTs with transforms of quantum states more tightly and clearly, using the IWOP 
technique we can directly perform the integral in Eq. d500D 111 1611 



U 2 (n,z) — seen A exp 



— ; — \z\ 2 + alal tanh AH — ( z*al — zal ) scch A 

2(1 + /i 2 ) 1 1 12 2 V 2 11 



x exp 



(a\ai + a^a^ InsechA 



z ai - za 2 
exp 1 — ; — ; — - a\a 2 tanh A 



1 + [i? 



(509) 



where we have set fi = e , sechA = , tanh A = ^ 2+1 , and we have used the operator identity 
e ga a _ . ex p w e g _ ^ a t a j . _ j n particular, when z = 0, U 2 z = 0) becomes to the usual normally 



o 1 



ordered two-mode squeezing operator S 2 (n). From Eq.( |5Q9D it then follows that 



U 2 1 00) = sechAexp 

1 r 



Oi - iz 2 ) (zi + iz 2 ) + + 

h a[a 2 tanh A 



+ - 



2(l + /i 2 ) 

(zi — iz 2 ) a\ — [z\ + iz%) a\ sech A \ |00) 



(510) 



Substituting Eqs. (f505l) , fl506t and tanh A = cosh A = into Eq. (f510l) yields 



1 I 1 2 tr-y — -y u 

e --M -^u^ru 2 ((i,z 1 ,z 2 )\Q0) 



1 + K 



exp 



I7I + K \a\" K<7 + 7 + 7* — KCT* + t t K — 1 1 i™x 1 x 

±LL -, ^ + -Oi + 4 + 44 ? 00) = cr,7)„ . 

2(k + 1) 1 + k 1 1 + k 1 « + 1 J 7 1 ' ,/K 



then the CWT of Eq.( 1508D can be further expressed as 



(511) 



(512) 



Using normally ordered form of the vacuum state projector |00) (00 1 =: e a i ai a 2 a2 : , and the 
IWOP method as well as Eq. (151 ID we have 



l°">7> re « (o-,7l = 



4k 



K<X* + 7* 



exp 



K+ 1 1 + K 

At- 1 



1 + K 



KC7 + 7 7 - KC7 K — 1 / + + \ + + 

H a\ + — a 2 H a[al, + aia 2 - a{a x - a 2 a 2 

1 + K 1 + K K + 1V / 

(ax + 4-7) (4 + a 2 - 7*) 



(i + K.r 



exp 



1 + K 



k (ai 



1 ~ a 2 — a j ya[ — a 2 — a 

1 + K 



= A,, (ct,7,k) 



(513) 



Now we explain why (a, 7, k) is the entangled Husimi operator. Using the formula for converting 
an operator A into its Weyl ordering form H37H 



.4 = 4 



d 2 ad 2 S3 



(-a, -(3\ A |a,/3) ' exp{2 (a*ai - a\a + f3*a 2 - a\(3 + a\a\ + a\a 2 ^j}' , (514) 



where the symbol ; ; denotes the Weyl ordering, |/3) is the usual coherent state, substituting Eq. fl513D 
into Eq. J514D and performing the integration by virtue of the technique of integration within a Weyl 
ordered product of operators, we obtain 



16k 



(1 + kY 



d 2 ad 2 f3 



,t 



-a, — /J| : exp < — 



(ai +4-7) (4 + a 2 - 7*) 



1 + K 



k {ai — a 2 — j ya[ — a 2 — a 

1 + K 



exp{2 ^a*ai — a{a + /3*o,2 — a|/? + a\a\ + 4 a 2^} 



4 exp <j — k I ai — a 2 — er j \ a[ — a 2 — a* J - ( ai + a 2 — 7 ] ( a} + 02 — 7 



(515) 



n o 



where we have used the integral formula 

J ^ CX p(c|z| 2 +^ + ^*) =-ie-?,Re(C)<0. (516) 

Eq. (1515l) is the Weyl ordering form of \a, 7) (er, 7] . Then according to Weyl quantization scheme 
we know the Weyl ordering form of two-mode Wigner operator is given by 

A w (a, 7) = 6 (ai — a\ — crj S (a\ — a 2 — a*^j 6 (a± + a 2 - jj 5 (a\ + 02 — 7*) [ , (517) 

thus the classical corresponding function of a Weyl ordered operator is obtained by just replacing 

a\ — a\ — > a' , a\ + a\ — > 7', i.e., 



4 exp <y— k yai — a 2 — erj [a[ — a 2 — 0* J [ai + a 2 — 7 J yi\ + a 2 — 7 

^4cxp(-«| ( 7'-a| 2 -i|y-7| 2 l, (518) 



K 

and in this case the Weyl rule is expressed as 



I°»7)kk ( CT '7l = 4y d 2 a'd 2 y S (ax - a\ - <rj S (a\ - a 2 - er*) 5 (ai +a\ -7) 
x <5 (a| + a 2 - 7*) [ exp |-/t \a' - af - i I7' - 7| 2 | 

= 4| rf 2 C r'd 2 7'A w (cr',7')exp|~ K |cr'- C r| 2 -i| 7 '-7| 2 |. (519) 

In reference to Eq. d497D in which the relation between the entangled Husimi function and the two- 
mode Wigner function is shown, we know that the right-hand side of Eq. Q519D should be just the 
entangled Husimi operator, i.e. 

l cr >7) KK ( a il\ = 4 J d 2 a'd 2 "f'A w (<r',7')exp |-k|<t' - erf - ^ | 7 ' - 7| 2 | = &h (<t,7>«0 , 



(520) 

K J 

thus Eq. fl504D is proved by combining Eqs. (1520D and (I512D . 

Thus we have further extended the relation between wavelet transformation and Wigner-Husimi 
distribution function to the entangled case. That is to say, we prove that the entangled Husimi 
distribution function of a two-mode quantum state is just the modulus square of the complex 
wavelet transform of I 2 with i/j (ry) being the mother wavelet up to a Gaussian function, i.e., 

2 2 

(ip\ A h (a, 7, k) \il>) = e~~l 7 l" J -^^e - ! 1 '! / 2 ip* ((77 — z) . Thus is a convenient approach for 

calculating various entangled Husimi distribution functions of miscellaneous quantum states. 



16 Symplectic Wavelet transformation (SWT) 

In this section we shall generalize the usual wavelet transform to symplectic wavelet transformation 
(SWT) by using the coherent state representation H117H . 



16.1 Single-mode SWT 

First we are motivated to generalize the usual wavelet transform, which concerns about dilation, 
to optical Fresnel transform (we will explain this in detail in section below), i.e. we shall use the 
symplectic-transformed — translated versions of the mother wavelet 

A,s- K (z) = [s{z-K)-r (z* - K * )] (521) 



as a weighting function to synthesize the original complex signal / (z), 



W*f (r, s;k) = J ^±f (z) W* s;K (z) , (522) 
d 2 z = dxdy, z = x + iy, 

this is named the symplectic-transformed — translated wavelet transform. One can see that the 
mother wavelet ip generates the other wavelets of the family ip* [s (z — k) — r (z* — k*)] through 
a translating transform followed by a symplectic transform, (r, s are the symplectic transform pa- 
rameter, |s| 2 — |r| 2 = l, k is a translation parameter, s, r and kgC), this can be seen more clearly 
by writing the second transform in matrix form 

^ +M( *Z" k . ),M=( _% 'J ), (523) 



where M is a symplectic matrix satisfies M T JM = J, J = ^ ^ 



I 
J' 

Symplectic matrices in Hamiltonian dynamics correspond to canonical transformations and keep 
the Poisson bracket invariant, while in matrix optics they represent ray transfer matrices of optical 
instruments, such as lenses and fibers. 



16.1.1 Properties of symplectic-transformed — translated WT 

It is straightforward to evaluate this transform and its reciprocal transform when / (z) is the complex 
Fourier exponentials, / (z) = cxp (z/3* - z*f3) , (note that zf3* - z*(3 is pure imaginary): 

W^f = y/7s J ^ exp (z/3* - z*p) f* [s(z-K)-r {z* - k*)\ 

= ^ exp [(z + «)/?*- (z* + k*) p] r (sz - rz*) 

= ^ cxp [(sV + rz'* + k) P* - (sz'* +r*z' + «*) p] 4>* (z') 

= exp [up* - K*p] y/s~ J ^ exp \z' (s* P* - r* p) - z'* (sp - rp*)] ^* [z') 
= Vs~exp[ K p*-K*P}<i>(s*p*-r*P), (524) 
where $ is the complex Fourier transform of ip* . Then we form the adjoint operation 

W; (W^f) (z) = Vs^j^ {W+f) (r, s; k)^[s(z - «) - r (z* - k*)} 

= \s\ $ (s*P* - r*P) J ^ cxp [(k + z)p* - [k* + z*) p] $ [-an + rn*} 

= \s\ exp (z/3* - z*P) $ (s*P* - r*P) 

x J ^ exp [-«' [s*p* - r*P) + k'* {sp - rp*)} V> («') 

= |s| exp (zp* - z*P) |$ (s*p* - r*P)\ 2 , (525) 

from which we have 

f W^iWjJ) (z)d 2 s f \<S>( s *p*-r*P)\ 2 
J —t- —2 = exp(z^ - z P) J d z s (526) 



n a 



so we get the inversion formula 



/ d 2 s |$ {s*/3* - r*f3)\ I \s\ 



Eq. G527D leads us to impose the normalization 

J d 2 s\<P(s*(3* -r*f3)\ 2 /\s\ = 1, (528) 
in order to get the wavelet representation 

/ (z) = [ d 2 sW; (W^f) (z) I \s\ 2 . (529) 



Then we can have a form of Parseval's theorem for this new wavelet transform: 
Preposition: For any / and /' we have 

d 2 nd 2 s f d 2 z 



W*f (r, a; «) W 4 f* (r, a; «) = j — / (z) /'* (z) 



(530) 



Proof: Let us assume F((3) and F'(f3) are the complex Fourier transform of / (z) and /' (z) respec- 
tively, 

F(J3) =J^f (*) exp - z*/3) (531) 
recall the convolution theorem defined on complex Fourier transform, 

d 2 zf (a — z,a* — z*) f (z) 

j £ z J fl F (p) e i<*'-**)P-{«-*)P* J ^-F , (fi')expW-zP») 
d 2 pd 2 f3'F'(j3')F(j3)e a 'P- a P'6{p - p) S {(i* - &*) 
d 2 pF(/3)F'{(3)e a * l3 - al3 * (532) 

so from ([5321 and (152TT) , ([5221 we see that W^f (r, s;k) = J (z) ip* s . K (z) can be considered 
as a convolution in the form 

d 2 zf (z) ip* [s (z — k) — r (z* — k*)] 
= J d 2 pF{(3)<$>* (s(3-r/3*)exp(nP* -k*/3) (533) 
It then follows from (I532D that 



W^f (r, s; k) W^f* (r, s; «) d 2 K 
x F'* (/?')$' (a/3' -r/3'*)S(j3-j3')d(j3* -j3'*) 

= \s\ J ^-F(/3)F'*((3) |$ (s(3 - r/3*)\ 2 , (534) 



Therefore, using J528D we see that the further integration yields 

<fs r 

W^f (r, s; k) Wipf* (r, s; rc) d 2 n 



±l F (f3)F'*(P) J ^ |$(s/3-r/3*)| 2 



(535) 



which completes the proof. 

Theorem: From the Proposition J53QD we have 



W^f (r, s; k) Vv, s;k (z) 2 = f (z) , 



(536) 



that is, there exists an inversion formula for arbitrary function / (z). In fact, in Eq. (I522D when we 
take / (z) = <5 (z — z') , then 



(537) 



Substituting A537D into (I535D we obtain (15366 . 



16.1.2 Relation between W^/ (r, s; ft) and optical Fresnel transform 

Now we explain why the idea of W^f (r, s; ft) is originated from the optical Fresnel transform. We 
can visualize the symplectic-transformed — translated wavelet transform in the context of quantum 

mechanics, letting / (z) = (z\ /), (z\ is the coherent state, \z) = exp [zal — z*a\ = y ^* J V |0) * s 

the vacuum state in Fock space, then Eq. (I521D can be expressed as 

Wj,f (r, s;k) = ^ I —r [s{z-K)-r (z* - «*)] / (z) 



s —r 
—r* s* 



Z — ft 
Z* — ft* 



<*!/) 



(i/>\F 1 (r,a ) K)\f) 



where F^ r ' s ' K ^ is defined as 



Fx (r, s,k) — y/s [ \sz — rz*) (z + k\ 

J 7T 



(538) 



(539) 



and 



.To know the explicit form of F\ (r, s, ft), we employ the 



s —r \ I z 
-r* s* ) \ z* 

normal ordering of the vacuum projector |0) (0| =: cxp (— a^a) : and the IWOP technique to perform 
the integration in Q539D , which leads to 



Fl (r,s,ft) = exp(- re + — re* 2 + - — ) exp (-— a u - — 

y ' ' ; M 1 1 8s 8s* ; V 2s* 2s* V 



s| 3 + |r| 2 ot 



exp 



t M, 1 

a ' a + — In — 

2 1 s* 



(540) 



The transformation matrix element of Fi 0", s, re = 0) in the coordinate representation \x) is just the 
kernel of optical diffraction integration j92D (Fresnel transform), this explains our motivation to 
introduce W^f (r, s; k). In particular, when r* = r = sinh A, s = s* = cosh A, F\ (r, s, ft = 0) reduces 
to it the well-known single-mode squeezing operator exp[^ (a 2 — a^ 2 )] which corresponds to dilation 
in the usual WT. 



n r 



16.2 Entangled SWT 

In the above subsection, the mother wavelet is gained through a translating transform followed by 
a symplectic transform. This motivation arises from the consideration that symplectic transforms 
are more general than the dilated transform, and are useful in Fresnel transform of Fourier optics, 
e.g. ray transfer matrices of optical instruments, such as lenses and fibers in matrix optics, while in 
quantum optics symplectic transforms correspond to single-mode Fresnel operator (or generalized 
SU(1,1) squeezing operator). 

Recalling that in section 9 we have introduced the 2-mode entangled Fresnel operator which is a 
mapping of classical mixed transformation (z, z') — > (sz + rz'*,sz' + rz*) in 2-mode coherent state 
\z, z 1 ) representation onto quantum operator F 2 (r, s), thus we are naturally led to develop the SWT 
in (I522D to the so-called entangled SWT (ESWT) H118H for signals g (z, z') defined in two complex 
planes, 



W^g (r, s; k, k') = 



d 2 zd 2 z' 



g(z,z') 4>*,s;k,k' 0> z ') > 



here 



r,s\k.k' (^5 Z ) 



^ - s*cj)[s [z-k)+r (V* - k'*) ,s(z' ~ k') +r(z* - k*)} , 



(541) 



(542) 



is used as a weighting function to synthesize the signal g (z,z') regarding to two complex planes. One 
can see that the mother wavelet <j) generates the family </>* [s [z — k) + r (z'* — k'*) ,s (z' — k') + r (z* - 
through a translating transform followed by an entangled symplectic transform. We emphasize that 
this transform mixes the two complex planes, which is different from the tensor product of two 
independent transforms (z, z') -> [s (z - k) — r (z* - k*) , s (z' — k') — r (z 1 * — k'*)} given by J538D . 
The new symplectic transform can be seen more clearly by writing it in matrix form: 



/ z 

z* 
z' 
z'* 



k 
k* 
k' 
k'* 



( 



M 



z 
z' 
z'* 



k 

k* 
k' 
k'* 



M 







r 



s* 



s 
r 



r 
s 



(543) 



where M is symplectic satisfying M T JM = J, J = 



1 




, X is the 2x2 unit matrix. 



For Eq. J541D being qualified as a new wavelet transform we must prove that it possesses funda- 
mental properties of the usual wavelet transforms, such as the admissibility condition, the Parseval's 
theorem and the inversion formula. When g (z, z') is the complex Fourier exponential, 



gi (z, z') = exp (z/3* - z* (3 + z'j* - z'* 7 ) , 
according to ( |54H -( f542l ) we evaluate its ESWT 

•d 2 zd 2 z' 



(544) 



2 exp (z/3* - z*/3 + z'j* - z'*j) 4>* r ^ k , [z, z') 
sz + rz , sz + rz 



'za z 



x exp [(z + k) P* — (z* + k*) P + (z 1 + k') j* — (z 1 * + fc'*) 7] . (545) 
Making the integration variables transform sz + rz'* — > w, sz' + rz* — > w', Eq. Q545D becomes 

'■wd 2 w' 

— 5 — 9 {w,w ) 



W m = s exp [k(3* -k*/3 + k'j* - k'*j) 

x exp [w {s*/3* + r* 7 ) - w* (s/3 + r 7 *) + w' (s* 7 * + r*/3) - w'* (s 7 + r/3*)} , (546) 
the last integration is just the complex Fourier transform (CFT) of <jf , denoting it as <£>* , we have 

W+gx = s exp {kf3* - k* p + fc' 7 * - k'*j) $* (s*/3* + r* 7 , s* 7 * + r*/3) . (547) 



Then we form the adjoint operation of Q547D , 

w;(w m ){z,z>) 

= s* / / ^—{(W^g) (r, a; k, k')}</> [s {z - k) + r (z'* - k'*) ,s(z'~ k') + r (z* - k*)] 



2 - , / i"d 2 kd 2 k' 



\s\ z $* (s*(3* + r* 7 , s* 7 * + r* (3) / / g— exp (jfe/3* - fc*/3 + k'j* ~ k'*-y) 



7T 

x <j) [s (z - k) + r 0'* - k'*) ,s(z'- k') + r [z* - k*)] 
= \s\ 2 $* (s*/3* + r* 7 , s* 7 * + r*/3) / / ^—<p [sk - rk'*,-sk' - rk* 



x exp [(fc + z) p* - [k* +z*)p + {k' + z') 7* - (A; + z'*) 7] 

[ [ d 2 vd 2 v' 

= \s\ $* (s*/3* + r*7, 5*7* + r*/3) exp (z/3* ~ z* j3 + z'7* - z'*j) / / - — 4> (v, v') 



x exp [-v (s*/3* + r* 7 ) + v* (s/3 + r 7 *) - t/ (s* 7 * + r* (3) + v'* (s 7 + r/3*)} , (548) 

where the integration in the last line is just the CFT of 4> (comparing with ( 1546t ), thus J548D leads 
to 

W; {W^gi) (z, z') = \s\ 2 exp (z/3* - z*/3 + z'7* - z'* 7 ) |$ {s*/3* + r* 7 , s* 7 * + r*/3)| 2 . (549) 
From Eq. d549D we have 

= exp (z£* - z*/3 + z' 7 * - z'* 7 ) 

x / d 2 s |$ (s*/3* + r* 7 , s* 7 * + r*/3)| 2 / |s| 2 , (550) 



ffi (z,z ) = — i3i ^ ^ ^ y — : r — — : oNl2 2 - (551) 



which together with (1544D lead to 

J d 2 s |$ (s*/3* + r*7, s*7* + r*^)| 2 / |s| 
Eq. (I551D implies that we should impose the normalization 

J d 2 s |$ (s*/3* + r* 7 , s* 7 * + r*/3)| 2 / |s| 2 = 1, (552) 

such that the reproducing process exists 



g x (z, z') = J d 2 sW; {W m ) (z, z') I ■ (553) 

Q552D may be named the generalized admissibility condition. Now we can have the corresponding 
Parseval theorem: For any g and g' we have 

d 2 kd 2 k'd 2 s f f 2 ,2 / / r\ 1* 



W^g (r, s; k, k') W; 9 ' (r, s; k, k') — A = / / d 2 zd 2 z'g (z, z') g'* (z, z') . (554) 



Proof: Assuming F (/3, 7) and F' ((3, 7) be CFT of g (z, z') and g' (z, z'), respectively, 

F (/3, 7) = / / ^g (z, z 1 ) exp (zp* - z*f3 + z' 7 * - z'* 7 ) , (555) 



recalling the corresponding convolution theorem 



i2 j2 I i * * I II* l*\ I I l\ 

a za z g (a — z,a — z ;a — z ,a — z ) g (z,z ) 

d 2 /3d 2 jF (/3, 7) F' (/3, 7) exp [a* ft - a(3* + a'* 7 - c/7* ) , 



(556) 



so from Eqs. J54TT ) and (I555l) - (l556l) we see that W^g (r, s; k, k') = J J d ^f z ' g (z, z') 4>* >s . kjk i (z, z') 
can be considered as a convolution in the form (noting that the CFT of cj)* is <&*, see ( |546| )-( f547l )) 



d 2 zd 2 z'g (z, z 1 ) (f>* [s(z~k) + r (z'* - k'*) ,s(z'- k!) + r (z* - k*)} 

d 2 j3d 2 1 F (/3, 7) $* (s*P* + r* 7 , s* 7 * + r*/3) exp (k/3* - k* [3 + fc'7* - fc'* 7 ) . 
Using Eq. J557D we calculate 

W<j,g (r, s; k, k') W%g' (r, s; k, k') d 2 kd 2 k! 

= \s\ 2 JJJJ d 2 pd 2 1 d 2 p'd 2 1 'F (/3, 7) $* (s*p* + r* 7 , s*7* + r*/3) 
x F'* (/?', 7') $ (s*/T + r* 7 , s* 7 * + r*/3) 5 (0 - /3') (5 (/?* - /?'*) 5 (7 - 7') S (7* - 7'*) 



(557) 



= H 2 J J d 2 f3d 2 1 F (0, 7) F'* (p, 7) |$ (s*p* + r* 7 , S * 7 * + r*/3)| 2 . 
As a consequence of (I552D and J558D the further integration yields 

J J^JJ W *9 (r, *; A, fc') (r, s; fc, fc') d 2 fcd 2 £;' 

d 2 pd 2 1 F(f3,- f )F'* ((3,j) J d 2 s\<S>(s*(3* +r* 7 , s* 7 *+r*/3)| 2 
d 2 /M 2 7 F 08, 7) F* 08, 7) = // rf 2 ^d 2 z' 5 (z, z') g'* (z, z') , 



(558) 



which completes the proof. 

Inversion Formula: From Eq. Q554D we have 



9 0, = / / / W^g (r, s; k, k') 4> r , S ;k,k' 0, z') 



d 2 kd 2 k'd 2 s 



(559) 



(560) 



that is, there exists an inversion formula for g (z, z') which represents the original signal g (z, z') 
as a superposition of wavelet functions (j> rjS ;k,k' , with the value of entangled wavelet transform 
W<j,g (r, s; k, k') serving as coefficients. In fact, in Eq. Q541D when we take 

g (z, z') =5(z-u)5(z* - u*) 6 (z' - v!) 5 (z'* - u'*) , (561) 

then 

' ' (562) 



W^g (r, s; k, k') = —<j)*,s;k,k' ( M > u ') 



Substituting (I561I) - (I562I) into Q559D , we obtain (15601) . We can visualize the ESWT in the context of 
quantum mechanics, letting g (z, z') = (z, z'\ g) and using Eqs. d542D - j543l ). Eq. J541D is expressed 
as 



W^g (r, s; k, k') = s 



s(z-k)+r (z'* - k'*) ,s(z' - k') +r(z* - k*)} g (z, z') 

M Z * z ,J* \(z,z'\g) = ^\F 2 (r, S ;k,k l )\g) (563) 
v z'* -k'* J 



nn 



where F 2 (r, s; k, k') is defined as 



F 2 (r,s;k,k') = sjj \sz + rz'* , sz' + rz*) (z + k,z' + k'\ 

\sz + rz'* , sz' + rz*) = \sz + rz'*) 1 <Ei \sz' + rz*) 2 . 



(564) 



When k = and k' = 0, F 2 (r, s; k = kl = 0) is just the 2-mode Fresnel operator. 

Thus, we have extended the SWT of signals in one complex plane to ESWT of signals defined in 
two complex planes, the latter is not the tensor product of two independent SWTs, this generalization 
is inevitable, since it resembles the extending from the single-mode squeezing transform (or Fresnel 
operator) to the two-mode squeezing transform (or entangled Fresnel operator) in quantum optics. 



16.3 Symplectic-dilation mixed WT 



Next we shall introduce a new kind of WT, i.e., symplectic-dilation mixed WT H119H . Recalling that 
in Ref. 1112011 we have constructed a new entangled-coherent state representation (ECSR) \a,x), 



| a, x) — exp 



1 2 1 , ,2 / a \ t 
-x - - \a\ + (x + -)a[ 



+(z-f)4- J(4 + 4) 2 



2 
|00) 



(565) 



which is the common eigenvector of the operator (Xi + X 2 ) /2 and a\ — a 2 , i.e., (a\ — a 2 ) \a, x) = 
a\a,x) and \{X\ + X 2 ) \a,x) = -^x\a,x) , where Xi — ^( a i + a \) is the coordinate operator, 
= 1, 2). \a, x) is complete, 



00 dx fd 2 a 



- , \<*,x) (a,x\ = 1, 
and exhibits partly non-orthogonal property (for a) and orthonormal property (for a;), 

(a', x' \a, x) 

1 i 1 

5 (x' — x) , 



(566) 



'7T exp 



-\{\a\ 2 + \a'\ 2 ) + \aa'* 



(567) 



so |a, x) possess behavior of both the coherent state and the entangled state. An interesting question 
is: Can we introduce a new kind of continuous WT for which the \a,x) representation underlies? 
The answer is affirmative. Our motivation of this issue comes from the mixed lens-Fresnel transform 
in classical optics 1112111 (see (12501) below) . 

By synthesizing (1392D and (15221) and in reference to (1566D we propose the mixed WT for g (a, x) 

(a = ai + ia 2 ): 



W^g (s, r, k; a, b) 



dx f d z a 



g(a,x) V£ r ,« ;a ,& («>^) 



(568) 



where d 2 a = daida2, the family of mother wavelet tp involves both the the symplectic transform of 
a and the dilation-transform of x, 



i>s,r,K-a.,b (01,X) 



s (a — k) — r (a* —«*),■ 



- b 



(569) 



Letting g (a, x) = (a, x\ g) , then J568D can be expressed as quantum mechanical version 

W^g (s, r, k; a, b) = (ip\ U (s, r, n; a, b) \g) , (570) 



r\r\ 



where U (s, 7% Ki« cL, b) is defined as 

U (s, r, ft; a, 6) 



dx /" d^a 
|a| 7_oo n/tt 7 2?r 

x — b 



sa — ra , 



(a + k, x\ . 



(571) 



U (s, r, k = 0; a, b = 0) is just the generalized squeezing operator, which causes a lens-Fresnel mixed 
transform. 

For Eq. 05681 ) being qualified as a new WT we must prove that it possesses fundamental prop- 
erties of the usual WTs, such as the admissibility condition, the Parseval theorem and the inversion 
formula. It is straightforward to evaluate the transform ( 15681 ) and its reciprocal transform when 

g (a, x) is the exponential gi (a, x) — exp (a*f3 — a(3* — ipx) , 



Wm = J±_ e K*f>-KP*-i P b I" da [da 



x ip* (sa — ra* , — )e° 
a 



J 2tt 



(572) 



Making the integration variables transform sa — ra* — > w, - — > x', leading to d 2 a — > d 2 w and 
dx -> |a| dx', (15721 ) becomes 



W^gi = ^/s~\^<t>* (s*(3* -r*p, ap) e "'f , -«r-if 
where $* is just the Fourier transform of tp* , 



(573) 



$* (s*/3* - r *p, ap) 



dx f d z w 



x e 



-co V n J 2n 

(sj3— r/3* ) — w(s* 0* — r* /3) — iapa; 



(574) 



Then we perform the adjoint WT of ( 1568D , using Q569D and (|573T ) we see 

W;(W m )(a,x) 



db { d 2 K 



0T J 2lT 



-W^gi 



x -0 



s{a — k) — r {a* — K*) , 



x — b 



-7= 

-co V 7 *" 



dV 

~2T 



erf ( SK < - rn'* , 6') 



,i|a| 5l (a,x)|$(s*/3*-r*/3, ap)| 5 



(575) 



From Eq. Q575D we obtain 



/ ^ -,w;(w m )(a,x) 

J-00 a J \s\ 

= 9l (a,x) r n / ttI^^'-^ a ^ 

l a l J |s| 



which leads to 



9i (a,x) 



I- 00 ^I^WUW^g 1 )(a,x) 



(576) 



(577) 



Eq. Q577D implies that we should impose the normalization 
such that the reproducing process exists 



(578) 



da f d s 

9l (a,x)= I W ;(W m )(a,x). (579) 

) a J s 



J578D may be named the generalized admissibility condition. Now we can have the corresponding 
Parseval theorem: For any g and g' we have 



dad& f d 2 Kd 2 s TTr , , ,\ 

— — W^g (s, r, k; a, b) W^g (s, r, k; a, 6) 



a" 



dx / d 2 ag (a, a;) g 1 * (a, x) . 



(580) 



Proof: Assuming F (/3,p) and F' (0,p) be the Fourier transforms of g (a, x) and g' (a, x), respectively, 



F(fi,p) 



dx f d 2 a 



>2lt J 7T 



g{a,x)e 



a (3* — a* {3-\-ipx 



(581) 



In order to prove J580D , we first calculate W^g (s, r, k; a, b). In similar to deriving Eq. dl84D , using 
(??), Q180D and the inversion formula of J581D we have 



W^g(s,r, ft; a, b) 
s 



30 



dp /• d 2 ^ 



x e 



a* /3—aj3* — ±px 



's a 



dp 



s (a — k) — r (a* —«;*), 



-F(/3,p) 



dx 


/■d 2 a 




/ 2tt 




x — b 


a 



, \/27r 7 7T 

x $* (s*/3* - r*£, ap) e «*/3-«/3*-ip6. 



(582) 



It then follows 



d& / d 2 nW^,g (s, r, «; a, 6) W^g' (s, r, k; a, 6) 

OO t/ 

|as| f°° dpdp' / d 2 /3d 2 /3'^(/3,p)F'*(/3',p') 



x $* (s*£* - r*/3, ap) $ (s*/T* - r*/3', ap') 



db f d 2 n 



2TT / 7T 



K*(0-^')-«(^*-/3'*)+i(p'-p)l) 



as 



dp / d 2 /3F(/3,p)F'* (/3,p)|$( s */3*-r*/?,ap)| 2 . 
Substituting ( [5831 into the left-hand side (LHS) of J5801 ) and using ( f578l we see 
LHS of (12341) = / dp / d 2 /3F (/3,p) F'* (/3,p) 

J — OO J 



(583) 



^ / ii|<r( S */r-r*/?,ap)| 2 
a ./ s 



dp / d 2 /3F(/3,p)F'* (/3,p). 



(584) 



Thus we complete the proof of Eq.( 1580D . 
Inversion Formula: From Eq. Q580D we have 

dadfe f d 2 Kd 2 s 



(a,x) 



2tt Isl 



W^g (s, r, k; a, b) ^ s , r ,K-,a.,b (a, x) 



(585) 



that is the inversion formula for the original signal g (a, x) expressed by a superposition of wavelet 
functions ^ s ,r,K;&,b (a, x) , with the value of continuous WT W^g (s, r, k; a, b) serving as coefficients. 
In fact, in Eq. ( 15681 ) when we take g (a, x) = 5 (a — a') 5 (a* — a'*) 5 (x — x') , then 



W^g(s,r, K;a, 6) 



1 



(586) 



Substituting ( 15861 ) into J5841 ) yields J5851 ). 

We can visualize the new WT W^g (s, r, k; a, b) in the context of quantum optics. Noticing that 
the generalized squeezing operator U {s,r,n = 0;a, b = 0) in J571D is an image of the combined 
mapping of the classical real dilation transform x — >• x/a, (a > 0) and the classical complex symplectic 
transform (a, a*) — > (sa — ra* ,s*a* — r*a) in \a,x) representation, one can use the technique of 
integration within normal product of operators to perform the integration in J571D to derive its 
explicit form (see Eq. (15) in Ref. [?]). The transform matrix element of U (s, r, n = 0; a, b = 0) in 
the entangled state representation | -77) is 



(r 7 |C/(s,r,K = 0;a,6 = 0)|j7 / ) 
da; 



d 2 a 



(V 



sa — ra 



(a,x\ i]') 



(587) 



In Fock space \rj = 771 + 1772) is two-mode EPR entangled state in 
Then using Q565I ) and ([29]), we obtain 



(77 \a,x) 



V2 



exp 



-A 



Substituting A588D into (I587D and using ([85]), we obtain 

(r ) \U(s,r,K = 0;&,b = 0)\r l ') 



-7= & W2 ~ V2/&) 



1 



x exp 



j]\a — \rj 2 x 



(588) 



^ {Atf - 2 VlV [ + Dr,l] 



(589) 



which is just the kernel of a mixed lens— Fresnel transform, i.e., the variable r]i of the object expe- 
riences a generalized Fresnel transform, while 772 undergoes a lens transformation. Thus, based on 
I a, x) we have introduced SDWT which involves both the real variable dilation-transform and com- 
plex variable symplectic transform, corresponding to the lens-Fresnel mixed transform in classical 
optics. 



17 Fresnel-Hadamard combinatorial transformation 

In the theoretical study of quantum computer, of great importance is the Hadamard transform. This 
operation is n Hadamard gates acting in parallel on n qubits. The Hadamard transform produces an 
equal superposition of all computational basis states. From the point of view of Deutsch-Jozsa quan- 
tum algorithm, the Hadamard transform is an example of the N = 2" quantum Fourier transform, 
which can be expressed as [122] 

|j) = i^ e2 " Jfc/2 " |fc> - (590) 
v 1 k=0 



no 



Now the continuous Hadamard transform, used to go from the coordinate basis \x) to the momentum 
basis, is defined as[ 123H 

1 f°° 

$\z) = -=- \ dy exp (2ixy/a 2 ) \y) , (591) 

where a is the scale length. $ is named Hadamard operator. Using the completeness of dx\x) (x\ ■■ 
1, we have 

3 = i / / dxd y ex P C 2 **!//" 2 ) ll/> (*l ■ (592) 

The above two transforms (Fresnel transform and Hadmard transform) are independent of each 
other, an interesting question thus naturally arises: can we combine the two transforms together? To 
put it in another way, can we construct a combinatorial operator which play the role of both Fresnel 
transform and Hadmard transform for two independent optical modes? The answer is affirmative, 
in this section we try to construct so-called Fresnel-Hadmard combinatorial transform. 



17.1 The Hadamard-Fresnel combinatorial operator 



Based on the coherent-entangled representation \a,x), and enlightened by Eq. ( 1871 ) and J591D we 
now construct the following ket-bra integration H124H 



U = 



_ [ [[dxdyex-p(2ixy/a 2 )\sa — ra*,y)(a,x\, (593) 

Vttct JttJJ 



we name U the Hadamard-Fresnel combinatorial operator. 

Substituting Eq.( |591D into Eq.( |593t , and using the two-mode vacuum projector's normally or- 
dered form 1 00) (00| = : exp [— afai — a 2 ] ■ as well as the IWOP technique we get 



U 



'ira 



J^A(z,z*) J J dxdyB{x,y)e C : 



(594) 



where 



(a+ + a+) + (oi + a 2 ) 2 
C = a{ai — a2, 



B(x, y) = exp 



y 2 + x 2 , , _i_ \ 2ixy 
h y (a{ + aj j + x (ai + a 2 ) H 5- 



and 



A (z, z*) = exp 



\sz — rz*| 2 + \z\ 2 sz — rz* f i i\ z* (ai — a 2 ) 



( a i - 4) 



they are all within the normal ordering symbol :: . Now performing the integration over dxdy within 
:: and remembering that all creation operators are commute with all annihilation operators (the 
essence of the IWOP technique) so that they can be considered c-number during the integration, we 
can finally obtain 



U 



1 4y/TTO 



s* Vc 4 + 4 
1 



1 r / a\ — a\ 



cr 4 - 4 ( a\+a\ 



1 



a\ — a\a\ — a 2 



1 r* / ai — a 2 
2^ I, V2 



V2 V2 

2 _4 



V2 
Aia 2 



2 (a 4 + 4) ^ V2 

a\ + a\ a% + a 2 
V2 V2 



a — 4 / a\ + a 2 



2 (a 4 + 4) V 



(595) 



which is the normally ordered form of Hadamard-Fresnel combinatorial operator. 



r\ a 



17.2 The properties of Hadamard-Fresnel operator 

Note 

a\ — a 2 a\ + a\ _ 

"1 



V2 ' s/2 



and 



ai — a 2 a\ — a\ 



V2 ' v/2 



ai + a 2 a\+ a 2 



(596) 



(597) 



ai ^ 2 can be considered a mode independent of another mode ai ^ 2 , thus we have the operator 
identity 

A 



exp 



/ (4 ±4) (°i ±fl 2) 



= : exp[^ (e 2 f - l) (4 ± a\ ) {a x ± a 2 )] : 



Using J598D we can rewrite Eq. fl597D as 



(598) 



u = u 2 u 1 = u 1 u 2 , 



(599) 



where 



(7i = , exp 



4 + a\ 



2 (a 4 + 4) I 



exp 



a! + al ai + a 2 , 4icr 2 
In ■ 



V2 V2 ((T 4 + 4) 



exp 



cr 4 - 4 / a x + a 2 



2 (a 4 + 4) V V2 



(600) 



and 



U 2 = exp 



1 a 2 



r I a\ — a. 
2? 



exp 



a} - flj ai - fl2 1 \ 1 



exp 



_r_ I o-i — a 2 



(601) 



while £/ 2 is the Fresnel operator for mode ai ^ 2 , U\ is named the Hadamard operator for mode 

ai+Q2 
%/2 • 

It then follows 



V2 



V2 ■ 2 



V2 



V2 ' 



V2 V2 



ai — a 2 at — at 
t + s- 



72 



V2 ' 



(602) 



from which we see the Hadamard-Fresnel combinatorial operator can play the role of Fresnel trans- 
mation for ai ^ 2 ■ Physical 

In a similar way, we have 



formation for ai ^ 2 ■ Physically ai ^ 2 and ai ^ 2 can be two output fields of a beamsplitter. 



u a 1 + a 2u _ 1 = Ui a 1 + a lu _ 1 



V2 

,t , J 



V2 



U 



a i + a 2 U - 1 = u 1 ^^ur 1 = 



at + at 



V2 



V2 



Aio 2 
1 

Aio 2 



( CT 4 + 4) £i+3_ (fT . 1 
v2 

al + at 



; ,\ 4 + 4 



V2 

-(- 4 + 4)^ + (a 4 -4)^ 



(603) 



which for the quadrature = ( o$ + a\ ) / a/2, -Pi = \ a i — a \) /V%i, (* = 1,2), leads to 



c/ ^ i +^2 [/ _ 1 = ^ (Pi+P2)) [/(Pi + p 2)L/ -i = __L^i+^ 5 (604) 

from which we see that the Hadamard-Fresnel combinatorial operator also plays the role of exchang- 
ing the total momentum — average position followed by a squeezing transform, with the squeezing 
parameter being ^- . 

The mutual transform in ( 16041 ) can be realized by 

e <f(«I«i+«4»a) (Xi + X 2 ) e -if («I«i+4«2) = P x + p 2 , (605) 

e if (ala 1+ a+a 2 ) ^ + g -tf (a} 01 +a| a 2 ) _ _ ^ + ^ (ggg-) 



while the two-mode squeezing operator is 52 = cxp In ^ ( a|a 2 — a x a 2 



, therefore 



[/l = S^V* (607) 

From Eq. 05991 ) and Eq. ( 16071 ), we see that the Hadamard-Fresnel combinatorial operator can be 
decomposed as 

U = U 2 S 2 ; 1 e i %( a i ai+a z a2 ) = 5 2 " 1 e i *W Ql+a 2 a2 )[/ 2 . (608) 

It can be also seen that U is unitary, U + U = UU + = 1. 

In this section, we have introduced the Fresnel-Hadamard combinatorial operator by virtue of 
the IWOP technique. This unitary operator plays the role of both Fresnel transformation for mode 
ai ^ 2 and Hadamard transformation for mode ai ^ 2 , respectively, and the two transformations are 
combinatorial. We have shown that the two transformations are concisely expressed in the coherent- 
entangled state representation as a projective operator in integration form. We also found that the 
Fresnel-Hadamard operator can be decomposed as U 2 S 2 1 e l ^^ aiCL1+a2a2 \ a Fresnel operator U 2 , a 
two-mode squeezing operator S^ 1 and the total momentum-average position exchanging operator. 
Physically, ai ^ 2 and ai ^ 2 can be two output fields of a beamsplitter. If an optical device can be 
designed for Fresnel-Hadamard combinatorial transform, then it can be directly applied to these two 
output fields of the beamsplitter. 



In summary, although quantum optics and classical optics are so different, no matter in the math- 
ematical tools they employed or in a conceptual view (quantum optics concerning the wave-particle 
duality of optical field with an emphasis on its nonclassical properties, whereas classical optics works 
on the distribution aqnd propagation of the light waves), that it is a new exploration to link them 
systematically. However, In this review, via the route of developing Dirac's symbolic method we have 
revealed some links between them by mapping classical symplectic transformation in the coherent 
state representation onto quantum unitary operators (GFO), throughout our discussion the IWOP 
technique is indispensable for the derivation. We have resorted to the quantum optical interpreta- 
tion of various classical optical transformations by adopting quantum optics concepts such as the 
coherent states, squeezed states, and entangled states, etc. Remarkably, we have endowed com- 
plex fractional Fourier transform, Hankel transform with quantum optical representation-transform 
interpretation. Our formalism, starting from quantum optics theory, not only provides quantum me- 
chanical account of various classical optical transformations, but also have found their way back to 
some new classical transformations, e.g. entangled Fresnel transform, Fresnel-wavelet transform, 
etc, which may have realistic optical interpretation in the future. As Dirac predicted, functions that 
have been applied in classical optical problems may be translated in an operator language in quan- 
tum mechanics, and vice-versa. We expect that the content of this work may play some role in 
quantum states engineering, i.e., optical field states' preparation and design. 

Once the correspondence in this respect between the two distinct fields is established, the power 
of Dirac's symbolic method can be fully displayed to solve some new problems in classical optics, 



e.g., to find new eigen-modes of some optical transforms; to extend the research region of classical 
optics theoretically by introducing new transforms (for example, the entangled Fresnel transforms), 
which may bring attention of experimentalists, who may get new ideas to implement these new 
classical optical transformations. 
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